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Abstract 


This  final  report  describes  the  progress  to  date  on  six  individual  projects. 

1.  “The  Effect  of  Working  Fluid  Inventory  on  the  Performance  of  Revolving  Helically- 
Grooved  Heat  Pipes”  is  related  to  an  experimental  analysis  of  working  fluid  inventory 
and  its  effects  on  the  capillary  limit  of  grooved  heat  pipes. 

2.  “Fully-Developed  Laminar  Flow  in  TVapezoidal  Grooves  with  Shear  Stress  at  the 
Li  quid- Vapor  Interface”  concerns  the  effects  of  geometry  and  interfacial  shear  stress 
on  the  pressure  drop  in  trapezoidal  grooves. 

3.  “Fully-Developed  Laminar  Flow  in  Sinusoidal  Grooves”  is  similar  to  the  above,  except 
that  the  sinusoidal  grooves  typically  seen  in  micro-etching  are  examined. 

4.  “Micro  Capillary  Pumped  Loop  Testing”  deals  with  the  design,  manufacture  and  cal¬ 
ibration  of  an  experimental  facility  for  testing  capillary  pumped  loops  embedded  in 
silicon  wafers.  A  description  of  the  current  experimental  setup  and  the  details  of  the 
cahbration  of  the  laser  controller  and  the  calibration  of  the  thermocouples  are  provided. 
The  second  project  is  entitled, 

5.  “Design  and  Testing  of  a  Thermodynamic  Filling  Station  for  Miniature  Heat  Pipes,” 
which  concerns  a  novel  filhng  procedime  for  miniature  heat  pipes.  The  experimental 
setup  and  procedure  is  described,  which  involves  a  vacuum  chamber  within  an  environ¬ 
mental  chamber.  A  thermodynamic  analysis  of  the  process  is  provided,  and  the  results 
are  compared  to  the  experimental  data.  The  final  project  of  this  report  is  entitled, 

6.  “Acceleration  Testing  of  Three  Raytheon  Heat  Pipes  using  the  AFRL/PRPS  Cen¬ 
trifuge  Table”.  The  objective  of  this  project  was  to  determine  the  performance  of  heat 
pipes  under  various  acceleration  conditions,  which  were  imposed  by  the  centrifuge  ta¬ 
ble  located  in  Bldg.  71-B  H-Bay  in  Area  B  of  Wright- Patterson  AFB.  Information 
on  the  experimental  setup,  experimental  procedure,  and  an  analysis  of  the  results  axe 
provided. 
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1  The  Effect  of  Working  Fluid  Inventory  on  the  Performance  of  Revolving 

Helically-Grooved  Heat  Pipes 

1.1  Abstract 

The  results  of  a  recently  completed  experimental  and  analytical  study  showed  that  the 
capillary  limit  of  a  helically-grooved  heat  pipe  (HGHP)  was  increased  significantly  when  the 
transverse  body  force  field  was  increased.  This  was  due  to  the  geometry  of  the  helical  groove 
wick  structure.  The  objective  of  the  present  research  was  to  experimentally  determine  the 
performance  of  revolving  helically-grooved  heat  pipes  when  the  working  fluid  inventory  was 
varied.  This  report  describes  the  measurement  of  the  geometry  of  the  heat  pipe  wick  struc¬ 
ture  and  the  construction  and  testing  of  a  heat  pipe  filling  station.  In  addition,  an  extensive 
analysis  of  the  uncertainty  involved  in  the  filling  procedure  and  working  fluid  inventory  has 
been  outlined.  Experimental  measurements  include  the  maximum  heat  transport,  thermal 
resistance  and  evaporative  heat  transfer  coefficient  of  the  revolving  helically-grooved  heat 
pipe  for  radial  accelerations  of  |or|  =  0.0,  2.0,  4.0,  6.0,  8.0,  and  10.0-g  and  working  fluid  fills 
of  G  =  0.5,  1.0  and  1.5.  An  existing  capillary  limit  model  was  updated  and  comparisons 
were  made  to  the  present  experimental  data. 

1.2  Introduction 

Helically-grooved  heat  pipes  (HGHPs)  have  potential  applications  in  the  thermal  man¬ 
agement  of  rotating  equipment  such  as  aircraft  alternators,  large-scale  industrial  electric 
motors,  and  spinning  satellites.  In  two  recent  studies  (Klasing  et  al.,  1999;  Thomas  et  ah, 
1998),  the  performance  of  revolving  HGHPs  was  investigated.  It  was  found  that  the  capil¬ 
lary  limit  increased  with  the  strength  of  the  acceleration  field  perpendicular  to  the  heat  pipe 
axis.  In  order  to  move  HGHPs  closer  to  application,  knowledge  must  be  gained  concern¬ 
ing  the  sensitivity  of  the  capillary  limit  to  working  fluid  fill  amount,  since  variations  in  the 
fill  amount  are  inevitable  during  the  manufacture  of  these  devices.  Very  few  studies  were 
available  concerning  the  effect  of  working  fluid  fill  on  the  performance  of  axially-grooved 
heat  pipes,  but  those  found  have  been  outlined  below.  In  addition,  synopses  of  the  two 
aforementioned  studies  on  revolving  HGHPs  have  also  been  provided. 
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Brennan  et  al.  (1977)  developed  a  mathematical  model  to  determine  the  performance 
of  an  axially-grooved  heat  pipe  which  accounts  for  liquid  recession,  liquid-vapor  shear  in¬ 
teraction  and  puddle  flow  in  a  1-g  acceleration  environment.  The  model  considered  three 
distinct  flow  zones:  the  grooves  unaffected  by  the  puddle,  the  grooves  that  emerge  from 
the  puddle,  and  the  grooves  that  are  submerged  by  the  puddle.  The  model  for  the  puddle 
consisted  of  satisfying  the  equation  of  motion  for  the  puddle  and  the  continuity  equation 
at  the  puddle-groove  interface,  and  was  solved  by  a  fourth-order  Runge-Kutta  integration 
method  with  self-adjusting  step  sizes.  The  assumptions  made  by  the  model  for  the  puddle 
were  uniform  heat  addition  and  removal  with  a  single  evaporator  and  a  single  condenser 
section,  and  one-dimensional  laminar  flow  in  the  puddle.  The  transport  capability  of  the 
grooves  unaffected  by  the  puddle  and  the  grooves  extending  beyond  the  puddle  were  approx¬ 
imated  by  a  closed-form  solution  with  laminar  liquid  and  vapor  flow.  The  working  fluids 
used  for  the  experiment  were  methane,  ethane  and  ammonia.  Brennan  et  al.  (1977)  stated 
that  the  mathematical  model  agreed  well  with  the  experimental  data  for  ideally  filled  and 
overfilled  heat  pipes,  but  some  differences  were  noted  for  underfilled  heat  pipes.  In  general, 
it  was  found  for  ideally  filled  heat  pipes  the  predicted  transported  heat  was  higher  than  that 
measured.  Also,  this  discrepancy  was  more  significant  for  lower  operating  temperatures.  In 
addition,  it  was  found  diming  the  experiments  that  the  maximum  transported  heat  increased 
with  fill  volume. 

Vasihev  et  al.  (1981)  performed  a  series  of  experiments  on  an  aluminum  axially-grooved 
heat  pipe  which  was  overfilled  and  ideally  filled.  The  width  and  height  of  the  grooves  were 
w  =  0.123  mm  and  h  =  0.7  mm,  respectively,  with  an  overall  heat  pipe  length  of  Lt  — 
80.0  cm.  The  working  fluids  were  acetone  and  ammonia.  Vasiliev  et  al.  showed  that  the 
temperature  difference  from  the  evaporator  to  the  adiabatic  regions  increased  at  a  much 
slower  rate  with  increasing  overfills.  This  was  attributed  to  a  thin  film  of  Liquid  emerging 
from  the  overfill  pool  wetting  the  upper  grooves.  Vasiliev  et  al.  stated  that  this  thin  film 
was  hfted  over  the  grooves  by  capillary  forces  due  to  microroughness  on  the  groove  surface. 
A  mathematical  model  was  developed  for  low  temperature  axially-grooved  heat  pipes  to 
estimate  heat  pipe  performance  for  0-g  and  1-g  applications.  The  mathematical  model  was 
a  set  of  boundary-value  problems  applied  to  each  groove  and  was  solved  by  a  numerical 
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iteration  method.  The  model  was  based  on  pressure  balance  equations  and  mass  continuity 
written  for  a  single  groove.  The  temperature  of  the  vapor  in  the  adiabatic  region  was  an 
input  parameter,  and  the  vapor  pressure  gradient  was  assumed  to  be  one-dimensional.  In 
addition,  the  liquid- vapor  shear  stress  was  assumed  to  be  constant,  and  the  starting  liquid 
film  thickness  was  of  the  same  order  of  magnitude  as  the  groove  microroughness.  Very  good 
agreement  was  reported  between  the  mathematical  model  and  experimental  transported  heat 
results  for  ideally  filled  and  overfilled  heat  pipes  under  gravity. 

Thomas  et  al.  (1998)  presented  experimental  data  obtained  from  a  helically-grooved 
copper  heat  pipe  which  was  tested  on  a  centrifuge  table.  The  heat  pipe  was  bent  to  match 
the  radius  of  curvature  of  the  table  so  that  uniform  transverse  (perpendicular  to  the  axis  of 
the  heat  pipe)  body  forces  field  could  be  applied  along  the  entire  length  of  the  pipe.  The 
steady-state  performance  of  the  curved  heat  pipe  was  determined  by  varying  the  heat  input 
{Qin  —  25  to  250  W)  and  centrifuge  table  velocity  (radial  acceleration  |ar|  =  0.01  to  10-g).  It 
was  found  that  the  capillary  limit  increased  by  a  factor  of  five  when  the  radial  acceleration 
increased  from  |ar|  =  0.01  to  6-g  due  to  the  geometry  of  the  helical  grooves.  A  model  was 
developed  to  calculate  the  capillary  limit  of  each  groove  in  terms  of  centrifuge  table  angular 
velocity,  the  geometry  of  the  heat  pipe  and  the  grooves,  and  the  temperature-dependent 
working  fluid  properties.  The  agreement  between  the  model  and  the  experimental  data  was 
satisfactory. 

Klasing  et  al.  (1999)  developed  a  mathematical  model  to  determine  the  operating  limits 
of  a  revolving  hehcally-grooved  straight  heat  pipe.  The  capillary  limit  calculation  required  an 
analysis  of  the  total  body  force  imposed  by  rotation  and  gravity  on  the  liquid  along  the  length 
of  the  helical  grooves.  The  boiling  and  entrainment  limits  were  calculated  using  methods 
described  by  Faghri  (1995).  It  was  found  that  the  capillary  limit  increased  significantly  with 
rotational  speed  due  to  the  helical  geometry  of  the  heat  pipe  wick  structure.  The  maximum 
heat  transport  was  found  to  be  a  function  of  angular  velocity  and  tilt  angle  from  horizontal. 
In  addition,  a  minimum  value  of  angular  velocity  was  required  to  obtain  the  benefits  of  the 
hehcal  groove  geometry. 

The  first  objective  of  the  present  study  was  to  determine  the  sensitivity  of  the  performance 
of  revolving  HGHPs  to  the  working  fluid  fill  amount.  This  required  a  precise  knowledge  of 
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the  geometry  of  the  heat  pipe  and  helical  grooves.  In  addition,  a  precision  filhng  station  was 
constructed  and  calibrated  to  determine  the  uncertainties  involved  in  the  filling  procedure. 
The  copper-ethanol  heat  pipe  was  tested  on  a  centrifuge  table  at  Wright-Patterson  AFB 
(AFRL/PRPG)  to  determine  the  capillary  limit,  thermal  resistance  and  evaporative  heat 
transfer  coefficient  for  fiU  ratios  of  G  =  0.5,  1.0  and  1.5,  and  radial  accelerations  of  |ar|  = 
0.01,  2.0,  4.0,  6.0,  8.0  and  10.0-g.  The  second  objective  of  the  present  study  was  to  improve 
the  existing  analytical  capillary  limit  model  developed  by  Thomas  et  al.  (1998)  using  the 
above-mentioned  geometric  measurements  and  by  using  improved  equations  for  the  working 
fluid  properties. 

1.3  Determination  of  Heat  Pipe  Working  Fluid  Inventory 


The  objective  of  this  analysis  was  to  determine  the  working  fluid  inventory  of  a  HGHP, 
which  consists  of  the  mass  of  liquid  in  the  grooves  and  the  mass  of  vapor  in  the  vapor 
space.  Since  the  heat  pipe  is  a  closed  container  under  saturation  conditions,  the  total  mass 
of  working  fluid  in  the  heat  pipe  is  given  by 

mt  =  rnv+‘mi  = - 1- — —  (1) 

Vv  vi 

where  G  =  Vi/Vgr  is  the  ratio  of  the  volume  of  Hquid  to  total  groove  volume.  The  volume  of 
the  vapor  space  is 

Ks  =  +  Vgr{l  -  G)  (2) 

The  second  term  in  eqn.  (2)  accounts  for  the  increase  or  decrease  in  the  vapor  space  volume 
when  the  parameter  G  is  varied.  The  total  volume  of  the  grooves  is 


kgr  —  Lg[Ngj-Ag[ 


(3) 


A  cross-sectional  view  of  a  typical  helical  groove  in  the  experimental  test  article  is  shown 
in  Fig.  1 .  The  cross-sectional  area  of  the  trapezoidal  groove  accounts  for  the  differing  side 
angles. 


^gr 


wh  +  -/i^(tan  9i  +  tan  62) 


The  total  length  of  each  groove  is 
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The  radius  of  the  helix  is  given  by 


2  (■^vs  T 


(6) 


The  helical  pitch  is  the  distance  through  which  the  hehx  makes  one  revolution  around  its 
radius. 


P 


2x(s  —  Si) 


(7) 


{(f)  -  <j){) 

The  helix  angle  (f)  corresponds  to  s,  which  is  the  distance  traveled  along  the  centerHne  of  the 
heat  pipe. 

In  order  to  calculate  the  working  fluid  inventory  for  the  HGHP,  measurements  of  the 
appropriate  geometric  parameters  were  made.  In  addition,  an  extensive  uncertainty  analysis 
was  performed  to  determine  the  uncertainties  of  both  the  measured  and  calculated  variables 
used  in  finding  the  working  fluid  inventory. 

The  physical  variables  given  in  eqn.  (4)  for  the  cross-sectional  area  of  the  grooves  have 
been  measured.  A  sample  of  the  HGHP  container  was  set  in  an  epoxy  resin  mold,  polished, 
and  examined  under  a  microscope  with  50x  magnification.  Computer  software  was  used  to 
make  bitmap  pictures  of  ten  different  grooves  and  a  microscopic  caHbration  scale.  These 
pictures  were  then  analyzed  to  determine  the  geometric  values  shown  in  Fig.  1.  Since  the 
corners  at  the  top  of  the  land  between  grooves  were  not  weU  defined,  a  special  procedure 
was  established  to  determine  the  geometry  of  the  grooves.  First,  lines  were  drawn  along 
the  bottom  and  sides  of  each  groove.  Then,  a  Hne  was  drawn  across  the  bottom  of  the 
land  between  grooves,  as  shown  in  Fig.  1.  This  hne  was  then  transposed  to  the  top  of  the 
land.  The  intersections  between  this  line  and  the  lines  along  the  sides  of  the  groove  were 
defined  as  the  upper  corners  of  the  groove.  Note  that  the  two  fines  along  the  land  tops  are 
at  different  angles  due  to  the  radius  of  curvature  of  the  heat  pipe  container.  The  angles  6i 
and  02,  and  the  height  and  width  of  the  groove  h  and  w  were  found  using  a  bitmap  picture 
of  the  microscopic  calibration  scale  as  described  by  Castle  (1999). 

An  optical  comparator  was  used  to  determine  the  vapor  space  diameter  of  the  heat  pipe 
container  sample.  The  cross  hairs  of  the  optical  comparator  were  carefully  aligned  with  the 
top  of  the  land  between  grooves  on  the  left  edge  of  the  pipe.  The  comparator  table  was  then 
moved  until  the  land  tops  on  the  right  edge  of  the  pipe  were  aligned  with  the  cross  hairs. 
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The  diameter  of  the  heat  pipe  vapor  space  was  the  distance  of  the  table  movement. 

The  hehcal  groove  pitch  was  found  using  a  vertical  milling  machine  and  an  angular 
displacement  transducer.  The  heat  pipe  container  material  was  originally  1  m  long.  Approx¬ 
imately  one-half  was  used  to  form  the  heat  pipe,  and  the  other  half  was  used  to  determine 
the  pitch.  The  rotation  angle  (0  —  ^i)  and  the  corresponding  distance  along  the  centerline 
of  the  heat  pipe  (s  —  Si)  has  been  found  as  shown  in  Fig.  2(a).  A  heat  pipe  holding  de¬ 
vice  was  constructed  from  two  angle  aluminum  uprights  mounted  to  the  table  of  a  vertical 
milhng  machine.  Precision  alignment  blocks  were  attached  to  the  undersides  of  the  uprights 
to  engage  one  of  the  grooves  in  the  milling  machine  table  for  improved  alignment.  Nylon 
bushings  were  placed  in  the  uprights  to  center  both  the  heat  pipe  container  and  the  shaft, 
which  was  concentric  with  the  heat  pipe  container.  A  small  pin  was  made  from  a  1.58  mm 
(0.0625  in)  dowel  pin,  where  one  end  was  ground  to  0.26  mm  to  fit  in  the  base  of  the  helical 
groove.  This  sprung  pin  was  set  in  a  hole  in  the  shaft  where  it  engaged  one  of  the  grooves,  as 
shown  in  Fig.  2(6).  An  angular  displacement  transducer  was  mounted  onto  another  piece  of 
angle  aluminum.  A  vertical  6.35  mm  (0.25  in)  dowel  pin  was  placed  in  the  angle  aluminum 
to  align  with  the  angular  displacement  transducer  shaft.  The  dowel  pin  was  held  by  a  col¬ 
let  installed  in  the  milling  machine  spindle  in  order  to  fix  the  location  of  the  displacement 
transducer.  The  shaft  of  the  transducer  was  linked  to  the  shaft  within  the  heat  pipe  by 
three  set  screws.  As  the  milhng  machine  table  moved  the  pipe  over  the  stationary  shaft, 
the  pin  followed  the  hehcal  groove,  causing  the  shaft  to  rotate.  The  angular  displacement 
transducer  measured  this  rotation.  A  multimeter  was  used  to  measure  the  output  voltage 
of  the  angular  displacement  transducer.  The  distance  of  the  table  movement  was  (s  —  Si), 
which  was  read  from  the  milhng  machine  display  unit.  The  transducer  output  voltage  was 
measured  over  10  cm  lengths  for  ten  different  groups.  Backlash  errors  were  avoided  by  not 
reversing  the  table  movement  while  taking  data.  The  pitch  was  calculated  using  eqn.  (7)  at 
a  point  in  the  center  of  each  10  cm  length.  An  average  of  88  values  were  used  to  calculate 
the  hehcal  pitch. 

Using  the  analysis  given  by  Miller  (1989),  the  root-sum-square  uncertainties  for  the 
groove  cross-sectional  area,  hehcal  pitch,  hehx  radius,  groove  length,  groove  volume,  vapor 
space  volume,  and  total  mass  of  the  working  fluid  inventory  have  been  calculated.  The 
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Table  1:  The  geometric  variable  values  associated  with  the  working  fluid  inventory. 


Measured  Values 

h 

0.03831  ±  0.00076  cm 

w 

0.03445  ±  0.0010  cm 

0i 

15.44°  ±  0.91° 

l^l 

13.80°  ±  0.96° 

1.359  ±  0.005  cm 

m 

43.8  ±  0.084  cm 

Calculated  Values 

E9 

1.703  X  10-=*  ±  6.0  X  10-^  cm^" 

p 

135.8  ±  5.9  cm 

Th 

0.6992  ±  0.0025  cm 

Lgr 

43.82  ±  0.84  cm 

3.73  ±  0.13  cm=* 

measured  and  calculated  imcertainties  for  all  geometric  variables  presented  are  shown  in 
Table  1. 

A  hterature  survey  was  completed  to  determine  the  specific  volumes  of  ethanol  vapor  and 
hquid  at  various  saturation  temperatures,  as  shown  in  Fig.  3.  This  information  was  needed 
to  determine  the  total  mass  and  xmcertainty  of  the  working  fluid  inventory  rrit  ±  ArUf.  While 
existing  texts  report  these  properties  (Faghri,  1995;  Peterson,  1994,  Lide  and  Kehiaian,  1994; 
Carey,  1992;  Schlunder,  1983;  Ivanovskii  et  al.,  1982),  it  was  found  that  most  simply  referred 
to  previous  sources.  Therefore,  the  data  shown  in  Fig.  3  represent  information  gathered  from 
primary  sources  that  cannot  readily  be  traced  further.  In  Fig.  3(o),  the  available  data  for  the 
specific  volume  of  liquid  in  the  range  of  Tgat  =  0  to  100°C  are  relatively  scattered.  Vargaftik 
(1975)  stated  that  the  ethanol  used  was  96%  pure  by  volume,  with  water  making  up  most 
of  the  other  4%.  Ethanol  is  agressively  hygroscopic,  so  special  procedures  are  required  for 
further  purification  as  outlined  by  Timmermans  (1950)  concerning  anhydrous  ethanol.  Since 
the  data  by  Timmermans  (1950)  and  TRC  (1983)  are  nearly  coincident,  it  is  beheved  that 
the  data  reported  by  TRC  (1983)  are  also  for  anhydrous  ethanol.  Dunn  and  Reay  (1978) 
do  not  provide  information  concerning  purity.  Therefore,  the  Vargaftik  (1975)  data  and  the 
Dunn  and  Reay  (1978)  data  have  been  discarded  in  Fig.  3(o).  In  Fig.  3(6),  the  deviation 
of  the  Dunn  and  Reay  (1978)  data  for  the  specific  volume  of  vapor  is  significant.  Therefore, 
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Table  2:  The  calculated  total  mass  of  the  working  fluid  inventory. 


G 

mt  (g) 

Amt  (g) 
Calculated 

Amd  (g) 
Filling 
Station 

0.5 

1.47 

±  3.6% 

±  5.0% 

1.0 

2.92 

±  3.7% 

±  2.9% 

1.5 

4.38 

±  3.6% 

±  1.9% 

the  Dunn  and  Reay  (1978)  data  has  been  discarded  in  Fig.  3(6).  Polynomial  curve  fits  from 
0  <  Tsat  <  100°C  have  been  obtained  for  the  data  shown  in  Figs.  3(a)  and  3(6)  for  the 
specific  volumes  of  liquid  and  vapor  ethanol.  These  curve  fits  have  been  evaluated  at  room 
temperature  to  determine  the  proper  values  to  be  used  in  the  uncertainty  analysis,  since  the 
heat  pipe  was  filled  at  room  temperature.  Information  concerning  the  uncertainty  of  the 
original  data  was  not  available.  Therefore,  the  uncertainties  of  these  properties  have  been 
estimated  to  be  the  maximum  variance  of  the  data  from  the  curve  fits  (Avi  =  3.5  x  lO'® 
m^/kg,  Avv  =  0.39  m^/kg).  The  specific  volumes  of  liquid  and  vapor  ethanol  (m^/kg)  as 
functions  of  saturation  temperature  (°C)  are  shown  below  for  the  range  0  <  Tgat  <  100°C 

vi  =  exp(ao  +  OiTgat  +  +  a4rg'*at)/1000  (8) 

v-u  —  exp(6o  +  6iTsat  +  62T^t  +  +  ^4T^t)/1000  (9) 


where  the  coefficients  are 


Oo  =  0.2153  6o  =  10.35 


ai  =  1.049  X  10-^ 

02  =  -1.345  X  10“® 
03  =  2.025  X  10-® 

04  =  -5.474  X  10-11 


=  -6.375  X  10-2 

62  =  1.735  X  10-^ 

63  =  5.714  X  10-'^ 

64  =  -6.003  X  10-® 


The  total  mass  of  the  working  fluid  inventory  nit  and  the  associated  uncertainty  Am^  for 
the  range  of  fill  values  are  given  in  Table  2. 
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1.4  Heat  Pipe  Filling  Station 


A  filling  station  has  been  constructed  which  is  capable  of  placing  a  low-temperature 
working  fluid  (i.e.,  water,  ethanol,  methanol)  into  a  heat  pipe  without  also  introducing 
ambient  air  (Fig.  4).  The  station  consisted  of  a  manifold  of  valves  and  interconnecting 
stainless  steel  tubing,  a  working  fluid  reservoir,  a  dispensing  burette,  a  vacuum  pump,  and 
a  container  of  compressed  dry  nitrogen  gas.  Previous  experience  with  filling  stations  showed 
that  long  runs  of  horizontal  tubing  could  cause  significant  filling  errors  due  to  vapor  bubbles 
within  the  tubing.  To  address  this  problem,  the  manifold  was  constructed  such  that  the 
interconnecting  tubing  runs  were  very  short  (on  the  order  of  2  cm).  In  addition,  the  tubes 
which  intersect  the  main  vertical  tube  between  valves  2  and  5  (Fig.  4)  were  offset  from  each 
other  and  ran  at  a  diagonal  from  the  main  tube.  Again,  the  purpose  of  this  design  was  to 
reduce  the  possibility  of  vapor  bubbles  adhering  to  the  tubing  walls,  thus  causing  errors  in 
the  fill  amount.  However,  it  is  likely  that  some  vapor  still  does  adhere  to  the  tubing,  so 
certain  procedures  were  carried  out  during  fiUing  to  eject  as  much  vapor  as  possible.  For 
instance,  the  1  psig  relief  valve  over  valve  1  was  cycled  on  and  off  several  times.  In  addition, 
valves  2  and  5  were  cycled  on  and  off  while  noting  the  meniscus  displacement  within  the 
dispensing  burette.  If  the  meniscus  was  displaced  more  than  0.06  cm^,  vapor  was  probably 
trapped  within  the  valve.  The  valve  in  question  was  then  cycled  until  the  bubble  was  ejected. 

To  fill  the  heat  pipe,  the  container  was  first  evacuated  to  a  pressure  of  10~®  Torr  using 
a  turbomolecular  vacuum  pump.  The  sealed  pipe  was  then  connected  to  the  filling  station 
at  valve  5.  The  working  fluid  was  frozen  and  thawed  repeatedly  to  reduce  the  amount  of 
dissolved  air  within  the  fluid.  The  entire  filling  station  was  then  evacuated  by  a  roughing 
pump,  except  the  working  fluid  reservoir.  After  evacuation,  the  liquid  working  fluid  was 
drawn  up  into  the  dispensing  burette  and  into  all  interconnecting  tubing.  After  noting  the 
height  of  the  meniscus,  the  desired  amount  of  working  fluid  was  metered  into  the  heat  pipe 
by  carefully  opening  the  heat  pipe  fill  valve  8.  The  difference  in  height  of  the  hquid  column 
was  related  to  the  dispensed  mass  of  working  fluid. 

During  initial  testing  of  the  filling  station,  it  was  found  that  the  mass  of  working  fluid 
dispensed  into  the  heat  pipe  container  was  different  than  what  was  indicated  by  the  dis- 
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Table  3:  HeUcally-grooved  heat  pipe  specifications. 


Working  fluid 

Ethanol 

Working  fluid  charge 

rrit  =  1.47,  2.92  and  4.38  g 

Evaporator  length 

Le  =  15.2  ±  0.16  cm 

Adiabatic  length 

La  =  8.2  ±  0.16  cm 

Condenser  length 

La  =  15.2  ±  0.16  cm 

Tube  outside  diameter 

Do  =  1.588  ±  0.005  cm 

Tube  wall  thickness 

ty,  =  0.0757  cm 

Radius  of  curvature 

R  =  1.22  m 

Wall/wick  materials 

Copper 

Wick  structure 

Helical  grooves 

Number  of  Grooves 

-  50 

Heater  element 

Nichrome  heater  tape 

Fill  valve 

Nupro  B-4HW  bellows  valve 

Calorimeter 

1/8  in.  OD  coiled  copper  tubing 

pensing  burette.  Therefore,  a  rigorous  calibration  of  the  filling  station  was  undertaken  to 
determine  a  correlation  between  the  change  in  volume  read  by  the  dispensing  burette  and 
the  change  in  mass  of  a  receiving  burette  attached  at  valve  5,  which  was  measured  using  a 
precision  scale.  The  total  uncertainty  of  the  working  fluid  inventory  dispensed  by  the  heat 
pipe  filhng  station  l^rrid  is  given  in  Table  2. 

1.5  Experimental  Setup 

The  purpose  of  the  experiment  was  to  examine  the  steady-state  performance  of  a  helically- 
grooved  copper-ethanol  heat  pipe  under  various  heat  inputs  and  transverse  body  force  fields 
using  a  centrifuge  table  located  at  Wright-Patterson  AFB  (AFRL/PRPG).  Specifically,  the 
amount  of  working  fluid  was  varied  (G  =  0.5,  1.0  and  1.5)  to  determine  the  effects  of  un¬ 
der/overfilling  on  the  capillary  limit,  thermal  resistance  and  evaporative  heat  transfer  coef¬ 
ficient  of  the  HGHP.  To  ensure  uniform  radial  acceleration  fields  over  the  length  of  the  heat 
pipe,  the  pipe  was  bent  to  match  the  radius  of  curvature  of  the  centrifuge  table  {R  =  1.22  m). 
Physical  information  concerning  the  heat  pipe  is  given  in  Table  3.  It  should  be  noted  that 
the  total  helix  angle  was  very  small:  Each  groove  rotated  through  an  angle  of  approximately 
2.03  rad  (116  arc  degrees)  over  the  length  of  the  pipe.  The  heat  pipe  was  mounted  to  a 
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platform  overhanging  the  edge  of  the  horizontal  centrifuge  table.  This  allowed  the  heat  pipe 
to  be  positioned  such  that  the  radius  of  curvature  was  equivalent  to  the  outermost  radius 
of  the  centrifuge  table.  Insulative  mounting  blocks  were  used  to  ensure  that  the  heat  pipe 
matched  the  prescribed  radius  as  closely  as  possible.  The  horizontal  centrifuge  table  was 
driven  by  a  20-hp  dc  motor.  The  acceleration  field  near  the  heat  pipe  was  measured  by  a 
triaxial  accelerometer.  The  acceleration  field  at  the  centerline  of  the  heat  pipe  radius  was 
calculated  from  these  readings  using  a  coordinate  transformation. 

A  pressure-sensitive  nichrome  heater  tape  with  an  aluminized  backing  was  uniformly 
wound  around  the  circumference  of  the  evaporator  section  for  heat  input.  Power  was  supplied 
to  the  heat  pipe  evaporator  section  by  a  power  supply  through  power  slip  rings  to  the  table. 
While  the  current  reading  could  be  made  directly  using  a  precision  ammeter,  the  voltage 
across  the  electric  heater  had  to  be  measured  on  the  rotating  table  because  of  the  voltage  drop 
between  the  control  room  and  the  table.  Therefore,  the  voltage  at  the  heater  was  obtained 
through  the  instrumentation  slip  ring  assembly  and  read  by  a  precision  multimeter. 

The  calorimeter  consisted  of  a  length  of  1/8  in.  OD  copper  tubing  wound  tightly  around 
the  condenser  section.  The  size  of  the  tubing  was  chosen  to  be  small  to  minimize  the  effects 
of  acceleration  on  the  performance  of  the  calorimeter.  Thermal  grease  was  used  between  the 
heat  pipe  and  the  calorimeter  to  decrease  contact  resistance.  Type  T  thermocouples  were 
inserted  through  brass  T-branch  coimectors  into  the  coolant  inlet  and  exit  streams,  and  a 
high-resolution  digital  flow  meter  was  used  to  measure  the  mass  flow  rate  of  the  coolant 
(50%  by  mass  ethylene  glycol/water  mixture).  The  mass  flow  rate  was  controlled  using  a 
high-pressure  booster  pump,  which  aided  the  low-pressure  pump  in  the  recirculating  chiller. 
The  percentage  of  ethylene  glycol  was  measured  periodically  during  testing  using  a  precision 
hydrometer  to  ensure  that  the  mixture  did  not  change.  The  temperature  of  the  coolant 
was  maintained  at  a  constant  setting  by  the  recirculating  chiller.  Coolant  was  delivered  to 
the  centrifuge  table  via  a  double-pass  hydraulic  rotary  coupling.  The  mass  flow  rate  was 
constant  for  all  experiments.  Values  of  the  specific  heat  of  ethylene  glycol/water  mixtures 
were  obtained  from  ASHRAE  (1977),  which  were  in  terms  of  percent  ethylene  glycol  by 
weight  and  temperature.  The  average  temperature  between  the  calorimeter  inlet  and  outlet 
was  used  to  evaluate  the  specific  heat.  The  specific  heat  did  not  vary  appreciably  during 
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testing  since  it  is  a  weak  function  of  temperature. 

Heat  pipe  temperatures  were  measured  by  Type  T  surface-mount  thermocouples,  which 
were  held  in  place  using  Kapton  tape.  Mounting  locations  for  the  thermocouples  are  shown 
in  Fig.  5.  A  short  unheated  length  next  to  the  evaporator  end  cap  was  instrumented  with 
thermocouples  specifically  for  accurate  thermal  resistance  measurements.  In  addition,  groups 
of  four  thermocouples  were  arranged  around  the  circumference  of  the  heat  pipe  at  stations 
in  the  evaporator  section  for  local  heat  transfer  coefficient  information.  Temperature  signals 
were  conditioned  and  amplified  on  the  centrifuge  table.  These  signals  were  transferred  off 
the  table  through  the  instrumentation  slip  ring  assembly,  which  was  completely  separate 
from  the  power  sfip  ring  assembly  to  reduce  electronic  noise.  Conditioning  the  temperature 
signals  prior  to  leaving  the  centrifuge  table  eliminated  difficulties  associated  with  creating 
additional  junctions  within  the  shp  ring  assembly.  Temperature  and  acceleration  signals 
were  collected  using  a  personal  computer  with  data  logging  software.  Since  a  shortage  of 
thermocouple  channels  existed  on  the  centrifuge  table,  a  series  of  three  electrical  relays  were 
engaged  to  read  one  set  of  thermocouples,  and  disengaged  to  read  the  other  set. 

Since  the  heat  pipe  assembly  was  subjected  to  air  velocities  due  to  the  rotation  of  the 
table  (up  to  11  m/s  =  25  mi/hr),  efforts  were  made  to  reduce  convective  heat  losses  from 
the  exterior  of  the  heat  pipe.  A  thin-walled  aluminum  box  was  fabricated  to  fit  around  the 
heat  pipe.  Ceramic  wool  insulation  was  placed  inside  the  box  and  around  the  heat  pipe 
through  three  small  doors  on  the  top  of  the  box.  This  insulation/box  arrangement  provided 
an  effective  barrier  to  convective  losses  from  the  heat  pipe  to  the  ambient. 

The  helicaUy-grooved  copper-ethanol  heat  pipe  was  tested  in  the  following  manner.  The 
recirculating  chiller  was  turned  on  and  allowed  to  reach  the  setpoint  temperature,  which  was 
measured  at  the  calorimeter  inlet.  The  centrifuge  table  was  started  from  the  remote  control 
room  at  a  slow  constant  rotational  speed  to  prevent  damage  to  the  power  and  instrumentation 
slip  rings.  In  this  case,  the  radial  acceleration  was  less  than  |or|  <  0.01-g.  In  all  cases,  the 
centrifuge  table  rotated  in  a  clockwise  direction  as  seen  from  above.  Power  to  the  heater 
was  applied  (Qjn  =  10  W)  and  the  heat  pipe  was  allowed  to  reach  a  steady-state  condition. 
The  power  to  the  heater  was  then  increased  to  Qm  =  20  W  and  again  the  heat  pipe  was 
allowed  to  reach  a  steady-state  condition.  This  was  repeated  until  the  maximum  allowable 
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Table  4:  Maximum  uncertainties  of  measured  and  calculated  values. 


Measured  Values 

Coolant  mass  flow  rate 

Arfic  =  ±  0.05  g/s 

Heater  voltage 

AV  =  ±  0.5  V 

Heater  current 

A/  =  ±  0.1  A 

Radial  acceleration 

AOr  =  ih  0.1-g 

Calorimeter  inlet  temperature 

ATin  =  ±  0.07  K 

Calorimeter  outlet  temperature 

ATout  =  ±  0.08  K 

Evaporator  end  cap  temperature 

ATeec  =  ±  0.09  K 

Condenser  end  cap  temperature 

ATcec  =  ±  0.11  K 

Calculated  Values 

Heat  input 

See  Fig.  6 

Heat  transported 

AQt  =  ±  3.2  W 

Thermal  resistance 

See  Fig.  7 

Heat  transfer  coefficient 

See  Figs.  12  and  13 

evaporator  temperature  was  reached  (T-u,, max  =  100°C).  After  all  data  had  been  recorded  the 
power  to  the  heater  was  turned  off,  and  the  heat  pipe  was  allowed  to  cool  before  shutting 
down  the  centrifuge  table. 

Using  the  analysis  given  by  Miller  (1989),  the  uncertainties  for  aU  of  the  measured  and 
calculated  values  for  the  experimental  data  are  presented  in  Table  4. 

1.6  Results  and  Discussion 

The  objective  of  this  experiment  was  to  determine  the  steady-state  performance  of  a 
revolving  helically-grooved  heat  pipe  as  a  function  of  the  working  fluid  inventory.  The  heat 
input,  radial  acceleration  and  working  fluid  fill  were  varied  as  follows:  Qjn  =  10  to  180  W, 
|dr|  =  0.01  to  10-g,  and  G  =  0.5,  1.0  and  1.5.  Thermocouples  on  the  inboard,  outboard, 
top,  and  bottom  sides  of  the  heat  pipe  (Fig.  5)  were  used  to  determine  the  axial  and 
circumferential  temperature  distributions.  Typical  steady-state  temperature  distributions 
for  the  heat  pipe  for  G  =  1.0  at  |ar|  =  0.01-g  are  shown  in  Fig.  6.  For  low  power  input  levels, 
the  temperature  distribution  was  imiform.  As  the  power  input  increased,  the  temperatures 
within  the  evaporator  and  the  short  unheated  section  adjacent  to  the  evaporator  increased 
significantly,  indicating  a  partial  dryout  situation.  Since  the  coolant  temperature  and  flow 
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rate  were  constant  for  all  tests,  the  adiabatic  and  condenser  temperatures  increased  slightly 
with  input  power.  Figure  7  shows  the  thermal  resistance  versus  transported  heat  over  the 
entire  range  of  radial  acceleration  for  each  fill  level.  In  Fig.  7(a)  the  thermal  resistance  was 
quite  high,  which  indicates  that  the  heat  pipe  was  partially  dried  out  for  G  =  0.5,  even  at 
the  lowest  power  input  levels.  However,  the  thermal  resistance  decreased  significantly  as  the 
radial  acceleration  increased,  showing  that  the  capillary  pumping  abiHty  of  the  hefical  grooves 
increased.  For  G  =  1.0  and  1.5,  the  thermal  resistance  decreased  and  then  increased  with 
transported  heat  when  dryout  commenced.  The  G  —  1.5  fill  tests  showed  dryout  occurring 
only  for  la^l  =  0.01  and  2.0-g.  Dryout  was  not  reached  for  G  =  1.5  with  |ar|  =  4.0,  6.0,  8.0 
and  10.0-g  due  to  reaching  the  maximum  allowable  heater  temperature.  The  capillary  limit 
was  considered  to  be  reached  when  the  thermal  resistance  began  to  increase. 

Thomas  et  al.  (1998)  presented  a  mathematical  model  which  predicted  the  capillary  limit 
of  a  helicaUy-grooved  heat  pipe  subjected  to  a  transverse  body  force.  This  model  accounted 
for  the  geometry  of  the  heat  pipe  and  the  grooves  (including  helix  pitch),  body  force  field 
strength,  and  temperature-dependent  working  fluid  properties.  This  model  was  updated 
to  include  the  improved  measurements  of  the  wick  geometry  and  working  fluid  properties. 
Castle  (1999)  provides  a  detailed  comparison  of  the  present  model  and  that  given  by  Thomas 
et  al.  (1998).  The  capillary  fimit  versus  radial  acceleration  is  given  in  Fig.  8  for  various 
working  temperatures  with  the  Thomas  et  al.  (1998)  model  and  the  present  model.  The 
capillary  limit  increased  steadily  with  radial  acceleration  and  working  temperature.  The 
present  model  shows  a  significantly  lower  prediction  for  the  capillary  limit  when  compared  to 
the  Thomas  model  due  to  the  improved  geometric  measurements  and  working  fluid  property 
equations. 

Figure  9  shows  a  comparison  of  the  experimental  data  and  present  analytical  model  for 
the  capillary  limit  of  a  revolving  hefically-grooved  heat  pipe.  In  the  legend  of  this  figure. 
No  Dryout  and  Partial  Dryout  refer  to  the  experimental  data.  No  attempt  was  made  to 
maintain  a  constant  adiabatic  temperature  during  the  experiments.  Therefore,  the  working 
fluid  temperature  in  the  model  was  set  to  the  adiabatic  temperature  found  experimentally. 
For  G  =  0.5,  the  heat  pipe  operated  successfully  only  for  |ar|  >  8.0-g.  In  Fig.  9(6)  (G  = 
1.0),  the  capillary  fimit  increased  significantly  with  radial  acceleration.  With  the  heat  pipe 


14 


overfilled  by  50%  {G  =  1.5),  the  capillary  limit  increased  dramatically,  showing  the  effect 
that  overfilling  has  on  performance.  The  agreement  of  the  analytical  model  was  very  good  for 
G  =  1.0  as  expected.  For  G  =  0.5,  the  model  overpredicted  the  experimental  data  because  it 
was  assumed  that  the  grooves  were  completely  filled.  For  G  —  1.5,  the  model  underpredicted 
the  data  due  to  the  assumption  that  no  hquid  communication  occurred  between  the  grooves. 

Temperatures  within  the  evaporator  section  are  shown  in  Figs.  10  and  11  for  la^l  =  0.01-g 
and  10.0-g,  respectively.  In  general,  the  temperatm-es  within  the  evaporator  increased  with 
transported  heat.  In  addition,  the  wall  temperatures  decreased  with  G  for  a  given  heat  trans¬ 
port  due  to  the  fact  that  more  grooves  were  active.  The  temperatures  along  the  length  of  the 
evaporator  section  can  be  tracked  by  examining  the  case  for  G  =  1.0.  Near  the  evaporator 
end  cap,  the  temperatures  departed  those  for  G  =  1.5  at  approximately  =  15  W  (Fig. 
10(a)).  At  a  =  92.1  mm  (Fig.  10(6)),  this  departure  was  delayed  until  approximately  Qt  = 
25  W,  and  at  a;  =  168  mm  (Fig.  10(d)),  the  data  for  G  =  1.0  and  1.5  were  nearly  coincident. 
This  behavior  shows  that  the  grooves  were  essentially  full  near  the  adiabatic  section,  and 
proceeded  to  dry  out  closer  to  the  evaporator  end  cap,  as  expected.  Dryout  for  the  G  =  1.5 
case  can  be  seen  in  Fig.  10(a)  where  the  temperatures  converged  to  nearly  the  same  value 
around  the  circumference.  It  should  be  noted  that  the  temperatures  around  the  circumfer¬ 
ence  were  relatively  imiform  for  |dr|  =  0.01-g.  Evaporator  temperatures  for  |dr|  =  10.0-g 
are  shown  in  Fig.  11.  In  comparison  to  |dr|  =  0.01-g,  the  evaporator  temperatures  were  in 
general  lower  due  to  the  improved  pumping  ability  of  the  helical  grooves  under  increased 
radial  acceleration.  In  addition,  the  temperatures  tended  to  overlap  over  a  greater  range  of 
heat  transport  values.  In  contrast  to  |dr|  =  0.01-g,  the  evaporator  temperature  variation 
was  greater  around  the  circumference  at  higher  Qt,  but  no  pattern  was  distinguishable  in 
the  data. 

Local  heat  transfer  coefficient  data  versus  heat  transport  is  shown  in  Figs.  12  and  13 
for  |or|  =  0.01-g  and  10.0-g.  Overall,  the  values  for  he  were  very  low  for  G  =  0.5  due  to 
the  fact  that  most  of  the  grooves  were  dried  out.  As  the  percent  fill  increased  from  G  = 
0.5  to  G  =  1.0,  the  heat  transfer  coefficient  increased  significantly.  For  |arl  =  0.01-g  (Fig. 
12),  he  increased  and  then  decreased  with  transported  heat.  This  trend  was  also  reported  by 
Vasihev  et  al.  (1981)  for  an  aluminum  axially-grooved  heat  pipe  with  acetone  as  the  working 
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fluid.  For  G  —  1.0  and  1.5,  the  heat  transfer  coefficient  near  the  evaporator  end  cap  (Fig. 
12(a))  decreased  until  all  of  the  values  around  the  circumference  converged.  Closer  to  the 
adiabatic  section,  the  heat  transfer  coefficient  values  around  the  circumference  had  not  yet 
converged,  showing  these  portions  to  still  be  active.  For  |or|  =  10.0-g  (Fig.  13),  the  values  of 
he  were  signiflcantly  more  uniform  around  the  circumference  and  along  the  axial  direction, 
even  during  a  dryout  event  (G  =  1.0,  Fig.  13(a)).  In  addition,  the  heat  transfer  coefficient 
seems  to  be  more  constant  with  respect  to  the  transported  heat  compared  to  |ar-|  =  0.01-g. 
During  the  experiments,  the  heat  pipe  working  temperature  was  not  constant,  which  resulted 
in  changes  in  the  specific  volume  of  the  Hquid  and  vapor  of  the  working  fluid.  Since  the  heat 
pipe  was  filled  at  room  temperature,  it  was  important  to  quantify  the  potential  effects  of  the 
change  in  volume  of  liquid  in  the  grooves  with  temperature.  Figure  14  shows  the  variation 
of  the  percentage  of  groove  volume  occupied  by  hquid  G  with  saturation  temperature  for 
the  three  fiU  amounts  over  the  range  of  working  temperatures  seen  in  the  experiments.  The 
maximum  percent  difference  was  2.7%,  which  was  not  deemed  to  be  significant. 

1.7  Conclusions 

The  effect  of  fluid  inventory  on  the  steady-state  performance  of  a  helically-grooved 
copper-ethanol  heat  pipe  has  been  examined  both  experimentally  and  analytically.  It  was 
found  that  the  capillary  Hmit  increased  and  the  thermal  resistance  decreased  significantly 
as  the  amount  of  working  fluid  within  the  heat  pipe  increased.  In  addition,  the  evaporative 
heat  transfer  coefficient  was  found  to  be  a  strong  function  of  the  fiU  amount.  The  updated 
analytical  model  was  in  very  good  agreement  with  the  experimental  capillary  limit  results 
for  G  =  1.0.  However,  the  analytical  model  overpredicted  the  capillary  limit  data  for  G  = 
0.5  and  underpredicted  the  data  for  G  =  1.5. 

1.8  Nomenclature 

a  adiabatic  length  near  the  evaporator  end  cap,  m 

Ur  radial  acceleration,  m/s^ 

Ae  surface  area  in  the  evaporator  section,  ‘KrlLe,  m^ 

^gr  cross-sectional  area  of  a  groove,  m^ 
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b  adiabatic  length  near  the  condenser  end  cap,  m 

Do  tube  outside  diameter,  m 

jDvs  diameter  of  the  heat  pipe  vapor  space,  m 
G  Ratio  of  liquid  volume  to  total  groove  volume,  Vi/Vgr 

h  groove  height,  m 

he  local  heat  transfer  coefficient  in  the  evaporator  section,  Qt/Ae{T.u,  —  Ta 

W/m^-K 

I  heater  current,  A 

La  adiabatic  length,  m 

Lc  condenser  length,  m 

Le  evaporator  length,  m 

Lgr  hehcal  groove  length,  m 

Lt  total  heat  pipe  length,  m 

rrid  mass  of  working  fluid  dispensed  by  the  Ailing  station,  kg 

mi  mass  of  liquid,  kg 

mt  total  mass  of  working  fluid  inventory,  kg 

mass  of  vapor,  kg 

the  coolant  mass  flow  rate,  kg/s 

Ngr  number  of  grooves 

p  helical  pitch,  m 

Qcap  capillary  Hmit,  W 

Qin  heat  input  at  the  evaporator,  W 

Qt  heat  transported,  mcCp{ToMt  —  7i„),  W 

rh  radius  of  the  helix,  m 

r,,  radius  of  the  heat  pipe  vapor  space,  m 

R  radius  of  curvature,  m 

Rth  thermal  resistance,  (Teec  -  Tcec)/Qt,  K/W 

s  coordinate  along  the  centerHne  of  the  heat  pipe,  m 

tw  tube  wall  thickness,  m 

Ta  adiabatic  temperature,  K 
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Tcec  condenser  end  cap  temperature,  K 

'I'eec  evaporator  end  cap  temperature,  K 

Tin  calorimeter  inlet  temperature,  K 

Tout  calorimeter  outlet  temperature,  K 

Tsat  saturation  temperature,  K 

Tyj  outer  wall  temperature,  K 

vi,  Vy  specific  volume  of  liquid  and  vapor,  m^/kg 

V  heater  voltage 

Vgc  volume  of  the  grooves,  m^ 

Ks  vapor  space  volume,  m^ 

w  width  along  the  bottom  of  the  groove,  m 

X  distance  from  the  evaporator  end  cap,  m 

A  uncertainty 

^1,  02  angles  from  the  sides  of  the  groove  to  vertical,  rad 
(j)  helix  angle,  rad 


Subscripts 

cap  capillary 

gr  groove 

I  hquid 

max  maximum 

V  vapor 

vs  vapor  space 
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Figure  1:  Photomicrograph  of  the  helical  groove  geometry. 
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1 .  Angular  Transducer  4,  Shaft 

2.  Milling  Machine  Spindle  5.  Heat  Pipe  Container 

3 .  Milling  Machine  Collet  6.  Milling  Machine  Table 


Figure  2:  Schematic  of  the  helical  pitch  measurement  technique:  (a)  Major  components;  (b) 
Cross-sectional  view  of  sprung  pin  engaging  a  hehcal  groove. 
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Tsat  (°C) 

Figure  3:  Specific  volume  of  ethanol  versus  temperature:  (a)  Saturated  liquid;  (6)  Saturated 
vapor. 
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Figure  4:  Schematic  of  the  heat  pipe  filling  station. 
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Figure  6:  Steady-state  temperature  distributions  for  |orl  =  0.01-g,  G 
(b)  Outboard;  (c)  Top;  {d)  Bottom. 
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Figure  7:  Thermal  resistance  versus  heat  transport:  (a)  G  =  0.5;  (6)  G  =  1.0;  (c)  G  =  1.5. 
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Figure  8:  Capillary  limit  versus  radial  acceleration  comparison  of  present  model  and  Thomas 
et  al.  (1998). 
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Figure  9:  Comparison  of  present  model  and  experimental  capillary  limit  data  versus  radial 
acceleration:  (a)  G—  0.5;  (6)  G  =  1.0;  (c)  G  =  1.5. 


28 


(W/m^K)  h,  (W/m^K)  he  (W/m^K)  he  (W/m^K) 


2500 
2000 
1500 
1000 
500 
0 

0  20  40  60  80  li 

2500 
2000 
1500 
1000 
500 
0 

0  20  40  60  80  100 

2500 
2000 
1500 
1000 
500 
0 

0  20  40  60  80  11 

2500 
2000 
1500 
1000 
500 
0 

0  20  40  60  80  100 

Qtm 

Figure  12:  Heat  transfer  coefficients  within  the  evaporator  section  versus  transported  heat 
for  |ar|  =  0.01-g:  (a)  x  =  54.0  mm;  (6)  x  =  92.1  mm;  (c)  x  =  130  mm;  {d)  x  =  168  mm. 
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Figure  14:  Ratio  of  liquid  volume  to  total  groove  volume  versus  saturation  temperature. 
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2  Fully-Developed  Laminar  Flow  in  Trapezoidal  Grooves  with  Shear  Stress  at 

the  Liquid- Vapor  Interface 

2.1  Abstract 

This  paper  discusses  the  behavior  of  liquid  flowing  in  a  groove  with  a  trapezoidal  cross 
section.  For  fully  developed  laminar  flow,  the  conservation  of  mass  and  momentum  equations 
reduce  to  the  classic  Poisson  equation  in  terms  of  the  liquid  velocity.  A  finite  difference 
solution  was  employed  to  determine  the  mean  velocity,  volumetric  flow  rate,  and  Poiseuille 
number  (Po  =  /Re)  as  functions  of  the  groove  aspect  ratio,  groove-half  angle,  meniscus 
contact  angle  and  imposed  shear  stress  at  the  flquid-vapor  interface.  Comparisons  with 
existing  solutions  for  fully  developed  flow  in  rectangular  ducts  and  rectangular  and  triangular 
grooves  are  provided.  The  volumetric  flow  rate  in  a  groove  in  which  the  fill  amount  varies 
is  discussed.  A  semi-analytical  solution  and  a  two-point  numerical  solution  for  the  mean 
velocity  in  a  groove  are  presented  and  used  to  determine  the  capillary  limit  for  a  revolving 
helically-grooved  heat  pipe.  The  effects  of  interfacial  shear  stress  and  groove  fill  ratio  on 
heat  pipe  performance  are  investigated. 

2.2  Introduction 

Internally  grooved  ducts  are  used  in  process  equipment  to  improve  heat  transfer  during 
condensation  and  evaporation.  In  some  cases,  such  as  a  refrigeration  cycle  evaporator  or 
condenser  with  internal  grooves,  the  vapor  flow  is  cocmrrent  with  respect  to  the  liquid  flow. 
In  axially  grooved  heat  pipes,  the  vapor  flow  is  coimtercurrent  to  the  hquid  flow.  The 
interfacial  shear  stress  due  to  the  cocurrent  or  countercurrent  vapor  flow  contributes  to  the 
liquid  pressure  drop,  which  can  significantly  affect  the  heat  transfer  capacity  of  the  grooved 
surface.  The  objective  of  the  present  research  was  to  numerically  model  the  flow  of  hquid  in 
trapezoidal  grooves  using  a  finite  difference  approach  in  order  to  provide  accurate  information 
on  the  effects  of  cocurrent  and  countercurrent  vapor  flow  on  the  pressure  drop  in  the  hquid. 
This  geometry  was  chosen  due  to  the  fact  that  rectangular  and  triangular  grooves  are  special 
cases  of  the  trapezoidal  groove,  thus  making  the  analysis  as  general  as  possible. 

DiCola  [1]  solved  the  conservation  of  momentum  equation  for  the  laminar  flow  of  liquid 
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in  rectangular  grooves  with  a  uniform  shear  stress  imposed  at  the  Hquid-vapor  interface 
using  separation  of  variables.  Unfortunately,  this  manuscript  is  no  longer  available  in  the 
open  literature.  Schneider  and  DeVos  [2]  provided  the  exact  solution  determined  by  DiCola 
[1],  along  with  an  expression  for  the  friction  factor  which  approximates  the  exact  solution 
to  within  1%  by  using  the  first  term  of  the  infinite  series  solution.  This  expression  was 
used  by  Schneider  and  DeVos  to  determine  the  nondimensional  heat  transport  capacity  of 
axially-grooved  heat  pipes.  Upon  examination  of  the  DiCola  equation,  it  is  obvious  that  the 
rectangular  groove  is  completely  full;  i.e.,  the  meniscus  contact  angle  is  (f>  =  90°. 

Ayyaswamy  et  al.  [3]  used  the  Galerkin  boimdary  method  to  solve  the  Poisson  equation 
to  determine  the  fluid  velocity  in  triangular  grooves.  In  this  study,  interfacial  shear  stress  was 
zero,  and  the  groove  half-angle  and  contact  angle  were  varied  from  0  =  5  to  60°  and  0  =  0.1° 
to  the  full  groove  condition  (9  +  ^  =  90°).  The  results  were  presented  graphically  and  in 
tabular  form,  which  included  the  cross-sectional  area,  mean  velocity,  hydraulic  diameter,  and 
Poiseuille  number.  It  was  found  that  the  Poiseuille  number  increased  monotonically  with 
meniscus  contact  angle. 

Ma  et  al.  [4]  determined  the  Poiseuille  number  for  the  flow  of  liquid  in  triangular  grooves 
with  liquid-vapor  frictional  interaction.  The  groove  half-angle  ranged  from  20  <  0  <  60°  and 
the  meniscus  contact  angle  was  varied  from  0  =  0  to  60°.  The  conservation  of  momentum 
equation  was  transformed  into  the  Laplace  equation  in  terms  of  hquid  velocity,  which  was 
solved  using  separation  of  variables  and  linear  superposition.  Difficulties  were  encountered 
with  respect  to  application  of  the  liquid-vapor  interface  boundary  condition,  since  the  liquid- 
vapor  interface  velocity  is  an  unknown  function  of  the  vapor  velocity.  The  methodology  used 
to  overcome  this  difficulty  required  an  experimentally  determined  coefficient.  A  dimension¬ 
less  liquid-vapor  interface  flow  number  was  introduced  to  account  for  the  interfacial  shear 
stress.  This  value  determined  the  relative  velocities  of  the  Hquid  and  vapor  at  the  liquid- 
vapor  interface,  which  impacted  the  magnitude  of  the  PoiseuiUe  number.  It  was  found  that 
the  friction  factor  increased  with  the  interface  flow  number  and  contact  angle.  Results  from 
the  experiment  executed  by  Ayyaswamy  et  al.  [3]  for  no  Hquid-vapor  shear  stress  showed  an 
excellent  comparison  with  the  analytical  solution  over  the  ranges  of  groove  half-angle  and 
contact  angle  mentioned  above. 
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Peterson  and  Ma  [5]  performed  a  follow-up  analysis  to  that  of  Ma  et  al.  in  which 
the  Poiseuille  number  for  flow  in  triangular  grooves  was  determined  using  the  Nachtsheim- 
Swigert  iteration  scheme  and  a  numerical  solution  of  a  two-point  boundary  value  problem. 
This  allowed  the  velocity  distributions  in  both  the  liquid  and  vapor  to  be  solved  as  a  coupled 
problem.  Channel  angles  of  9  =  10°,  20°,  30°  and  40°  and  meniscus  contact  angles  of  0  = 
0°,  20°,  40°  and  60°  were  analyzed.  It  was  determined  that  the  dimensionless  Hquid-vapor 
flow  number  (previously  defined  by  Ma  et  al.)  was  a  function  of  the  vapor  velocity  and  the 
fluid  properties  of  the  flquid  and  vapor.  An  experiment  was  devised  to  validate  the  results  of 
the  numerical  model.  It  was  shown  that  increasing  the  countercurrent  vapor  flow  increases 
the  friction  factor  of  the  liquid.  The  agreement  between  the  experimental  data  and  the 
numerical  model  was  quite  good. 

Romero  and  Yost  [6]  analyzed  the  flow  of  flquid  in  a  triangular  groove  with  no  shear  stress 
at  the  flquid- vapor  interface.  In  particular,  the  flow  from  a  sessile  drop  into  the  groove  was 
of  interest.  A  nonlinear  partial  differential  diffusion  equation  was  presented  which  described 
the  time-dependent  height  of  flquid  in  the  groove  in  terms  of  the  groove  geometry,  meniscus 
contact  angle,  and  fluid  properties.  A  simplified  similarity  solution  was  presented  for  the 
region  which  was  far  from  the  sessile  drop.  A  full  similarity  solution  was  also  shown  which 
accounted  for  the  region  near  the  fluid  droplet.  It  was  found  that  the  wetting  front  position 
was  proportional  to  where  the  diffusion  coefficient  D  was  related  to  the  groove 

geometry,  flquid  viscosity,  and  flquid- vapor  surface  tension. 

Lin  and  Faghri  [7]  modeled  the  flow  of  flquid  in  the  triangular  grooves  of  a  rotating 
miniature  heat  pipe.  A  correlation  for  the  friction  factor  was  provided  in  terms  of  the  shear 
stress  at  the  liquid-vapor  interface.  The  laminar  flow  in  the  triangular  groove  was  solved 
using  a  finite  element  technique  for  side  lengths  ranging  from  hV  1  +  tan^  9  =  0.2  to  0.65 
mm  and  flquid- vapor  shear  stress  riv  =  7.7  x  10~^  to  0.055  N/m^  for  a  groove  half-angle  of 
9  —  20°  and  meniscus  contact  angle  0  =  30°.  A  regression  analysis  was  used  to  represent 
the  data  to  within  ±2.8%. 

Khrustalev  and  Faghri  [8]  analyzed  the  fully  developed  laminar  flow  of  flquid  and  vapor 
in  miniature  heat  pipes  using  a  finite  element  solution.  In  particular,  the  case  in  which  the 
vapor  velocity  was  high  and  the  cross-sectional  areas  of  the  vapor  and  flquid  were  comparable 
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was  of  interest.  It  was  assumed  that,  with  respect  to  the  vapor  flow,  the  liquid  velocity  at 
the  liquid-vapor  interface  was  zero.  For  the  hquid  flow,  the  shear  stress  at  the  interface 
was  equal  and  opposite  to  that  of  the  vapor.  This  meant  that  the  velocity  gradient  in  the 
liquid  was  related  to  that  of  the  vapor  via  a  ratio  of  absolute  viscosities.  The  momentum 
equation  for  the  vapor  domain  was  first  solved  to  determine  the  shear  stress  distribution  at 
the  liquid-vapor  interface.  Then  the  momentum  equation  for  the  Uquid  domain  was  solved 
using  the  shear  stress  information  from  the  vapor  solution.  The  results  presented  were  for  a 
specific  heat  pipe  geometry  that  matched  a  previous  experimental  study.  It  was  found  that 
the  shear  stress  at  the  Uquid-vapor  interface  was  not  uniform,  being  greater  near  the  point 
of  contact  with  the  solid  groove  waU.  This  effect  was  more  significant  for  smaller  values  of 
meniscus  contact  angle.  In  addition,  the  shear  stress  at  the  interface  became  more  uniform 
as  the  vapor  space  became  more  restricted. 

Kolodziej  et  al.  [9]  analyzed  the  gravity-driven  flow  of  liquid  in  a  triangular  groove  with 
no  shear  at  the  liquid-vapor  interface.  The  shape  of  the  liquid-vapor  interface  was  determined 
in  terms  of  the  Bond  number  and  meniscus  contact  angle.  Starting  with  the  Young-Laplace 
relation,  a  nonlinear  boundary-value  problem  for  the  liquid-vapor  interface  shape  was  solved. 
The  flow  field  was  then  solved  for  the  friction  factor  using  the  boundary  collocation  method. 
The  range  of  parameters  was  as  follows:  groove  half-angle  0  =  5  to  70°,  meniscus  contact 
angle  0  =  5  to  45°,  and  Bond  number  Bo  =  pgh^/o^  —  0.001,  0.01,  0.1,  1.0  and  5.0,  where  h 
is  the  vertical  distance  from  the  liquid-vapor  interface  to  the  bottom  of  the  groove.  It  was 
found  that  the  Bond  number  had  a  significant  effect  on  the  friction  factor  of  the  flow. 

The  objective  of  the  present  study  was  to  determine  the  mean  velocity,  volumetric  flow 
rate,  and  Poiseuille  number  for  the  flow  of  hquid  in  trapezoidal  grooves.  The  effect  of  vapor 
flowing  over  the  Hquid-vapor  interface  was  accounted  for  by  relating  the  liquid  velocity 
gradient  to  the  friction  factor  of  the  vapor.  This  approach  assumed  that  the  hquid  velocity 
did  not  affect  the  vapor  velocity;  i.e.,  the  vapor  velocity  at  the  interface  was  zero.  In 
addition,  the  variation  of  the  shear  stress  along  the  hquid-vapor  interface  was  neglected 
[8].  The  conservation  of  momentum  equation  was  solved  using  Gauss-Seidel  iteration  with 
successive  over-relaxation.  The  analysis  was  vahdated  using  the  results  of  several  previous 
studies  including  the  flow  of  hquid  in  rectangular  and  triangular  grooves  with  hquid-vapor 
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interaction.  The  mean  velocity,  volumetric  flow  rate,  and  Poiseuille  number  are  presented 
in  terms  of  the  groove  aspect  ratio  (0  <  ^  <  1.5),  groove  half-angle  {0  <6  <  60°),  meniscus 
contact  angle  (0  <  0  <  90°  —  6),  and  dimensionless  shear  stress  at  the  liquid- vapor  interface 
(—0.45  <  <  5.0).  The  results  were  used  to  determine  the  effects  of  groove  fill  amount  on 

the  capillary  limit  of  the  revolving  helically-grooved  heat  pipe  studied  by  Castle  et  al.  [1]. 
The  predictions  of  the  improved  capillary  limit  model  were  compared  to  the  experimental 
data  obtained. 

2.3  Mathematical  Model 


A  constant  property  liquid  flows  steadily  in  a  trapezoidal  groove  as  shown  in  Fig.  31. 
A  meniscus,  which  is  assumed  to  be  circular,  comprises  the  liquid-vapor  interface.  For  fully 
developed  laminar  flow  with  no  body  forces,  the  dimensionless  conservation  of  mass  and 
momentum  equations  reduce  to  [11] 

d'^v* 

On  the  groove  walls,  the  no-slip  condition  is  in  effect. 

'  0<x*  <  (3,  z*  =  0 

v*  =  0:  <  (3  <  X*  <  P  +  tan.9^  z*  =  {x*  —  P)  cot  9  for  9>0  (11) 

X*  =  P,  0  <  ^*  <  1  for  6>  =  0 

At  the  line  of  symmetry,  the  velocity  gradient  is  zero  in  the  x*  direction 

^  0  :  x*  =  0,  0  <  z*  <  (1  +  d*)  -  i?*  (12) 


where 


d*-=R*\  I 


'  P  +  tan^' 

.  'r*  . 


The  dimensionless  radius  of  curvature  is  given  by 


P  sec  0  for  d  =  0 
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On  the  liquid-vapor  interface,  a  uniform  shear  stress  in  the  y  direction  is  imposed. 
dv* 


dn* 


rl:  0  <  x*  <  yS  +  tan^,  z*  =  {l  +  (f)  -  y/R*^  - 


(15) 


The  dimensional  liquid-vapor  interface  shear  stress  can  be  cast  in  terms  of  the  friction  factor 
of  the  vapor. 

Pv  {Vy) 


Tlv  =  < 


Pv  {Vv) 


fv  for  cocurrent  flow 


fv  for  countercurrent  flow 


(16) 


The  Poiseuille  number  of  the  Hquid  in  the  groove  is  given  by 

D*'2 

Po  ==  /Re  ^ 

^  2v* 

The  dimensionless  hydraulic  diameter  is 


(17) 


Dl 


l3-^{p  +  tan9){l  +  d*)  -  R*^  cos  ^ 


{P  +  sec  9)  ^  ior  9  +  (f>  <  7r/2 


2  {2P  -F  tan  9)  (P  +  sec  9)  ^  ior  9  +  (f)  =  7r/2 

X 

The  mean  velocity  is  defined  as 

2  rp+tan  0  nz* 


2  rp+tan  V  rz 

/ 

Jo  Jo 


V*  dz*  dx* 


(18) 


(19) 


where  the  dimensionless  cross-sectional  area  is  given  by 


P+  {P  +  tan  0)  (1  -t-  d*)  —  R*"^  cos 


R* 


ioT  9  +  <p  <  7r/2 


(20) 


2P  +  tan  9  for  9  +  (f)  =  7r/2 
2.4  Numerical  Model 


The  elliptic  Poisson  equation  given  in  eqn.  (44)  with  mixed  boundary  conditions  [eqns. 
(45),  (46)  and  (50)]  was  solved  using  Gauss-Seidel  iteration  with  successive  over-relaxation 
and  a  central  differencing  scheme  [12].  The  convergence  criteria  for  the  iterative  solution 
was  set  to  e  =  10“^  for  each  case.  A  grid  independence  check  was  made  in  which  the  number 
of  grids  in  each  direction  was  doubled.  When  the  value  for  the  Poiseuille  number  did  not 
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change  by  more  than  3%,  grid  independence  was  considered  to  be  reached.  The  convergence 
criteria  was  then  reduced  by  an  order  of  magnitude  while  maintaining  the  highest  number  of 
grids.  If  the  Poiseuille  number  did  not  change  by  more  than  2%,  the  solution  was  considered 
to  be  independent  of  both  grid  size  and  e.  Otherwise,  a  grid  independence  check  was  made  at 
the  smaller  value  of  e  until  a  converged  solution  was  obtained.  In  fact,  the  grid  independence 
for  423  of  the  446  data  points  reported  was  less  than  1%  [13]. 

The  numerical  model  was  tested  against  several  existing  solutions,  such  as  rectangular 
ducts,  triangular  grooves  without  interfacial  shear  stress,  and  rectangular  and  triangular 
grooves  with  shear  stress.  See  Lykins  [13]  for  details. 

Shah  [14]  determined  the  friction  factors  for  the  laminar  flow  within  ducts  of  various  cross 
sections  using  a  least-squares-matching  technique.  A  comparison  was  made  of  the  Poiseuille 
number  between  the  present  solution  and  those  given  by  Shah  [14]  and  Shah  and  London  [15] 
for  laminar  flow  in  a  family  of  rectangular  ducts  (0  =  0°,  0.01  <  P  <  1.0).  The  agreement 
was  excellent,  with  a  maximum  difference  of  0.9%. 

Ayyaswamy  et  al.  [3]  presented  the  friction  factors  obtained  for  laminar  flow  in  triangular 
grooves  using  the  Galerkin  method  of  solution.  Romero  and  Yost  [6]  derived  an  equation 
for  the  dimensionless  volumetric  flow  rate  of  liquid  in  a  triangular  groove  using  asymptotic 
methods  and  a  regression  analysis.  Kolodziej  et  al.  [9]  used  the  boundary  collocation  method 
to  solve  the  same  problem,  except  that  the  liquid-vapor  interface  was  not  assumed  to  be 
circular.  The  present  solution  was  compared  to  that  obtained  by  Ayyaswamy  et  al.  for 
6=5°  and  60°  for  the  full  range  of  meniscus  contact  angle  (0.1°  <  (p  <  90°  —  6).  The 
agreement  was  excellent  for  6  =  60°,  but  was  less  so  for  9  =  5°.  This  was  due  to  the  extreme 
narrowness  of  the  flow  field  for  this  case.  The  maximum  percent  differences  for  9  =  5°  and 
60°  were  3.7%  and  0.9%,  respectively.  In  comparison  to  the  results  by  Romero  and  Yost  [6], 
the  maximum  percent  differences  were  2.2%  for  6  =  5°  and  2.3%  for  6  =  60°.  The  agreement 
with  the  results  provided  by  Kolodziej  et  al.  for  the  lowest  value  of  Bond  number  presented 
(Bo  =  0.001)  was  less  satisfactory,  with  a  maximum  percent  difference  of  4.5%  for  9  =  5°, 
and  19.9%  for  6  =  60°.  This  may  be  due  to  the  approximate  nature  of  the  solution  by 
Kolodziej  et  al.,  which  was  in  terms  of  a  truncated  infinite  summation. 

DiCola  [1]  presented  the  solution  for  the  Poiseuille  number  for  the  laminar  flow  of  a 
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constant  property  fluid  within  a  rectangular  groove.  While  interfacial  shear  stress  at  the 
Uquid-vapor  interface  was  accounted  for,  the  groove  was  assumed  to  be  completely  fuU,  with 
a  meniscus  contact  angle  of  ^  =  90°.  The  comparison  between  the  equation  by  DiCola  and 
the  results  of  the  present  analysis  for  laminar  flow  in  a  family  of  rectangular  grooves  at  the 
fuU  groove  condition  (cf)  =  90°,  0.1  <  ^0  <  1.0,  Tj*  =  —0.1,  0.0,  and  1.0)  resulted  in  an 
excellent  agreement  with  a  maximum  percent  difference  of  2.3%. 

The  present  model  was  compared  to  the  correlation  presented  by  Lin  and  Faghri  [7], 
where  the  friction  factor  for  the  flow  of  liquid  in  triangular  grooves  with  liquid-vapor  shear 
was  presented.  Unfortunately,  not  enough  information  was  provided  by  Lin  and  Faghri  to 
precisely  determine  the  limits  of  applicability  for  their  equation.  Therefore,  the  correlation 
was  evaluated  over  a  fairly  wide  range  for  comparison  with  the  present  solution.  The  agree¬ 
ment  was  quite  good  between  0.075  <  — <  0.1,  where  the  maximum  percent  difference  in 
this  range  was  2.2%. 

2.5  Results  and  Discussion 
2.5.1  Parametric  Analysis 

A  numerical  study  has  been  completed  in  which  the  flow  field  in  a  trapezoidal  groove  has 
been  solved.  SpeciflcaUy,  values  of  the  mean  velocity,  Poiseuille  number,  and  volumetric  flow 
rate  are  reported  for  various  values  of  the  groove  aspect  ratio,  groove  half-angle,  meniscus 
contact  angle,  and  dimensionless  shear  stress  at  the  Uquid-vapor  interface.  Figure  16  shows  a 
contour  plot  of  the  dimensionless  mean  velocity  obtained  for  cocurrent  flow,  no  shear  stress, 
and  countercurrent  flow.  The  flow  behavior  changes  significantly  with  shear  stress,  where 
the  maximum  velocity  is  located  along  the  liquid-vapor  interface  for  T;*  =  5.0  and  0.0,  and 
within  the  interior  of  the  flow  field  for  =  —0.1.  For  countercurrent  flow  [Fig.  16(c)],  a 
region  of  reversed  flow  occurs  near  the  intersection  of  the  groove  wall  and  the  Uquid-vapor 
interface.  Depending  on  the  magnitude  of  the  countercurrent  shear  stress,  the  Uquid  velocity 
could  be  entirely  reversed. 

Figure  17  presents  the  mean  velocity  versus  shear  stress  at  the  Uquid-vapor  interface  for 
several  values  of  the  groove  half-angle.  The  range  of  the  meniscus  contact  angle  (0  <  ^  < 
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90°  —  6)  was  divided  equally  to  show  the  behavior  of  the  mean  velocity  with  0.  The  mean 
velocity  increases  linearly  with  shear  stress  since  the  flow  is  aided  by  .  In  addition,  v* 
increases  with  groove  half-angle  and  meniscus  contact  angle,  which  is  a  result  of  an  increase 
in  cross-sectional  area.  As  the  groove  half-angle  6  increases,  v*  becomes  more  sensitive  to 
the  meniscus  contact  angle  (j)  due  to  the  increased  length  of  the  perimeter  of  the  Hquid- vapor 
interface. 

Figure  33  shows  that  the  Poiseuille  number  decreases  dramatically  as  the  shear  stress  at 
the  liquid- vapor  interface  increases.  For  countercurrent  flow  (r^*  <  0),  the  Poiseuille  number 
is  a  strong  function  of  shear  stress  since  the  mean  velocity  approaches  zero.  In  addition, 
the  Poiseuille  number  decreases  with  increasing  meniscus  contact  angle  for  a  given  value  of 
shear  stress.  For  cocurrent  flow  (r^*  >  0),  the  Poiseuille  number  is  a  lesser  function  of  the 
shear  stress,  but  increases  signiflcantly  with  meniscus  contact  angle.  The  Poiseuille  number 
is  a  weak  function  of  the  groove  half-angle. 

The  volumetric  flow  rate  versus  shear  stress  for  various  meniscus  contact  angles  can  be 
seen  in  Fig.  19.  The  volumetric  flow  rate  and  mean  velocity  display  similar  trends.  The 
volumetric  flow  rate  is  a  linear  function  of  shear  stress,  and  increases  with  meniscus  contact 
angle  and  groove  half-angle.  The  flow  rate  is  slightly  more  sensitive  than  the  mean  velocity 
with  respect  to  the  meniscus  contact  angle.  The  mean  velocity  is  given  as  a  function  of  the 
interfacial  shear  stress  for  a  constant  meniscus  contact  angle  in  Fig.  20.  For  9  =  0  and 
30°,  the  mean  velocity  increases  and  then  decreases  with  groove  aspect  ratio.  This  point 
can  be  further  elucidated  in  Fig.  21,  which  presents  the  mean  velocity,  Poiseuille  number 
and  volumetric  flow  rate  for  0  =  30°  and  -q;  =  5.0.  As  mentioned  previously,  the  mean 
velocity  increases  and  then  decreases  with  P  ior  6  <  30°.  This  phenomenon  also  impacts 
the  Poiseuille  number  and  the  volumetric  flow  rate,  where  V*  attains  a  maximum  value  with 
respect  to  P  for  a  given  groove  half-angle  <  15°.  Figure  22  shows  v*,  Po,  and  V*  versus  9  for 
P  =  1.0  and  Tj*  =  5.0.  In  general,  these  functions  increase  signiflcantly  with  both  meniscus 
contact  angle  and  groove  half-angle,  except  for  the  Poiseuille  number  for  0  =  0°  and  (j)  <  20°. 
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2.5.2  Semi-Analytical  and  Two-Point  Numerical  Solutions  for  v* 

As  seen  in  Figs.  17  and  20,  the  mean  velocity  is  a  linear  function  of  the  imposed  shear  stress 
at  the  hquid-vapor  interface.  Since  a  direct  numerical  simulation  of  the  Hquid  flow  field  for 
a  number  of  values  of  the  shear  stress  is  computer  resource  intensive,  it  is  appropriate  to 
seek  a  semi-analytical  expression  for  v*.  Figure  23(a)  shows  the  definition  of  the  parameters 
involved,  where  the  mean  velocity  when  the  shear  stress  is  zero  (0,  Vq)  is  given  by  the 
numerical  solution.  The  value  for  the  liquid-vapor  shear  for  which  the  mean  velocity  is  zero 
(tj*  0,  0)  is  given  by  the  following  force  balance  analysis.  Figure  23(6)  shows  a  differential 
element  of  the  liquid  in  the  groove.  A  force  balance  between  the  pressure  drop  and  the  shear 
forces  at  the  hquid-vapor  interface  and  at  the  wall  results  in  the  following  relation. 

Py-^l  ~~  Py-i-dy-^l  “f  '^Iv-^lv  TujA-uj  —  0  (21) 

The  areas  over  which  the  shear  stresses  tiv  and  act  are  Ai^  =  Pivdy  and  A^  =  Pdy, 
respectively.  Using  these  areas  and  nondimensionahzing  gives 

a:  +  Tl^Pt,  -  =  0  (22) 

For  Poiseuille  flow  in  ducts  of  arbitrary  cross  section,  and  combined  Couette-Poiseuille  flow 

between  flat  plates,  the  shear  stress  at  the  wall  is  related  to  the  mean  velocity  of  the  fluid 
by  a  constant  [11].  In  the  present  analysis,  it  is  assumed  that  this  also  holds  for  the  flow  of 
liquid  in  a  trapezoidal  groove  with  an  imposed  shear  stress  at  the  liquid-vapor  interface. 

^  (23) 

It  should  be  noted  that  the  constant  G\  is  probably  a  function  of  the  groove  geometry 
and  meniscus  contact  angle.  However,  since  the  objective  of  this  analysis  is  to  determine 

the  liquid-vapor  shear  stress  when  the  mean  liquid  velocity  is  zero,  this  functionality  is 

unimportant. 

Substituting  this  relation  into  the  force  balance  equation  results  in  the  following  expres¬ 
sion  for  mean  velocity 

5?'=^(^r+T,;p,:)  (24) 


43 


where  the  perimeter  of  the  liquid-vapor  interface  is 


2f?*sin-^  for  ^  +  0  <  7r/2 

2(/3  +  tan^)  for  6  +  (f>  =  7r/2 


The  mean  velocity  is  zero  when  the  shear  stress  at  the  liquid-vapor  interface  is 


(25) 


(26) 


Figures  24(a)  and  24(6)  show  the  results  of  eqn.  (26).  The  numerical  results  shown  in  Figs. 
17  and  20  were  extrapolated  to  determine  the  values  for  shear  stress  at  the  hquid-vapor 
interface  when  v*  =  0.  The  prediction  given  by  eqn.  (26)  is  quite  good  given  the  simplicity 
of  the  closed-form  solution.  The  equation  for  the  normalized  mean  velocity  as  a  function  of 
the  shear  stress  is  given  by 

v'  =  v*/vq  =  1  —  r'  (27) 

where  r'  —  The  semi-analytical  solution  for  the  normalized  mean  velocity  is  shown 

in  Fig.  24(c)  with  the  corresponding  numerical  data  presented  in  Figs.  17  and  20.  Equation 
(27)  predicts  93%  of  the  data  to  within  ±30%  over  the  range  of  the  meniscus  contact  angle, 
groove  half-angle,  groove  aspect  ratio  and  liquid-vapor  shear  stress  examined  in  Figs.  17  and 
20. 

The  two-point  numerical  solution  of  u*  as  a  function  of  is  also  shown  in  Fig.  23(a). 
The  finite  difference  numerical  model  is  used  to  compute  the  mean  velocity  for  two  values 
of  hquid-vapor  shear  stress  [(0,  Ug)  and  u^].  The  equation  for  the  normalized  mean 
velocity  as  a  function  of  shear  stress  using  this  solution  is  given  by  eqn.  (27),  and  the  shear 
stress  when  the  mean  velocity  is  zero  is 


'^iv,o  — 


*lv^a 


{1  -  Vjvl) 


(28) 


2.5.3  Effect  of  Groove  Fill  Ratio 


Figure  25(a)  shows  the  case  when  liquid  evaporates  from  a  trapezoidal  groove.  Initially,  the 
groove  is  full  with  0  ±  0  =  90°.  The  contact  angle  decreases  until  the  minimum  meniscus 
contact  angle  for  the  particular  sofid-liquid  combination  is  reached.  Past  this  point,  the 
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meniscus  detaches  from  the  top  of  the  groove  and  recedes  until  the  lowest  part  of  the  meniscus 
reaches  the  bottom  of  the  groove  [16].  When  the  thickness  of  the  liquid  film  at  the  bottom  of 
the  groove  is  on  the  order  of  several  hundred  Angstroms,  forces  due  to  London-van  der  Waals 
interactions  with  the  surroimding  hquid  and  sohd  molecules  induce  instabihties  in  the  fluid 
[17].  These  instabilities  cause  the  liquid  in  the  groove  to  bifurcate  into  two  separate  flows  in 
the  corners  of  the  groove,  which  are  each  equivalent  to  the  flow  in  a  triangular  groove.  The 
hquid  in  the  two  corners  of  the  groove  wiU  continue  to  recede  until  it  is  depleted. 

The  dimensions  of  the  grooves  analyzed  by  Castle  et  al.  [1]  were  used  to  determine  the 
volumetric  flow  rate  of  ethanol  in  a  trapezoidal  copper  groove  as  a  function  of  the  amount  of 
hquid  in  the  groove.  The  geometric  values  of  the  parametric  analysis  depicted  in  Fig.  25(o) 
are  given  by  Lykins  [13].  Faghri  [7]  gives  =  7°  for  a  receding  meniscus  of  ethanol  on  copper. 
In  terms  of  the  present  analysis,  as  the  hquid  recedes  into  the  groove,  the  groove  aspect  ratio 
P  increases.  In  addition,  after  bifurcation  occurs  P  =  0,  and  the  groove  half-angle  9  changes 
to  that  of  the  corner  of  the  groove.  A  relation  for  the  point  at  which  bifurcation  occurs 

is  provided  where  the  liquid  is  assumed  to  bifurcate  when  the  lowest  part  of  the  meniscus 

'  « 

actually  reaches  the  bottom  of  the  groove.  The  radius  of  curvature  at  the  bifurcation  point 
for  a  trapezoidal  groove  is 

=  (29) 


The  meniscus  contact  angle  at  bifurcation  as  a  function  of  the  groove  geometry  is 

-1 


(j)b  =  tan 


-1 


1 


sin0 


COS0  + 


2  sin  9{P  +  tan  9) 

1  —  {P  +  tan  9)^  J 


T-l 


—  COS0 


(30) 


The  radius  of  curvature  of  the  liquid-vapor  interface  is  shown  in  Fig.  26(o)  as  a  function 
of  the  groove  All  ratio.  When  the  groove  is  nearly  full,  the  radius  of  curvature  approaches 
infinity.  As  the  amount  of  Hquid  in  the  groove  decreases,  the  radius  of  curvature  is  relatively 
constant,  and  then  becomes  very  small  after  the  Hquid  bifurcates  into  the  corners  of  the 
groove.  Figures  26(6)  and  26(c)  show  two  flow  parameters,  Fi  and  F2,  which  allow  the 
presentation  of  the  mean  velocity  and  volumetric  flow  rate  before  and  after  bifurcation  on 
the  same  graph.  Both  flow  parameters  increase  monotonicaUy  with  groove  fiU  ratio,  as 
expected.  In  Fig.  26(c),  for  a  groove  fifl  ratio  of  Ai/Ag  =  0.158,  the  volumetric  flow  rate 
was  1%  of  that  for  the  fuU  groove  due  to  the  decrease  in  flow  area.  This  figure  shows  that 
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the  groove  was  essentially  shut  off  for  Ai/Ag  <  0.158,  which  was  prior  to  bifurcation.  Using 
the  functional  relationship  for  mean  velocity  given  by  eqn.  (24),  the  volumetric  flow  rate  is 

=  (V  +  <K)  (31) 

When  Tj*  =0,  V  oc  A^,  which  confirms  that  the  volumetric  flow  rate  should  decrease  rapidly 
with  decreasing  flow  area  for  .4;  <  1. 

2.5.4  Capillary  Limit  Analysis  for  a  Revolving  Helically-Grooved  Heat  Pipe 

Using  the  results  of  the  numerical  analysis,  the  capillary  hmit  prediction  for  a  revolving 
hehcaUy-grooved  heat  pipe  proposed  by  Thomas  et  al.  [3]  was  improved  by  accounting  for 
the  effects  of  working  fluid  fiU  amount  and  the  shear  stress  at  the  liquid-vapor  interface.  The 
improved  model  was  compared  to  the  experimental  data  collected  by  Castle  et  al.  [1],  who 
determined  the  capillary  hmit  of  a  revolving  helically-grooved  copper-ethanol  heat  pipe  for 
radial  accelerations  of  |d^|  =  0.01,  2.0,  4.0,  6.0,  8.0  and  10.0-g  and  groove  fill  ratios  of  Vi/Vg 
=  0.5,  1.0  and  1.5.  The  dimensions  of  the  heat  pipe  examined  by  Castle  et  al.  [1]  are  given 
in  Table  5.  A  pressiue  balance  within  the  heat  pipe  results  in  the  following  expression  for 
the  capillary  limit  ([7],  [20]) 


^Pcap  ,max  ^  Ap„  -I-  Api  +  Apbf 


(32) 


The  maximum  capillary  pressure  for  an  axial  groove  is 


Ap, 


cap, max 


<7 

Re 


(33) 


It  is  assumed  that  the  capillary  Hmit  occurs  when  the  liquid  bifurcates  into  the  corners  of 
the  grooves.  This  statement  is  based  on  the  results  of  the  variation  of  the  volumetric  flow 
rate  with  groove  fill  ratio  in  Fig.  26.  The  capillary  radius  when  the  hquid  bifurcates  (Rc  = 
0.02252  cm)  was  found  using  eqns.  (29)  and  (30). 

For  a  circular  cross  section  heat  pipe  with  uniform  heat  input  and  output  along  the 
lengths  of  the  evaporator  and  condenser,  respectively,  the  pressure  drop  in  the  vapor  is 


= 


SflyLe^Qt 

'K  Pyhi^Ry 


(34) 
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The  pressure  drop  in  the  Hquid  was  found  by  using  the  normahzed  mean  velocity  relation 
[eqn.  (27)]  rewritten  in  dimensional  form. 


dpi  = 


"»■  “r  » 


dy 


(35) 


In  a  heat  pipe,  the  liquid  flows  opposite  to  the  vapor  in  all  regions.  Therefore,  the  shear  stress 
at  the  liquid- vapor  interface  for  countercurrent  flow  was  used  [eqn.  (16)].  The  Poiseuille 
number  of  the  vapor  flow  was  modeled  as  laminar  flow  within  a  smooth  tube  with  a  circular 
cross  section  (Po„  =  16).  Substituting  these  relations  into  eqn.  (35)  gives 


dpi  =  - 


p-m 


dy 


(36) 


hfvQ  hiRvTi^  Q  ^ 

It  is  assumed  that  the  cross-sectional  area  of  the  hquid  is  constant  along  the  length  of  the 
groove.  For  a  constant  heat  flux  in  both  the  evaporator  and  condenser  sections,  eqn.  (36) 
can  be  integrated  to  determine  the  total  pressure  drop  in  the  hehcal  groove. 


A  „  T  {  Pl'^l.rnax 

^Pl  =  -f'eff  -T^=T-  - 


\  Kfvl  hiR^TlQ 
The  maximum  Hquid  velocity  in  a  groove  is 


\  L,  ) 


+1 


(37) 


Similarly,  the  maximum  vapor  velocity  is 

^v,inax  ■ 


PlAlhfg 


(38) 


Qt 


(39) 


irRyPuhfg 

where  the  total  heat  transported  by  the  heat  pipe  accounts  for  the  contributions  by  aU  of 
the  individual  grooves.  Using  the  above  relations,  the  Hquid  pressure  drop  in  a  groove  as  a 
function  of  the  transported  heat  is 


^  Teff  I  plQg  4^P>vQt 

'  hfg  yhfl^piAi  irhiRlpyT^lo 


+  1 


(40) 


The  body  forces  imposed  on  the  fluid  within  a  particular  groove  may  either  aid  or  hinder 
the  return  of  the  fluid  to  the  evaporator,  depending  on  the  groove  pitch  Lp  and  the  circum¬ 
ferential  location  of  the  starting  point  of  the  hehcal  groove  ([3],  [2]).  However,  even  if  the 
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body  force  hinders  the  return  of  the  fluid,  each  groove  contributes  to  the  heat  transported 
Qf  Therefore,  the  capillary  limit  equation  [eqn.  (32)]  was  first  solved  for  the  heat  trans¬ 
ported  by  each  individual  groove  Qg,  and  the  results  were  summed  to  determine  the  total 
heat  transport  Qt-  Since  the  pressure  drop  in  the  vapor  space  and  the  pressure  drop  in  each 
groove  were  functions  of  the  total  heat  transport,  eqn.  (32)  was  solved  iteratively.  The  body 
forces  due  to  acceleration  and  gravity  were  integrated  over  the  length  of  the  groove  to  find 
the  average  pressure  drop  [3]. 


Combining  the  above  relations,  the  general  expression  for  the  maximum  capillary  limit 
for  a  single  helical  groove  which  accounts  for  shear  stress  at  the  liquid-vapor  interface  and 
the  effect  of  groove  fiU  ratio  is  given  by 


A  closed-form  solution  for  the  capillary  limit  of  a  heat  pipe  with  straight  axial  grooves 
and  no  body  forces  can  be  derived  from  eqn.  (42). 


(j/ifg 

■^cTeff 


_ fM _ ^fJ’v 

NghfvoPiAi  -KPvRt 


Ro  \\ 


(43) 


Figure  27(a)  shows  the  results  of  the  closed-form  solution  [eqn.  (43)].  The  groove  geometry 
given  by  Castle  et  al.  [1]  was  used,  except  that  straight  axial  grooves  were  assumed  {Lp  oo) 
instead  of  helical  grooves.  Over  the  range  of  groove  fill  ratio  examined,  the  capillary  limit 
increased  with  Ai/Ag  by  more  than  three  orders  of  magnitude.  For  this  case,  the  semi- 
analytical  solution  and  the  two-point  numerical  solution  were  nearly  identical  due  to  the 
low  vapor  velocities,  and  hence  the  low  liquid-vapor  shear  stress.  This  point  is  further 
demonstrated  in  Fig.  27(a)  by  the  graph  indicated  by  “No  Shear,”  where  the  term  in  eqn. 
(43)  that  accounts  for  the  effect  of  Hquid-vapor  shear  on  the  liquid  was  dropped  by  allowing 
0  “OO-  The  capillary  hmit  decreases  when  shear  stress  is  accounted  for,  as  expected. 
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Figure  27(6)  shows  the  closed-form  solution  when  water  is  the  working  fluid.  In  this  case, 
the  difference  between  the  results  of  the  two-point  numerical  solution  and  the  no-shear 
solution  is  much  more  pronounced  due  to  the  significantly  higher  vapor  velocities  involved. 
The  agreement  between  the  semi-anal3d,ical  solution  and  the  two-point  numerical  solution  is 
quite  good.  The  semi-analytical  solution  offers  very  close  results  with  a  significantly  reduced 
amount  of  computer  resources  required. 

The  capillary  hmit  prediction  for  the  helically-grooved  heat  pipe  given  by  eqn.  (42)  is 
shown  in  Fig.  28  using  the  semi-analytical  solution  [eqn.  (26)]  for  jorj  =  0.0  and  10.0-g.  The 
capillary  fimit  heat  transport  increases  both  with  groove  fill  ratio  and  working  temperature. 
In  addition,  Qcap  increases  significantly  with  the  radial  acceleration  due  to  the  improved 
hquid  pumping  ability  of  the  helical  grooves  [3]. 

The  experimental  data  collected  by  Castle  et  al.  [1]  for  the  capillary  hmit  of  a  revolving 
hehcally-grooved  heat  pipe  versus  radial  acceleration  rates  are  shown  in  Fig.  29,  along  with 
the  predictions  of  the  present  semi-analytical  model  and  that  given  by  Castle  et  al.  [Ij. 
During  the  experiments,  the  working  temperature  was  not  held  constant,  so  the  present 
model  was  evaluated  at  the  saturation  temperature  reported  for  a  given  value  of  radial 
acceleration.  For  a  groove  fill  ratio  of  Ai/Ag  =  0.5  [Fig.  29(o)],  the  present  model  more 
closely  matches  the  experimental  data  than  the  model  by  Castle  et  al-  [1],  which  did  not 
account  for  the  groove  fill  ratio  or  liquid- vapor  shear  stress.  For  Ai/Ag  =  1.0  [Fig.  29(6)], 
the  present  model  significantly  overpredicts  the  experimental  data  and  the  previous  model, 
but  matches  the  trend  in  the  data  quite  weU,  given  the  simpficity  of  the  model. 

2.6  Conclusions 

A  numerical  study  has  been  concluded  where  the  mean  velocity,  Poiseuille  number,  and 
volumetric  flow  rate  of  Hquid  in  a  trapezoidal  groove  have  been  determined  as  functions 
of  groove  geometry,  meniscus  contact  angle  and  shear  stress  at  the  Hquid-vapor  interface. 
The  mean  velocity  and  volumetric  flow  rate  have  been  shown  to  be  Hnear  functions  of  shear 
stress,  and  the  Poiseuille  number  is  a  strong  function  of  the  shear  stress  for  countercurrent 
flow.  A  semi- analytical  solution  and  a  two-point  numerical  solution  for  the  mean  velocity 
were  presented  and  used  to  predict  the  capillary  Hmit  of  a  revolving  heHcaUy-grooved  heat 
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pipe  for  various  groove  fill  ratios.  Interfacial  shear  stress  due  to  countercurrent  flow  in  a 
heat  pipe  decreases  the  maximum  heat  transport.  For  cases  in  which  the  vapor  velocities  are 
high,  this  effect  is  more  pronounced.  The  groove  fill  ratio  was  shown  to  have  a  significant 
impact  on  heat  pipe  performance.  Underfilling  the  heat  pipe  examined  by  10%  resulted  in  a 
decrease  in  the  predicted  capillary  limit  by  approximately  17  to  20%  for  water  and  ethanol, 
respectively. 

2.7  Nomenclature 

Ag  cross-sectional  area  of  the  groove,  m^ 

Ai  cross-sectional  area  of  the  liquid,  m^ 

Ai-a  area  of  the  liquid- vapor  interface,  m^ 

Aui  area  of  the  groove  wall,  m^ 

A\  Ailh^ 

dr  radial  acceleration  vector,  m/s^ 

A  acceleration  vector  at  any  point  in  the  helical  groove,  m/s^ 

b  distance  from  the  liquid-vapor  interface  to  the  bottom  of  the  groove,  m 

Bo  Bond  number,  pgb^/cr 

d*  parameter  defined  in  eqn.  (48) 

D  diflfusion  coefficient 

Dh  hydraulic  diameter,  4Ai/P,  m 

Dl  Du/h 

Fi  mean  velocity  parameter,  hfv*,  m^ 

F2  volumetric  flow  rate  parameter,  hfV*.,  m^ 

/  friction  coefficient,  2TZ/inP 

g  acceleration  due  to  gravity,  m/s^ 

h  groove  height,  m 

/ifg  heat  of  vaporization,  J/kg 

hi  height  of  the  liquid  in  the  groove  at  the  wall,  m 

La  adiabatic  length,  m 

Lc  condenser  length,  m 
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Le 

LeS 

Lp 

Lt 


n 

n* 

N, 

P 

P 

Plv 

p* 

Ptv 

Po 

Qca.p 

Q. 

Qt 

R 

Rc 

Rh 

Rv 

R* 

Rl 

Re 

t 

Tsa.t 

V 

V 

^2, max 
'^Vjinax 

V* 


evaporator  length,  m 

effective  heat  pipe  length,  Le/2  +  La  +  ic/2,  m 

pitch  length,  m 

total  heat  pipe  length,  m 

coordinate  normal  to  the  Hquid-vapor  interface 

n/h 

number  of  grooves 
pressure,  N/m^ 
wetted  perimeter,  m 

perimeter  of  the  hquid-vapor  interface,  m 
P/h 

Pivlh 

Poiseuille  number,  /Re 

capillary  limit  heat  transport,  W 

heat  transfer  due  to  a  single  groove,  W 

total  heat  transported,  Yl^iQg,i,  W 

radius  of  curvature  of  the  hquid-vapor  interface,  m 

capiUary  radius,  m 

radius  of  the  helix,  m 

radius  of  the  heat  pipe  vapor  space,  m 

R/h 

dimensionless  radius  of  curvature  at  bifurcation 
Reynolds  number,  pvDhIli 
time,  s 

saturation  temperature,  K 
y-direction  velocity,  m/s 
mean  y-direction  velocity,  m/s 
maximum  mean  hquid  velocity,  m/s 
maximum  mean  vapor  velocity,  m/s 
Hv/h^{-dp/dy) 
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V* 

K 

v’ 

Vg 

Vi 

V 

V* 

w 

Wi 

X,  y,  2 
X* 

z* 

P 

e 

9 

fJ' 

P 

a 


Tlv 


T* 

Hv 


T* 

nv,o 


Iv.i 


Tu 


(t> 

00 

06 


dimensionless  mean  y-direction  velocity 
dimensionless  mean  y-direction  velocity  when  r^*  =  0 
dimensionless  mean  y-direction  velocity  when  =  0 
dimensionless  mean  y-direction  velocity  when  Tj*  = 
normalized  mean  velocity,  v*Ivq 
total  groove  volume,  m^ 
liquid  inventory  volume,  m^ 
volumetric  flow  rate,  vAi,  m^/s 
IjV  l[h\-dpldy)\ 

width  of  the  bottom  of  the  groove,  m 

width  of  the  liquid  in  the  groove,  m 

coordinate  directions 

x/h 

z/h 

groove  aspect  ratio,  wl2h 
convergence  criteria 
groove  half-angle,  rad 
absolute  viscosity,  Pa-s 
density,  kg/m^ 
surface  tension,  N/m 

normalized  shear  stress  at  the  liquid- vapor  interface, 
shear  stress  at  the  hquid- vapor  interface,  N/m^ 

Tiv/h{-dp/dy) 

shear  stress  at  the  liquid- vapor  interface  when  i;*  =  0 
shear  stress  value  in  Fig.  23{a) 
average  shear  stress  at  the  wall,  N/m^ 

T^lhi-dpldy) 
meniscus  contact  angle,  rad 
minimum  meniscus  contact  angle,  rad 
meniscus  contact  angle  at  bifurcation,  rad 
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Table  5:  Specifications  of  the  heat  pipe  test  article  examined  by  Castle  et  al.  [1]. 


Groove  Height,  h 

0.03831  cm 

Groove  Base  Width,  w 

0.03445  cm 

Groove  Half-Angle,  6 

14.62° 

Evaporator  Length,  Le 

15.2  cm 

Adiabatic  Length,  La 

8.2  cm 

Condenser  Length,  Lc 

15.2  cm 

Vapor  Space  Radius,  Rv 

0.6795  cm 

Number  of  Grooves,  Ng 

50 

Helical  Pitch  Length,  Lp 

135.8  cm 

Helix  Radius,  Rh 

0.6992  cm 

Figure  15:  Flow  of  liquid  in  a  trapezoidal  groove:  (o)  Coordinate  system;  (6)  Solution 
domain. 
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Figure  16:  Dimensionless  velocity  fields  for  laminar  flow  in  trapezoidal  grooves  (/3  =  1.0, 
(j)  =  10°,  9  =  45°):  (a)  =  5.0  {cocurrent  flow);  (b)  Til  =  0.0;  (c)  r,*  =  -0.1  (countercurrent 

flow). 
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Figure  18:  Po  versus  r*^  for  laminax  flow  in  trapezoidal  grooves  (/?  =  1.0):  (a)  0  =  0°;  {b) 
9  =  30°;  (c)  9  =  60°. 
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Figure  19:  V*  versus  r^*  for  laminar  flow  in  trapezoidal  grooves  {P  =  1.0):  (a)  9  =  0°;  {b) 
6  =  30°;  (c)  9  =  60°. 


Figure  24:  Numerical  and  semi-analytical  solution  comparison:  (a)  — _o  versus  0  for  /5  =  1.0 
(from  Fig.  17);  (b)  — T;*_o  versus  for  0  =  30°  (from  Fig.  20);  (c)  Normalized  mean  velocity 
versus  normalized  shear  stress  at  the  hquid- vapor  interface  (from  Figs.  17  and  20). 
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Figure  25:  Effect  of  groove  fill  ratio  on  liquid  flowing  in  a  trapezoidal  groove:  (o)  Parametric 
analysis  (to  scale,  h  —  0.03831  cm,  w  =  0.03445  cm,  9  =  14.62°);  (6)  Definition  of  variables 
prior  to  bifurcation  of  the  liquid;  (c)  Definition  of  variables  after  bifurcation  of  the  liquid. 
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Figure  26:  Effect  of  groove  fill  ratio:  (a)  Ptadius  of  curvature  of  the  liquid-vapor  interface; 
(6)  Mean  velocity  parameter;  (c)  Volumetric  flow  rate  parameter. 


67 


Figure  27:  Maximum  heat  transport  predicted  by  the  closed-form  solution  versus  groove  fill 
ratio  (Straight  axial  grooves,  no  body  forces,  Tgat  —  40°C):  (a)  Ethanol;  (b)  Water. 
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Figure  29;  Maximum  heat  transport  versus  radial  acceleration:  (o)  Ai/Ag 
=  1.0. 
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3  Fully-Developed  Laminar  Flow  in  Sinusoidal  Grooves 


3.1  Abstract 

The  flow  of  a  constant  property  fluid  through  a  sinusoidal  groove  has  been  analyzed.  A 
numerical  solution  of  the  conservation  of  mass  and  momentum  equations  for  fully  developed 
flow  is  presented.  The  mean  velocity,  volumetric  flow  rate  and  Poiseuille  number  are  pre¬ 
sented  as  functions  of  the  groove  geometry,  meniscus  contact  angle  and  shear  stress  at  the 
liquid-vapor  interface.  In  addition,  a  semi-analytical  solution  for  the  dimensionless  mean 
velocity  in  terms  of  the  dimensionless  shear  stress  at  the  meniscus  is  shown  to  agree  with 
the  numerical  data  quite  well. 

3.2  Introduction 

Isotropic  materials  such  as  quartz  glass  or  borosilicate  glass  can  be  chemically  etched 
to  form  micro-grooves  for  the  enhancement  of  evaporative  heat  transfer  in  chip-level  devices 
(Kirshberg  et  al.,  1999).  A  typical  etched  profile  in  glass  is  shown  in  Fig.  30,  where  the  profile 
is  smooth  instead  of  having  sharp  corners  seen  in  the  etching  of  crystalUne  materials  such  as 
silicon  (Maluf,  2000).  While  many  studies  have  been  performed  on  sharp-cornered  geometries 
[triangular  grooves  (Ayyaswamy  et  al.,  1974;  Ma  et  al.,  1994;  Romero  and  Yost,  1996;  Lin 
and  Faghri,  1997;  Kolodziej  et  al.,  1999),  rectangular  grooves  (DiCola,  1968;  Schneider  and 
DeVos,  1980;  Khrustalev  and  Faghri,  1999),  and  trapezoidal  grooves  (Hopkins  et  al.,  1999; 
Thomas  et  al.,  2000)],  very  little  information  is  available  in  the  open  literature  on  the  flow 
of  liquid  in  rounded-corner  geometries.  Stroes  and  Catton  (1997)  compared  the  capillary 
performance  of  triangular  and  sinusoidal  grooves  by  means  of  an  experimental  study.  Two 
sets  of  grooves  were  machined  into  stainless  steel  test  plates  such  that  the  cross-sectional 
areas  of  the  grooves  were  equal.  Strip  heaters  were  placed  under  the  plates  to  provide  heat 
input.  The  test  plates  were  placed  at  incHnation  angles  of  4  and  6°  and  ethanol  was  added 
to  the  grooves  until  the  liquid  reached  the  lands  of  the  grooves.  The  average  wetted  length 
of  each  set  of  grooves  was  recorded  as  the  heat  input  was  varied  from  0  to  25  W.  The  study 
showed  that  the  triangular  grooves  had  a  greater  capillary  pumping  abihty  compared  to 
the  sinusoidal  grooves  with  the  same  cross-sectional  area,  inclination  angle,  and  heat  input. 
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Stroes  and  Catton  postulated  that  this  was  due  to  the  axial  rate  of  change  of  the  radius 
of  curvature  of  the  meniscus.  Sinusoidal  grooves,  however,  could  dissipate  a  given  heat 
input  with  a  smaller  wetted  area  than  triangular  grooves  due  to  the  larger  wetted  perimeter 
found  with  sinusoidal  grooves.  The  objective  of  the  present  study  was  to  examine  the  fully- 
developed  laminar  flow  of  liquid  in  sinusoidal  grooves.  The  effects  of  countercurrent  and 
cocurrent  vapor  flow  over  the  liquid-vapor  interface  were  investigated  by  relating  the  liquid 
velocity  gradient  to  the  friction  factor  of  the  vapor.  The  variation  of  the  shear  stress  on 
the  hquid-vapor  interface  (Khrustalev  and  Faghri,  1999)  was  neglected,  and  the  hquid-vapor 
interface  was  assumed  to  be  circular  (Bo  1).  The  mean  velocity,  volumetric  flow  rate  and 
Poiseuille  number  were  determined  as  functions  of  the  interfacial  shear  stress,  the  meniscus 
contact  angle,  the  groove  aspect  ratio  and  the  amount  that  the  groove  was  filled. 

3.3  Mathematical  Model 


A  constant  property  liquid  flows  steadily  in  a  sinusoidal  groove  as  shown  in  Fig.  31.  A 
meniscus,  which  is  assumed  to  be  circular,  comprises  the  liquid-vapor  interface.  For  fully 
developed  laminar  flow,  the  conservation  of  mass  and  momentum  equations  reduce  to  the 
classic  Poisson  equation  in  dimensionless  form  (White,  1991) 


d'^v*  d\* 

dx*^  ^  dz*^ 


(44) 


On  the  groove  wall,  the  no-slip  condition  is  in  effect. 


0  <  a;*  <  Wil2, 


At  the  line  of  symmetry,  the  velocity  gradient  is  zero  in  the  x*  direction 
dv* 

—  =  0  :  a;*  =  0,  0  <  z*  <  hf  + 1?* 

dx*  ,  _  _  i 

The  dimensionless  radius  of  curvature  is  given  by 


1 

^  ^  {1  +  d*  tan0)^  * 
{d*  —  tan  0)^ 


where 


(45) 

(46) 

(47) 

(48) 
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The  maximum  value  for  the  meniscus  contact  angle  (j>  for  a  wetting  fluid  can  be  determined 
for  a  given  geometry  by  allowing  the  radius  of  curvature  to  approach  R*  —>■  oo. 

At  the  liquid-vapor  interface,  a  uniform  shear  stress  is  imposed  in  the  y  direction. 
dv* 


(49) 


=  r;„  :  0<x*  <  w;/2,  z*  =  h*i  +  R\  1 - 


W7 


VR*^  —  X 


*2 


(50) 


The  dimensional  shear  stress  at  the  liquid- vapor  interface  can  be  cast  in  terms  of  the  friction 
factor  of  the  vapor. 

Pv  (v^f 


riv 


fv 


for  cocurrent  flow 


^\2 


p.  K) 


(51) 


/t,  for  countercurrent  flow 


The  Poiseuille  number  of  the  liquid  in  the  groove  is  given  by 


Po  =  /Re  = 


2v* 


(52) 


The  dimensionless  hydraulic  diameter  for  the  flow  of  liquid  in  a  sinusoidal  groove  with  a 
circular  meniscus  is  =  AA*/P*,  where  the  dimensionless  cross-sectional  area  of  the  liquid 
is  given  by 


a:  (2h!  -  1)  -  R* 


R*  cos 


-1 


1- 


W7 


w7 


+ 


2R\ 
for  0  < 


2R* 


A*  =  ^(2/i;-l)  + 


\7r  / 


sm 


for  0  =  0n 


The  dimensionless  wetted  perimeter  is  given  by  the  following  integral  equation. 

w7/2 


P*  =  2  [ 

Jo 

The  mean  velocity  is  defined  as 


1  + 


TT  \  . 

-  1  SI 


sm 


-TT 


+  1 


dx* 


2  fz* 

V*  =  —  I  V*  dz*  dx* 

a;  Jo  Jo 


(53) 


(54) 


(55) 


(56) 
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Table  6:  Poiseuille  number  versus  sinusoidal  duct  aspect  ratio:  Comparison  of  the  present 
solution  with  that  given  by  Shah  (1975). 


Poiseuille  Nu 
Shah  (1975) 

mber,  Po 
Present 

1/4 

14.553 

14.479 

1/3 

14.022 

13.931 

1/2 

13.023 

12.935 

l/\/3 

12.630 

12.543 

2/3 

12.234 

12.148 

1 

11.207 

11.115 

2 

10.123 

10.061 

4 

9.743 

9.6373 

3.4  Numerical  Model 

The  elhptic  Poisson  equation  given  by  eqn.  (44)  with  mixed  boundary  conditions  [eqns. 
(45),  (46)  and  (50)]  was  solved  using  Gauss-Seidel  iteration  with  a  central  differencing  scheme 
and  successive  over-relaxation  (Anderson  et  ah,  1984).  The  convergence  criteria  for  the 
iterative  solution  was  set  to  e  =  10“^°  for  each  case.  A  grid  independence  check  was  made  in 
which  the  number  of  grids  in  each  direction  was  doubled.  When  the  value  for  the  Poiseuille 
number  did  not  change  by  more  than  1%,  grid  independence  was  considered  to  be  reached. 
The  convergence  criteria  was  then  reduced  by  an  order  of  magnitude  while  maintaining 
the  highest  number  of  grids.  If  the  Poiseuille  number  did  not  change  by  more  than  1%, 
the  solution  was  considered  to  be  independent  of  both  grid  size  and  e.  Otherwise,  a  grid 
independence  check  was  made  at  the  smaller  value  of  e  until  a  converged  solution  was  reached. 
The  integral  equation  for  the  wetted  perimeter  [eqn.  (55)]  was  integrated  numerically  since 
no  closed-form  solution  exists. 

The  numerical  model  was  tested  against  an  existing  solution  in  the  archival  hterature. 
Shah  (1975)  determined  the  friction  factors  for  the  laminar  flow  within  ducts  of  various  cross 
sections  using  a  least-squares-matching  technique.  Table  6  shows  the  comparison  of  the 
Poiseuille  number  between  the  present  solution  and  that  given  by  Shah  (1975)  for  laminar 
flow  in  a  family  of  sinusoidal  ducts.  The  agreement  is  excellent,  with  a  maximum  difference 
of  1.1%. 
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3.5  Results  and  Discussion 


A  numerical  study  has  been  completed  in  which  the  flow  of  liquid  in  a  sinusoidal  groove 
has  been  solved.  Figure  32  presents  contour  plots  of  the  dimensionless  hquid  velocity.  The 
maximum  liquid  velocity  increases  with  cocurrent  shear,  and  decreases  with  countercurrent 
shear,  as  expected.  For  countercurrent  vapor  flow,  a  portion  of  the  hquid  flows  in  the  —y 
direction,  which  is  opposite  to  the  direction  of  the  pressure  gradient.  This  flow  reversal 
shows  the  potential  of  the  vapor  shear  to  drive  the  mean  velocity  of  the  Hquid  to  zero,  or  to 
completely  reverse  the  flow,  depending  on  the  magnitude  of  the  pressure  gradient. 

Figure  33  and  Table  7  show  the  mean  velocity,  volumetric  flow  rate  and  Poiseuille  number 
versus  shear  stress  at  the  liquid-vapor  interface  for  several  values  of  the  meniscus  contact 
angle.  The  mean  velocity  increases  with  both  '’'iv  and  (f),  but  is  a  relatively  weak  function 
of  (f).  For  a  given  value  of  meniscus  contact  angle,  the  mean  velocity  is  linear  with  r^*  due 
to  an  overall  force  balance  on  the  hquid.  The  volumetric  flow  rate  also  increases  with 
and  (j),  but  is  a  much  stronger  function  of  <{>  due  to  the  increase  in  the  cross-sectional  area 
of  the  hquid.  The  Poiseuille  number  decreases  as  r*^  increases.  For  cocurrent  vapor  flow 
(t;*  >  0),  Po  decreases  steadily  with  Tj*  .  For  countercurrent  flow,  Po  increases  dramatically 
with  — T;*  due  to  the  mean  velocity  approaching  zero.  In  general,  Po  increases  with  0  due  to 
the  increase  in  the  hydraulic  diameter  of  the  hquid. 

Figure  34  and  Table  8  present  the  mean  velocity,  volumetric  flow  rate  and  Poiseuille 
number  versus  the  groove  flh  ratio  for  several  values  of  the  groove  aspect  ratio.  The  mean 
velocity  increases  monotonically  with  area  ratio  and  groove  aspect  ratio.  The  volumetric 
flow  rate  also  increases  with  wf/w*  and  P,  but  drops  off  rapidly  for  Wi/w*  <  0.5  due  to  the 
decrease  in  flow  area.  The  Poiseuihe  number  is  relatively  constant  for  this  case.  For  ^  >  0.5, 
Po  decreases  and  then  increases  with  area  ratio.  For  j3  =  0.25,  Po  increases  over  the  range 
of  Wi/w*  examined. 

3.6  Semi- Analytical  Solution  for  W 

As  seen  in  Fig.  33(a),  the  mean  velocity  is  a  hnear  function  of  the  imposed  shear  stress 
at  the  hquid-vapor  interface.  Since  a  direct  numerical  simulation  of  the  liquid  flow  field  for 
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Table  7:  Mean  velocity,  PoiseuiUe  number  and  volumetric  flow  rate  versus  shear  stress  at  the 
liquid- vapor  interface  for  various  values  of  meniscus  contact  angle  {P  =  0.5,  tuf/2  =  0.25,  P* 
=  1.15245). 


4> 

V* 

Po 

V* 

0° 

.0.1 

0.18878  X  10-^ 

27.088 

0.17394  X  10-^ 

0° 

-0.075 

0.26896  X  10-^ 

19.013 

0.24782  X  10-^ 

0° 

-0.05 

0.35104  X  10-^ 

14.567 

0.32344  X  10-“ 

0° 

-0.025 

0.43465  X  lO'-'^ 

11.765 

0.40048  X  10-^ 

0° 

0.0 

0.51353  X  10-^ 

9.9581 

0.47316  X  10-^ 

0° 

0.0625 

0.71345  X  10“=^ 

7.1676 

0.65738  X  10-“ 

0° 

0.125 

0.91525  X  lO--* 

5.5873 

0.84330  X  10-^ 

0° 

0.25 

0.13188  X  10-^ 

3.8774 

0.12152  X  10-^ 

0° 

0.375 

0.17224  X  10-* 

2.9689 

0.15870  X  10-^ 

0° 

0.5 

0.21260  X  10-^ 

2.4053 

0.19589  X  10-^ 

0° 

0.75 

0.29332  X  10-' 

1.7434 

0.27026  X  10-^ 

0° 

1.0 

0.37404  X  10-^ 

1.3672 

0.34464  X  10-^ 

0° 

2.0 

0.69691  X  10-^ 

0.73377 

0.64214  X  10-^ 

25“ 

-0.1 

0.32506  X  10-^ 

27.245 

0.39414  X  10-^ 

25“ 

-0.075 

0.41572  X  10-^ 

21.303 

0.50408  X  10-^ 

25“ 

-0.05 

0.50437  X  10-^ 

17.559 

0.61158  X  10-^ 

25“ 

-0.025 

0.59348  X  10-^ 

14.923 

0.71962  X  10-“ 

25“ 

0.0 

0.68259  X  10-^ 

12.974 

0.82768  X  10-“ 

25“ 

0.0625 

0.90536  X  10-^ 

9.7820 

0.10978  X  10-^ 

25“ 

0.125 

0.11229  X  10-1 

7.8870 

0.13616  X  10-^ 

25“ 

0.25 

0.15650  X  10-i 

5.6588 

0.18977  X  10-^ 

25“ 

0.375 

0.20072  X  10-^ 

4.4123 

0.24338  X  10-^ 

25“ 

0.5 

0.24493  X  10-^ 

3.6158 

0.29700  X  10-^ 

25“ 

0.75 

0.33336  X  10-* 

2.6566 

0.40422  X  10-‘‘ 

25“ 

1.0 

0.42179  X  10-^ 

2.0997 

0.51144  X  10-'“' 

25“ 

2.0 

0.77551  X  10-‘ 

1.1420 

0.94034  X  10"^ 

50“ 

-0.1 

0.45097  X  lO-'^ 

27.149 

0.64296  X  10-^ 

50“ 

-0.075 

0.54643  X  10-^ 

22.406 

0.77904  X  10-^ 

50“ 

-0.05 

0.64189  X  10-^ 

19.074 

0.91514  X  lO-'’ 

50“ 

-0.025 

0.73345  X  10-^ 

16.693 

0.10457  X  10-^ 

50“ 

0.0 

0.82746  X  10-^ 

14.796 

0.11797  X  10-^ 

50“ 

0.0625 

0.10625  X  10-^ 

11.523 

0.15148  X  10-^ 

50“ 

0.125 

0.12975  X  10-^ 

9.4363 

0.18498  X  10-^ 

50“ 

0.25 

0.17586  X  10-^ 

6.9621 

0.25072  X  10“^ 

50“ 

0.375 

0.22255  X  10-* 

5.5014 

0.31728  X  10-^ 

50“ 

0.5 

0.26922  X  10-^ 

4.5477 

0.38382  X  lO-'^ 

50“ 

0.75 

0.36256  X  10-^ 

3.3769 

0.51690  X  10-^ 

50“ 

1.0 

0.45590  X  lO-i 

2.6855 

0.64998  X  10-^ 

50“ 

2.0 

0.82926  X  10-^ 

1.4764 

0.11823  X  10-^ 

72.34° 

-0.1 

0.54127  X  10-^ 

28.188 

0.86146  X  10-^ 

72.34“ 

-0.075 

0.64152  X  10-^ 

23.783 

0.10210  X  10-^ 

72.34“ 

-0.05 

0.74177  X  10-^ 

20.569 

0.11805  X  10-^ 

72.34“ 

-0.025 

0.84202  X  10-^ 

18.120 

0.13401  X  10-^ 

72.34“ 

0.0 

0.94227  X  10-^ 

16.192 

0.14996  X  10-^ 

72.34“ 

0.0625 

0.11929  X  10-^ 

12.790 

0.18985  X  10-^ 

72.34“ 

0.125 

0.14354  X  10- * 

10.629 

0.22844  X  10-^ 

72.34“ 

0.25 

0.19309  X  10-^ 

7.9016 

0.30730  X  10-^ 

72.34“ 

0.375 

0.24264  X  10-i 

6.2880 

0.38618  X  10“^ 

72.34“ 

0.5 

0.29219  X  10-^ 

5.2217 

0.46504  X  lO-'^ 

72.34“ 

0.75 

0.39129  X  10-^ 

3.8992 

0.62276  X  10-^ 

72.34“ 

1.0 

0.49039  X  10-^ 

3.1112 

0.78048  X  lO-'-* 

72.34“ 

2.0 

0.88244  X  10-^ 

1.7290 

0.14044  X  10-^ 
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Table  8:  Wetted  perimeter,  mean  velocity,  PoiseuiUe  number  and  volumetric  flow  rate  versus 
groove  fill  ratio  for  various  values  of  groove  aspect  ratio  (r^*  =  0.0,  4>  =  0°). 


p 

w*ilw* 

P* 

V* 

Po 

V* 

0.25 

0.1 

0.0731191 

OEIiskiM 

8.4006 

WifcaEManiBa 

0.25 

0.223950 

9.6365 

0.70676  X  10-® 

0.25 

0.450806 

11.141 

0.25 

0.4 

0.734032 

12.059 

0.25 

0.5 

1.04707 

12.514 

0.25 

0.6 

1.36011 

12.849 

0.25 

0.7 

1.64333 

12.871 

■igbfelikmiM 

0.25 

0.8 

1.87019 

■igViefeBltiBi 

12.940 

0.25 

0.9 

2.02102 

13.346 

0.25 

0.999 

2.09364 

15.203 

0.5 

0.1 

0.114187 

9.4469 

0.5 

0.2 

0.288528 

8.4052 

■ifcaVifeCThM 

0.5 

0.531454 

■tlliiWtglliBi 

8.7511 

■iftiTittelUtiM 

0.5 

0.827649 

9.3706 

0.5 

0.5 

1.15245 

liMlilddBllBI 

9.9581 

0.5 

0.6 

1.47724 

■iWAfekglliBl 

10.398 

0.5 

0.7 

1.77344 

iiwyAkBUiM 

10.780 

0.5 

0.8 

2.01636 

11.282 

■ifcrfrfiVJSHiM 

0.5 

0.9 

2.19070 

0.5 

0.999 

2.30389 

■iHWgfekRlIiM 

0.75 

0.1 

0.160025 

■irttbki 

■iWlfeEKlIiM 

0.75 

0.2 

0.367435 

■ 

■ 

0.75 

0.3 

0.634966 

■ilikfeifcCTliBl 

8.4299 

0.75 

0.4 

0.951588 

8.5312 

■laMMKWBa 

0.75 

0.5 

1.29509 

8.8579 

0.75 

0.6 

1.63860 

ItHfcfeiabM 

9.2444 

0.75 

0.7 

1.95522 

9.7209 

■ifMigfellltiBa 

0.75 

0.8 

2.22275 

10.489 

0.75 

0.9 

2.43016 

11.885 

■ilirieWillltiai 

0.75 

0.999 

2.58868 

ligWriifeglbBl 

13.997 

1.0 

0.1 

0.207697 

ligbiiVJaitiM 

0.73744  X  10-^“ 

1.0 

0.2 

0.453879 

iiiifekMisnim 

1.0 

0.3 

0.752420 

■iirfifeggwniBa 

8.7481 

■mfcetlirailiBa 

1.0 

0.4 

1.09559 

WtVkkkltTmUWlM 

8.4926 

1.0 

0.5 

1.46369 

TOB^tlrgiriM 

8.5590 

■iI;gWrf;RiriBl 

1.0 

0.6 

1.83180 

8.8427 

1.0 

0.7 

2.17497 

■iMiWdJSibBl 

9.3068 

1.0 

0.8 

2.47351 

fSSSM 

■ilfcfelrifUlliM 

1.0 

0.9 

2.71969 

1.0 

0.999 

2.92539 

13.281 

■llaijlafeVvRUM 

77 


a  number  of  values  of  the  shear  stress  is  computer  resource  intensive,  it  is  appropriate  to 
seek  a  semi-analytical  expression  for  v*.  Figure  35(o)  shows  the  definition  of  the  parameters 
involved,  where  the  mean  velocity  when  the  shear  stress  is  zero  (vq)  is  given  by  the  numerical 
solution.  The  value  for  the  liquid-vapor  shear  for  which  the  mean  velocity  is  zero  (tj*  q)  is 
given  by  the  following  force  balance  analysis.  Figure  35(6)  shows  a  differential  element  of 
the  liquid  in  the  groove.  A  force  balance  between  the  pressure  drop  and  the  shear  forces  at 
the  hquid-vapor  interface  and  at  the  wall  results  in  the  following  relation. 

'Py^l  T  -  0  (b^) 

The  areas  over  which  the  shear  stresses  tu  and  act  are  =  Pivdy  and  —  Pdy, 
respectively.  Using  these  areas  and  nondimensionalizing  gives 

A*  +  r^P^v  -^P*  =  0  (58) 


For  PoiseuiUe  flow  in  ducts  of  arbitrary  cross  section,  and  combined  Couette-PoiseuiUe  flow 
between  flat  plates,  the  shear  stress  at  the  wall  is  related  to  the  mean  velocity  of  the  fluid 
by  a  constant  (White,  1991).  Therefore,  in  the  present  analysis,  it  is  assumed  that  this 
also  holds  for  the  flow  of  liquid  in  a  sinusoidal  groove  with  an  imposed  shear  stress  at  the 
liquid-vapor  interface. 

r*  =  Civ*  (59) 


It  should  be  noted  that  the  constant  Ci  is  probably  a  function  of  the  groove  geometry 
and  meniscus  contact  angle.  However,  since  the  objective  of  this  analysis  is  to  determine 
the  liquid-vapor  shear  stress  when  the  mean  Hquid  velocity  is  zero,  this  functionality  is 
unimportant.  The  perimeter  of  the  liquid-vapor  interface  is 

P,;  =  2ii*  sin->  (^)  (60) 


Substituting  these  relations  into  the  force  balance  equation  results  in  the  following  expression 
for  mean  velocity. 


1 

CiP* 


A*i  +  2t;^R*  sin-^ 


(61) 


The  mean  velocity  is  zero  when 


(62) 
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Figure  36(a)  shows  the  results  of  eqn.  (62).  The  numerical  results  shown  in  Fig.  33(o) 
were  extrapolated  to  determine  the  values  for  shear  stress  at  the  liquid- vapor  interface  when 
^  =  0.  Both  curves  indicate  that  Tj*  ^  increases  with  (j)  due  to  the  increasing  depth  of  liquid 
in  the  groove.  The  prediction  given  by  eqn.  (62)  is  quite  good  given  the  simplicity  of  the 
closed-form  solution.  The  equation  for  the  mean  velocity  as  a  function  of  the  shear  stress  is 
given  by 

v'  =  v*/vq  =  1  —  t'  (63) 

where  t'  —  Tj* The  semi-analytical  solution  for  the  normalized  mean  velocity  is  shown 
in  Fig.  36(6)  with  the  corresponding  numerical  data  presented  in  Fig.  33.  Equation  (63) 
predicts  94%  of  the  data  to  within  ±20%  over  the  range  of  meniscus  contact  angle  examined. 

3.7  Conclusions 

Based  on  the  results  of  the  numerical  model  of  the  flow  of  flquid  in  a  sinusoidal  groove, 
the  following  conclusions  have  been  made:  For  a  given  meniscus  contact  angle,  the  mean 
velocity  was  linear  with  imposed  shear  stress  at  the  liquid-vapor  interface.  The  volumetric 
flow  rate  in  the  groove  was  negligible  for  groove  fill  ratios  of  less  than  <  0.5.  The 

Poiseuille  number  was  a  strong  function  of  the  countercurrent  shear  stress.  A  semi-analytical 
expression  was  provided  to  approximate  the  mean  velocity  as  a  function  of  the  shear  stress 
at  the  liquid-vapor  interface. 

3.8  Nomenclature 

Ai  cross-sectional  area  of  the  liquid,  m^ 

a;  Ai/h'^ 

Aiv  area  of  the  liquid-vapor  interface,  m^ 

Aw  area  of  the  groove  wall,  m^ 

Bo  Bond  number,  pgz\/cr 

d*  parameter  defined  in  eqn.  (48) 

Dh  hydraulic  diameter,  AAifP,  m 

Dt  Dhih 

f  friction  coefficient,  2r/pu^ 
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9 

h 

hi 

h*i 

n 

n* 

P 

P 

P* 

Plv 

p* 

Po 

R 

R* 

Re 

V 

V 
V* 

V* 

v' 

V 

V* 

w 

w* 

WI 

wt 

X,  y,  z 

X* 

Zi 


acceleration  due  to  gravity,  m/s^ 
groove  height,  m 

height  of  liquid  in  the  groove  at  the  wall,  m 

hi/h 

coordinate  normal  to  the  hquid-vapor  interface 
n/h 

pressure,  N/m^ 
wetted  perimeter,  m 
P/h 

perimeter  of  the  liquid-vapor  interface,  m 

Piv/h 

Poiseuille  number,  /Re 

radius  of  curvature  of  the  meniscus,  m 

R/h 

Reynolds  number,  pvDh/ix 
y-direction  velocity,  m/s 
average  y-direction  velocity,  m/s 
liv/h'^i-dp/dy) 

dimensionless  average  y-direction  velocity 
dimensionless  average  y-direction  velocity  when  Tj*  =  0 

V*/Vq 

volumetric  flow  rate,  vAi,  m^/s 

IxV/lkH-dp/dy)] 

period  of  the  sinusoidal  groove,  m 

w/h 

width  of  the  Hquid  in  the  groove,  m 
wi/h 

coordinate  directions 
x/h 

distance  from  the  liquid-vapor  interface  to  the  groove  bottom,  m 
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e 

fJ> 

P 

a 

'^Iv 

'^Iv 

^1*0,0 

T 

r' 

(t> 

^max 


z/h 

groove  aspect  ratio,  w/2h 
convergence  criteria 
absolute  viscosity,  Pa-s 
density,  kg/m^ 
surface  tension,  N/m 

shear  stress  at  the  Hquid- vapor  interface,  N/m^ 

Ti^/h{-dp/dy) 

dimensionless  shear  stress  at  the  liquid-vapor  interface  when  u*  =  0 
average  shear  stress ,  Ai{—dp /dy)  /P,  N/m^ 
average  shear  stress  at  the  groove  wall,  N/m^ 
r;;/h{-dp/dy) 

meniscus  contact  angle,  rad 
maximum  meniscus  contact  angle,  rad 
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3.10  Appendix 

3.10.1  Equation  of  the  circular  meniscus 


To  determine  the  equation  of  the  circular  meniscus,  a  relationship  must  be  made  based  on 
the  slope  of  the  tangent  Line  at  the  sinusoidal  groove  wall,  which  is  given  by 


tan^ 


dz* 

dx* 


(64) 
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The  equation  of  the  sinusoid  for  the  groove  wall  is  given  by 


^  =-|l  +  cos^-7r|-  +  ljJ|  (65) 

Differentiating  eqn.  (65)  with  respect  to  x*  gives 

Therefore,  substituting  eqn.  (66)  into  eqn.  (64)  and  evaluating  x*  =  wf/2  results  in  the 
equation  of  the  slope  of  the  tangent  line  at  the  wall  of  the  sinusoidal  groove. 

,  ,  dz*  TT  .  r  f  Wi  A 1 

The  slope  of  the  tangent  hne  on  the  circular  meniscus  is  given  by 


tan0  = 


circle 


Prom  Fig.  37,  the  angle  6  can  be  found  as  a  function  of  angles  a,  (f),  and  ip 


a  =  TT  —  Ip 


0  +  ^  +  a  =  7r 
6  =  ip  —  (p 


tan0  = 


Taking  the  tangent  of  both  sides  of  eqn.  (71)  and  using  the  trigonometric  addition  formula 
for  tangent  gives 

tan  9  =  tan(^  —  (p)  (72) 

,  ,  .  tana±tanj0 

tan(Q;  ±  P)  =  - - Ar  (73) 

1  q^tanatan/? 

tan^-tan0 

tan0  =  ; - -  (74) 

l  +  tan^tan0 

The  general  equation  for  a  circle  with  the  z*  axis  as  a  line  of  symmetry  is 

+  {z*  -  z*f  =  R*^  (75) 

Taking  the  derivative  of  eqn.  (75)  with  respect  to  x*  gives 

dz* 

2x*  +  2{z*-z:)^^=0  (76) 
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Rearranging  eqn.  (76)  to  solve  for  dz*/dx*  gives 


dz*  _  X* 
dx*  (z*  -  z*) 


(77) 


By  substituting  x*  —  w* /2  and  z*  =  hf  into  eqn.  (77)  and  then  substituting  the  result  into 
eqn.  (68)  gives  the  slope  of  the  tangent  Hne  on  the  circular  meniscus 


dz* 

dx* 


circle 


Combining  eqns.  (74)  and  (78)  gives 


wr/2 

-  ^o) 


=  tan0 


(78) 


Wj  /2  tan  ^  —  tan  (f) 
{hi  —  z*)  l  +  tan^tan^ 

Rearranging  eqn.  (79)  to  solve  for  z*  gives 


(79) 


<  =  h*i  + 


2 


(1  +  tan  Ip  tan  (p) 


tan  Ip  —  tan  <p 


(80) 


where  z^  is  the  height  of  the  center  of  the  circle  on  the  z*  axis.  By  letting  d  be  the  distance 
between  the  center  of  the  circular  liquid-vapor  interface  and  the  height  of  the  liquid,  it  can 
be  shown  that 


Combining  eqns.  (80)  and  (84)  gives 


(81) 

(82) 

(83) 

(84) 


(1  +  tan '0  tan  0) 


tan'0  —  tan^ 


(1  4-  tan '0  tan  0) 


2 


tan^  —  tan0 


(85) 


(86) 


84 


where 


_  d*  —  tan  4> 
2  1  +  d*  tan  0 


rf*  =  tan^  =  -sin  ^-  +  1 


2  \^]  {^  + dHaxKj)) 


-ii)’- 


d*  —  tan  0 
^*\  2 

M  (1  +  d*  tan  (f>y 


\2  J  (d*- tan  0)2 

^*2  V  fi  1  roi) 

^  -VTJ  L  (>“tan0>)J 

The  equation  of  the  radius  of  curvature  of  the  meniscus  in  terms  of  meniscus  contact  angle 


R*=r^]  1+ 


(1  +  d*tan0)^ 
(d*  —  tan  0)2) 


3.10.2  Determination  of  the  maximum  meniscus  contact  angle 

The  maximum  meniscus  contact  angle  can  be  foimd  when  R*  oo  in  eqn.  (92).  It  can  be 
seen  that  R*  oo  as  (d*  —  tan  0)  — >  0.  Therefore 

d*  —  tan  0max  =  0  (93) 


Substituting  eqn.  (93)  into  eqn.  (88)  for  d*  results  in 

—  sin  -TT  f  ^  ~  =  0 

Solving  for  0niax  gives  the  equation  for  the  maximum  meniscus  contact  angle. 
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3.10.3  Determination  of  the  hydraulic  diameter 

The  hydraulic  diameter  of  the  flow  field  must  be  determined  so  that  the  Poiseuille  number 
can  be  calculated.  To  find  D\  =  it  can  be  shown  from  Fig.  38  that 


j4i  +  /I2  +  -^3  —  2"^/* 

^  ^  ^  -  .43 


(96) 


(97) 


The  area  Az  can  be  found  by  integrating  the  equation  for  the  sinusoidal  groove  wall,  eqn. 
(65)  over  0  <  a;*  <  w*/2. 

r^r/2  I 


rw^/2  ^  ^ 

^  Jo  2\ 


yis  =  — r  —  —  Sin 


1  +  cos 


P  . 


-X 


|+i)]}d.- 


fwf 


(98) 


(99) 


4  27r 

From  Fig.  38,  Ai  is  the  area  above  the  meniscus.  It  can  be  shown  that  ^41  is  one-half  of  the 
area  of  the  segment 


4 

=  —  cos 


-1 


R* 


dVR*^  -  dP 


where 


Thus,  the  area  Ai  is 


d=\  R*^ 


w. 


(100) 


(101) 


^1  = 


R* 


R*  cos 


-1 


'1- 


Wi 

m* 


2R> 


(102) 


Substituting  eqns.  (99)  and  (102)  into  eqn.  (97)  gives  the  cross-sectional  area  of  the  flow 
field 


wj 


A;=^{2h^-1)-R^ 


R*  cos 


-1 


T- 


wj 


■ 

^ —  sm 

X 


wl 


2R* 
for  (j)  <C  (ptnai 


^\di- 


2W 


(103) 


When  the  meniscus  contact  angle  approaches  the  maximum  value  0  0max,  the  area  Ai 
0.  Therefore,  the  cross-sectional  area  of  the  flow  field  for  this  case  is 


a;  =  w^h*i  -  2Az 


(104) 
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(105) 


The  equation  for  the  wetted  perimeter  has  no  closed-form  solution,  and  consequently  the 
equation  for  the  hydraulic  diameter  has  no  closed-form  solution.  The  reason  is  that  the 
integral  in  eqn.  (109)  is  elUptic,  and  must  be  solved  numerically. 
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Figure  30:  Grooves  chemically  etched  in  glass  (Courtesy  of  D.  Liepmann,  University  of 
California  at  Berkeley). 


88 


89 


-0.00193 


Figure  32:  Dimensionless  velocity  fields  for  laminar  flow  in  a  sinusoidal  groove  (/5  = 
w*il2  =  0.25,  (f)  =  25°):  (a)  =  2.0  (cocurrent  flow);  (6)  r^*  =  0.0;  (c)  = 

(countercurrent  flow). 


:  0.5, 
-0.1 


10 


Figure  33:  Variation  of  the  flow  variables  with  shear  stress  at  the  liquid-vapor  interface  for 
various  values  of  meniscus  contact  angle  {(3  =  0.5,  wf/2  =  0.25,  P*  =  1.15245):  (a)  Mean 
velocity;  (6)  Volumetric  flow  rate;  (c)  PoiseuiUe  number. 


Figure  34:  Variation  of  the  flow  variables  with  groove  fill  ratio  for  various  values  of  groove 
aspect  ratio  (r^*  =  0.0,  (f)  =  0°):  (a)  Mean  velocity;  (6)  Volumetric  flow  rate;  (c)  Poiseuille 
number. 


Figure  36:  Comparison  of  the  semi-analytical  solution  with  numerical  data  [P  =  0.5,  wlft  = 
0.25):  (o)  Countercurrent  vapor  shear  stress  required  for  v*  —  0;  (6)  Normalized  mean 
velocity  versus  normalized  shear  stress  at  the  liquid-vapor  interface. 
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Figure  37:  Definition  of  geometric  parameters. 


95 


Figure  38:  Determination  of  the  cross-sectional  area  of  the  Hquid. 
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4  Micro  Capillary  Pumped  Loop  Testing 


4.1  Abstract 

This  project  seeks  to  determine  the  operational  characteristics  of  a  micro-capillary  pumped 
loop.  This  thermal  management  device  is  to  be  directly  integrated  into  the  silicon  chips  used 
in  various  devices.  An  experimental  testbed  has  been  constructed  and  calibrated  while  the 
test  specimens  were  being  manufactured.  This  chapter  provides  details  concerning  the  ex¬ 
perimental  setup  and  preliminary  cahbration  information,  which  will  be  used  in  the  near 
future  when  a  fully  instrumented  micro-capillary  pumped  loop  will  be  tested. 

4.2  Experimental  Setup 

The  objective  of  these  experiments  was  to  determine  the  heat  transfer  capabilities  of 
micro-Capillary  Pumped  Loops  (/zCPLs)  which  are  formed  using  a  layer  of  borosihcate  glass 
anodically  bonded  to  a  silicon  wafer.  These  nCFLs  are  being  tested  to  determine  their 
suitability  for  the  thermal  management  of  microprocessors. 

Figure  39  shows  two  typical  yttCPL  designs.  Several  /iCPLs  were  manufactured  on  a  3.9- 
inch  dia.  sihcon  wafer  at  the  University  of  California  at  Berkeley  using  photolithographic 
techniques.  The  individual  /LfCPLs  were  cut  from  the  wafer  using  a  diamond  saw  on  a 
mini  vertical  mill.  The  /iCPLs  consisted  of  a  series  of  channels  etched  in  the  silicon  wafer. 
These  channels  made  up  part  of  the  evaporator  section,  the  condenser  section,  and  the 
interconnecting  vapor  and  liquid  hues.  In  addition,  a  hole  was  etched  completely  through 
the  sihcon  to  allow  the  installation  of  a  fill  tube,  which  would  in  turn  be  attached  to  a 
remote  working  fluid  reservoir.  A  thin  sheet  of  borosihcate  glass  was  bonded  to  the  front 
of  the  silicon.  This  glass  had  grooves  etched  into  the  sections  that  were  positioned  over  the 
evaporator  and  condenser  sections  of  the  fiCFL.  These  grooves  acted  as  a  wick  structure 
and  provided  the  capillary  pressure  that  was  required  for  circulation  of  the  working  fluid. 
In  an  actual  application,  the  glass  would  be  replaced  by  another  sheet  of  sihcon,  but  for 
experimental  purposes,  glass  was  used  to  allow  observation  of  the  working  fluid  while  the 
/iCPL  was  being  tested.  Another  reason  that  glass  was  used  was  that,  while  sihcon  is 
virtuaUy  transparent  to  the  CO2  laser  energy,  the  borosihcate  glass  absorbed  a  significant 
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portion  of  the  incident  laser  light.  This  heated  the  glass  and  provided  the  source  of  thermal 
energy  needed  for  testing  the  //CPLs. 

Figure  40  shows  an  overall  schematic  diagram  of  the  experimental  setup.  This  setup  pro¬ 
vided  a  means  to  input  heat  to  the  //CPL  evaporator  section  using  a  laser,  extract  heat  from 
the  condenser  section  using  a  calorimeter,  insulate  the  /iCPL  from  parasitic  heat  losses,  and 
record  temperature  data.  In  addition,  photographic  records  of  the  tests  could  be  collected. 

The  juCPL  wafer  was  enclosed  in  a  stainless  steel  vacuum  chamber  to  reduce  natural 
convective  heat  losses  from  the  surface.  The  chamber  could  be  evacuated  to  a  pressure 
of  less  than  10“^  Torr  (Varian  Turbo-V450  vacuum  pump),  which  was  measured  using  a 
thermocouple  gauge  and  a  more  accurate  ionization  gauge  (Varian  Model  0571-K2471-303). 
In  addition,  multi-layer  mylar  radiation  insulation  was  installed  on  the  interior  of  the  vacuum 
chamber  to  reduce  thermal  radiation  heat  transfer.  The  /iCPL  was  isolated  from  conductive 
heat  transfer  by  using  an  acryhc  holding  fixture,  as  shown  in  Fig.  41.  Grooves  were  machined 
in  the  upper  and  lower  acrylic  fixtures  to  hold  the  juCPL.  The  fixtures  were  in  turn  mounted 
to  an  aluminum  ring,  which  was  held  in  the  chuck  of  a  micrometer  stage.  This  stage  could 
rotate  the  aluminum  ring,  as  well  as  move  it  horizontally  and  vertically,  to  position  the  ^CPL 
in  front  of  the  ZnSe  window  in  the  end  of  the  vacuum  chamber. 

Heat  was  applied  to  the  /iCPL  evaporator  section  by  a  25- Watt  CO2  laser  (Synrad  Model 
J48-2W).  Since  the  wavelength  of  this  laser  (10.5  fim)  is  not  in  the  visible  range,  a  Helium- 
Neon  laser  (Uniphase  Novette  1508-0)  was  used  to  determine  the  location  of  the  CO2  laser 
spot.  Ahgnment  was  verified  by  placing  a  piece  of  paper  in  front  of  the  vacuum  chamber  and 
burning  it  slightly  with  the  CO2  laser.  After  the  CO2  laser  was  deactivated,  the  removable 
mirror  Ml  (Fig.  40)  was  then  installed,  and  the  HeNe  was  activated.  The  HeNe  laser  beam 
was  then  afigned  with  the  biun  spot  using  adjusting  screws  on  the  removable  mirror.  Once 
the  two  beams  were  aligned,  changing  the  angle  of  the  turning  mirror  M2  would  redirect  both 
laser  beams.  In  this  way,  the  CO2  laser  beam  could  be  aimed  onto  the  fxCPL  evaporator 
using  the  HeNe  and  the  adjusting  screws  on  the  turning  mirror  M2.  The  CO2  laser  beam 
was  turned,  focused  with  a  plano-convex  lens  (ZnSe,  1.5  inch  dia.,  30  ft  focal  length,  AR 
coated  for  9.2  -  10.6  jum),  and  passed  through  a  window  (ZnSe,  AR  coated  for  9.2  -  10.6 
lj,m)  into  the  vacuum  chamber.  The  focusing  lens  provided  a  means  to  control  the  diameter 
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of  the  beam  at  the  surface  of  the  test  specimen,  and  therefore  the  heat  flux  applied  to  the 
evaporator  of  the  yiiCPL.  The  CO2  laser  beam  was  directed  through  a  length  of  PVC  pipe  to 
reduce  the  possibility  of  foreign  objects  accidentally  being  placed  in  the  path  of  the  invisible 
beam,  which  would  have  caused  a  hazardous  situation.  There  were  also  beam  stops  made  of 
plexiglas,  which  is  opaque  to  CO2  laser  hght,  along  the  sides  of  the  beam  path. 

A  water-cooled  calorimeter  was  designed  and  manufactured  to  determine  the  heat  rate 
transported  by  the  /liCPL  (Fig.  42).  It  consisted  of  a  small  block  of  copper  with  stainless  steel 
tubing  silver-soldered  into  holes  which  had  been  drilled  into  one  face  of  the  block.  A  hole  was 
drilled  into  another  face  of  the  block,  which  simultaneously  opened  both  tubes.  This  face 
was  then  adhered  to  the  fiCPL  using  RTV  sealant  such  that  the  coolant  was  in  direct  contact 
with  the  glass  or  silicon.  The  vacuum  chamber  had  pass-through  fittings  to  allow  the  coolant 
to  flow  to  and  from  the  calorimeter.  Tygon  tubing  was  used  in  the  interior  of  the  vacuum 
chamber  to  connect  the  calorimeter  to  the  pass-through  fittings.  Type  E  thermocouples 
were  placed  in  the  coolant  inlet  and  outlet  lines  within  the  vacuum  chamber  as  close  to  the 
calorimeter  as  was  practical.  In  this  way,  convective  and  conductive  heat  losses  from  the 
tygon  tubing  and  along  the  length  of  the  thermocouple  wires  to  the  ambient  were  mitigated. 
As  shown  in  Fig.  40,  an  open  reservoir  was  used  to  deliver  coolant  (water)  at  a  constant 
mass  flow  rate  to  the  calorimeter.  The  coolant  line  from  the  reservoir  to  the  calorimeter  was 
heavily  insulated  to  reduce  temperature  swings  when  the  room  air  conditioner  cycled  on  and 
off.  The  mass  flow  rate  of  the  water  flowing  through  the  calorimeter  was  calculated  using  a 
stopwatch  and  a  Mettler  scale  (Model  PC  4400). 

Three  electrical  feedthroughs  were  mounted  on  the  top  of  the  chamber.  One  feedthrough 
supplied  electrical  power  to  the  interior  of  the  chamber.  The  other  two  were  five-pair,  Type 
E  (Chromel/Constantan)  thermocouple  feedthroughs.  The  thermocouples  that  were  used  to 
make  temperature  measurements  directly  on  the  different  points  of  the  /iCPL  were  attached 
to  the  backside  of  the  silicon  using  Kapton  Tape  or  vacuum  epoxy,  as  shown  in  Fig.  43. 

A  digital-capable  video  camera  (Panasonic  GP-KR222)  with  a  long-range  microscope  lens 
and  a  monitor  with  a  VCR  were  used  to  observe  the  /itCPLs  while  they  were  being  tested. 
The  camera  was  suspended  just  above  the  path  of  the  CO2  laser  beam  to  reduce  parallax. 
Image-  and  video-capturing  software  was  used  to  store  pictures  and  videos  directly  to  a  PC 
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during  transient  events  such  as  start-up,  boiling  and  dry-out. 


4.2.1  Calibration 

In  order  to  accurately  test  the  /.tCPLs,  the  amount  of  laser  power  delivered  to  the  wafer  had 
to  be  quantified.  The  thermocouples  used  for  sensing  the  temperature  of  the  wafer  and  the 
coolant  inlet  and  outlet  had  to  be  calibrated.  This  section  gives  a  detailed  description  of  the 
procedures  used  to  achieve  these  calibrations. 

Laser  Controller  Calibration 

The  CO2  laser  beam  was  attenuated  and  partially  reflected  as  it  passed  through  the  vac¬ 
uum  chamber  window,  and  partially  reflected  when  it  struck  the  //CPL  wafer.  Therefore,  it 
was  necessary  to  calibrate  to  determine  the  above-mentioned  losses.  Power  measurements 
were  taken  in  front  of  and  behind  the  vacuum  chamber  window.  This  provided  a  mea¬ 
sure  of  the  amount  of  attenuation/reflection  caused  by  the  window.  The  power  output  of 
the  laser  was  measured  by  a  calibrated  optical  power  meter/disk  thermopile  combination 
(Newport  Models  1815-C  and  818T-10).  A  small  blank  wafer  of  the  glass/sificon  sandwich, 
instrumented  with  a  water-cooled  calorimeter  inside  the  evacuated  chamber,  was  used  to 
determine  the  amount  of  power  actually  absorbed  by  the  wafer.  The  amount  of  laser  power 
incident  on  the  window  and  the  wafer,  as  well  as  the  amount  of  power  extracted  by  the 
calorimeter,  were  calibrated  to  the  laser  controller  setting  (LCS),  which  was  in  terms  of  the 
percentage  of  the  full  power  output  of  the  laser.  In  a  preliminary  setup,  a  ZnSe  beam  splitter 
was  used  with  the  power  meter  to  directly  measure  the  laser  power  output.  However,  it  was 
determined  that  this  arrangement  caused  the  amount  of  laser  power  delivered  to  the  wafer 
to  fluctuate  significantly,  due  to  the  dependence  of  the  reflectivity  of  the  beam  splitter  on 
the  wavelength  of  the  laser.  Line  hopping,  or  the  variation  of  the  wavelength  of  the  CO2 
laser  with  respect  to  time,  could  not  be  avoided,  so  the  beam  splitter  was  removed  from  the 
system. 

The  following  procedure  was  used  to  obtain  the  power  readings  in  front  of  and  behind 
the  chamber: 

1.  Turn  on  the  recirculating  chiller,  external  pump  and  power  supply  for  the  CO2  laser. 
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2.  Using  the  HeNe  laser,  position  the  power  meter  in  front  of  the  chamber  ( “fore”  posi¬ 
tion). 

3.  Remove  the  Ml  mirror  from  the  kinematic  mount. 

4.  Set  the  laser  controller  setting  to  1.0%. 

5.  Turn  on  the  CO2  laser. 

6.  After  a  short  warm-up  period,  collect  500  data  points  using  the  power  meter  in  front 
of  the  chamber. 

7.  Interrupt  the  beam  with  the  laser  aperture  shutter  on  top  of  the  laser. 

8.  Place  the  power  meter  behind  the  chamber  (“aft”  position). 

9.  Open  the  shutter  on  the  laser. 

10.  Collect  500  data  points  using  the  power  meter  behind  the  chamber. 

11.  Using  this  method,  alternate  between  the  fore  and  aft  positions  until  four  calibration 
points  (500  data  points  each)  have  been  collected  both  in  front  of  and  behind  the 
chamber. 

12.  Repeat  steps  4  through  11,  incrementing  up  through  the  desired  controller  settings. 

13.  Using  the  same  procedure,  start  at  the  maximum  LCS  value  and  work  back  down 
through  the  specified  settings.  This  will  produce  eight  calibration  points  at  each  loca¬ 
tion  (fore  and  aft)  for  each  laser  controller  setting. 

The  procedure  for  determining  the  amount  of  power  extracted  by  the  water-cooled 
calorimeter  on  the  blank  wafer  is  as  follows: 

1.  Turn  on  the  recirculating  chiller,  external  pump  and  power  supply  for  the  CO2  laser. 

2.  Position  the  power  meter  in  front  of  the  vacuum  chamber  to  validate  the  laser  power 
setting. 
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3.  Set  the  laser  controller  setting  to  1.0%. 

4.  Turn  on  the  CO2  laser. 

5.  After  a  short  warm-up  period,  collect  500  data  points  using  the  power  meter  in  front 
of  the  chamber. 

6.  Interrupt  the  beam  with  the  laser  aperture  shutter  on  top  of  the  laser. 

7.  Remove  the  power  meter. 

8.  Start  the  coolant  flow  rate  to  the  calorimeter. 

9.  Open  the  shutter  on  the  laser. 

10.  Adjust  the  coolant  flow  rate  so  the  temperature  increase  across  the  calorimeter  is 
approximately  2°C. 

11.  Allow  the  calorimeter  temperatures  to  reach  steady  state  (approximately  30  min.). 

12.  Take  8  sets  of  600  data  points  for  the  calorimeter  inlet  and  outlet  temperatures.  For 
each  set,  measure  the  mass  flow  rate  through  the  calorimeter. 

13.  Reinstall  the  power  meter  in  front  of  the  vacuum  chamber  and  collect  500  data  points 
for  the  laser  power  at  that  location. 

14.  Repeat  steps  3  through  13,  incrementing  up  through  the  desired  controller  settings. 

Figures  44(a)  and  44(6)  show  typical  traces  obtained  during  calibration  for  the  calorimeter 
temperature  difference  and  the  laser  power  versus  time,  respectively.  In  addition.  Fig.  44(c) 
shows  the  results  of  the  calibration  with  respect  to  the  laser  controller  setting.  It  should 
be  noted  that  the  borosilicate  glass  melted  at  LCS  =  30%,  as  indicated  by  the  brackets 
surrounding  this  data.  Tables  9 — 11  provide  the  actual  values  presented  in  Fig.  44(c).  The 
hneax  best-fit  equations  for  the  data  in  Fig.  44(c)  are  as  follows: 


^window  =  0.2874  *  LCS  {R^  =  0.9997) 


(110) 


Qout  =  0.2011  *  LCS  =  0.9974)  (112) 

The  linear  equation  for  the  power  absorbed  by  the  wafer  does  not  include  the  data  at  LCS 
=  30%.  In  comparing  eqns.  (Ill)  and  (112),  the  amount  of  incident  laser  power  absorbed 
by  the  borosilicate  glass  is  approximately  78%. 

Thermocouple  Calibration 

The  thermocouples  used  for  wall  temperature  measurements  and  calorimetry  were  cah- 
brated  using  a  recirculating  chiller  (Lauda  Model  RC20)  with  a  serial  interface  (Lauda  Model 
R61)  and  an  RTD  (Hart  Scientific  model  5691).  The  working  bath  of  the  chiller  was  filled 
with  PAO  oil  to  allow  for  more  stable  temperatures.  In  the  bath,  there  was  a  length  of 
copper  pipe  with  a  flat  sheet  of  copper  silver-soldered  over  one  end.  Copper  was  chosen  for 
its  high  thermal  conductivity.  This  sheet  of  copper  acted  as  a  cap  and  a  base.  The  other 
end  of  the  copper  pipe  was  left  open.  This  allowed  it  to  be  filled  with  water  and  immersed 
in  the  oil  bath,  keeping  the  open  end  of  the  tube  above  the  surface  of  the  oil.  The  RTD 
and  thermocouples  were  placed  in  the  water  in  this  copper  tube,  to  protect  them  from  the 
oil  in  the  bath.  The  thermocouples  were  attached  to  the  RTD,  and  the  RTD  was  supported 
using  a  lab  clamp  to  insure  that  it  did  not  contact  the  sides  or  the  base  of  the  copper  tube. 
The  target  temperature  of  the  bath  was  set  by  a  PC,  which  also  read  the  RTD  temperature 
and  received  information  on  the  thermocouple  temperatures  from  a  data  acquisition  unit 
(Hewlett  Packard  Model  HP34970A).  The  data  acquisition  and  temperature  settings  were 
controlled  by  an  Excel/ Visual  Basic  for  Applications  program.  A  test  program  was  written 
to  determine  the  length  of  time  necessary  for  the  RTD  and  thermocouple  readings  inside  the 
water-filled  copper  tube  to  reach  steady  state.  The  temperature  of  the  bath  was  set  to  38°C, 
and  allowed  to  stabilize  for  30  min.  The  program  then  set  the  target  bath  temperature  to 
40°C,  and  began  taking  data  from  the  bath  temperature  sensor,  the  RTD,  and  a  thermocou¬ 
ple.  Scans  were  recorded  continuously  for  approximately  2  hr  as  shown  in  Fig.  45.  It  was 
determined  that  the  readings  of  the  RTD  and  the  thermocouple  in  the  water  reached  steady 
state  after  approximately  30  min. 

After  the  time  taken  to  reach  steady  state  was  known,  the  actual  calibration  tests  could 
be  performed.  A  new  program  was  written  which  initially  set  the  bath  temperature  to  10°C, 
waited  for  60  min.,  and  then  took  100  data  points  of  5  scans  each.  The  60-min.  interval  was 
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decided  on  because  it  should  have  given  more  than  enough  time  for  the  readings  to  reach 
steady  state.  The  program  then  incremented  the  temperature  to  15°C,  waited  one  hr,  and 
recorded  data  again.  It  continued  in  this  manner  from  10°C  to  95°C,  in  5°C  increments. 
After  it  had  been  set  at  95°C  for  one  hr  and  had  taken  the  data  for  that  temperature,  it 
then  went  down  to  92.5°C.  From  there,  it  went  down  iri  the  same  manner  (5°C  intervals 
with  60-min.  time  intervals)  to  12.5°C.  When  the  program  was  finished,  there  were  500  data 
points  for  each  temperature  from  10°C  to  95°C  in  increments  of  2.5°C.  This  raw  data  was 
used  to  determine  the  calibration  curves  for  the  twelve  thermocouples.  The  temperature 
reading  of  each  thermocouple  was  plotted  against  the  reading  of  the  RTD.  A  linear  best-fit 
line  was  determined  in  the  form  Tactual  =  'mTieading  +  b,  as  shown  in  Fig.  46(a).  Next,  all  of 
the  temperature  measurements  for  a  specific  thermocouple  were  evaluated  by  its  respective 
best-fit  equation.  The  difference  between  the  value  from  the  equation  and  the  corresponding 
RTD  value  was  calculated  as  shown  in  Fig.  46(6).  This  difference  plus  the  stated  NIST- 
traceable  uncertainty  of  the  RTD  (±0.005°C)  was  the  final  uncertainty  for  any  subsequent 
temperature  measurements  for  a  given  thermocouple  (Table  12).  Thermocouples  211  and 
212  were  used  in  the  calorimeter.  Since  the  coolant  temperature  was  not  anticipated  to  rise 
above  40°C,  only  the  calibration  data  between  10°C  and  40°C  were  used  when  calibrating 
these  two  thermocouples.  There  were  some  outlying  data  points  that  were  disregarded  when 
reporting  the  maximum  errors  in  Table  12.  There  were  no  more  than  10  out  of  3500  of  these 
points  dropped  from  each  graph. 


4.2.2  Uncertainty  Analysis 


The  formula  used  for  finding  the  mass  flow  rate  of  coolant  flowing  through  the  calorimeter 
is 

(113) 

V 

The  root-sum-square  uncertainty  for  the  mass  flow  rate  is 
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(114) 
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The  amount  of  heat  removed  from  the  {iCPL  by  the  calorimeter  can  be  expressed  by 


Qout  =  mC'p(To„t-7k) 


(115) 


The  uncertainty  for  values  found  using  this  equation  is 


The  thermal  resistances  based  on  incident  laser  power  and  power  extracted  by  the 
calorimeter,  respectively,  are 

(117) 

Vin 


(120) 

4.3  Summary  and  Future  Work 

This  report  documents  the  status  of  the  project  as  of  the  end  of  Summer  2001.  A  great 
deal  of  the  experimental  setup  is  complete  and  has  been  thoroughly  cahbrated.  A  new  gener¬ 
ation  of  /.iCPLs  have  been  fabricated,  and  are  in  the  process  of  being  diced  and  instrumented. 
A  major  problem  in  the  past  has  been  making  sure  that  the  juCPL  is  hermetically  sealed.  To 
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address  this  problem,  the  vacuum-tightness  of  the  fiCFL  will  be  evaluated  using  a  helium 
leak  detection  unit.  The  /liCPL  will  then  be  charged  and  tested  in  the  vacuum  chamber  to 
determine  the  operating  characteristics  for  various  power  input  levels.  Problems  are  antic¬ 
ipated  with  respect  to  the  start-up  behavior.  The  interaction  between  the  jnCPL  and  the 
working  fluid  reservoir  is  imknown  at  present,  but  will  be  controllable  with  an  electrical 
heater  tape  mounted  to  the  reservoir. 
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Table  9:  Laser  calibration  data. 


LCS  (%) 

^window  i  Cr  (W) 

1 

0.295  ±  0.0034 

0.265  ±  0.011 

1 

0.292  ±  0.0025 

0.261  ±  0.0080 

1 

0.293  ±  0.0044 

0.253  ±  0.0026 

1 

0.290  ±  0.0063 

0.256  ±  0.0043 

1 

0.290  ±  0.0010 

1 

0.291  ±  0.0042 

0.271  ±  0.0043 

1 

0.294  ±  0.0023 

0.268  ±  0.0019 

1 

0.307  ±  0.0028 

0.287  ±  0.0036 

1.5 

0.411  ±  0.0077 

0.360  ±  0.012 

1:5 

0.408  ±  0.0060 

0.356  ±  0.0058 

1.5 

0.409  ±  0.0054 

0.362  ±  0.014 

1.5 

0.409  ±  0.0053 

0.358  ±  0.0099 

1.5 

0.408  ±  0.0012 

0.380  ±  0.0057 

1.5 

0.404  ±  0.0079 

0.380  ±  0.0081 

1.5 

0.396  ±  0.0019 

0.376  ±  0.0071 

1.5 

0.397  ±  0.0023 

0.370  ±  0.0057 

2 

0.529  ±  0.0011 

0.489  ±  0.0060 

2 

0.526  ±  0.0037 

0.490  ±  0.0041 

2 

0.531  ±  0.0047 

0.489  ±  0.0038 

2 

0.532  ±  0.0045 

0.489  ±  0.0049 

2 

0.545  ±  0.0052 

0.502  ±  0.0093 

2 

0.549  ±  0.0016 

0.495  ±  0.0084 

2 

0.547  ±  0.0025 

0.484  ±  0.0037 

2 

0.547  ±  0.0035 

0.479  ±  0.0029 

3 

0.894  ±  0.0089 

0.777  ±  0.0033 

3 

0.892  ±  0.0055 

0.795  ±  0.015 

3 

0.890  ±  0.0032 

0.795  ±  0.016 

3 

0.890  ±  0.0014 

0.810  ±  0.012 

3 

0.915  ±  0.0031 

0.811  ±  0.0049 

3 

0.921  ±  0.0027 

0.812  ±  0.0059 

3 

0.922  ±  0.0032 

0.809  ±  0.0032 

3 

0.921  ±  0.0037 

0.807  ±  0.0031 

4 

1.15  ±  0.0079 

1.05  ±  0.0054 

4 

1.14  ±  0.0086 

1.01  ±  0.013 

4 

1.14  ±  0.0082 

1.01  ±  0.0056 

4 

1.14  ±  0.012 

1.02  ±  0.019 

4 

1.14  ±  0.0018 

1.06  ±  0.0017 

4 

1.15  ±  0.0061 

1.06  ±  0.0015 

4 

1.16  ±  0.011 

1.08  ±  0.0040 

4 

1.17  ±  0.0027 

1.08  ±  0.0064 

0.123 

± 

9% 

0.120 

± 

9% 

0.126 

± 

8% 

0.124 

± 

9% 

0.123 

± 

7% 

0.122 

± 

8% 

0.118 

± 

7% 

0.113 

± 

8% 

0.186 

± 

8% 

0.183 

± 

6% 

0.177 

± 

7% 

0.177 

± 

6% 

0.173 

± 

6% 

0.173 

± 

6% 

0.172 

± 

6% 

0.166 

± 

6% 

0.303 

± 

8% 

0.305 

± 

7% 

0.310 

± 

7% 

0.309 

± 

7% 

0.315 

± 

7% 

0.307 

± 

7% 

0.316 

± 

7% 

0.323 

± 

6% 

0.580 

± 

8% 

0.593 

± 

8% 

0.585 

± 

9% 

0.589 

± 

7% 

0.586 

± 

9% 

0.588 

± 

7% 

0.585 

± 

8% 

0.594 

± 

8% 

0.811 

± 

8% 

0.796 

± 

8% 

0.811 

± 

8% 

0.790 

± 

8% 

0.803 

± 

8% 

0.808 

± 

9% 

0.789 

± 

8% 

0.820 

± 

9% 

Table  10:  Laser  calibration  data,  continued. 


LCS  (%) 

5 

1.49  ±  0.015 

1.33  ±  0.016 

1.01  ±  6% 

5 

1.48  ±  0.0031 

1.32  ±  0.0024 

1.01  ±  6% 

5 

1.48  ±  0.0030 

1.34  ±  0.015 

1.03  ±  7% 

5 

1.49  ±  0.013 

1.36  ±  0.0058 

1.02  ±  5% 

5 

1.50  ±  0.0024 

1.39  ±  0.0018 

1.01  ±  5% 

5 

1.49  ±  0.0066 

1.38  ±  0.0059 

1.03  ±  5% 

5 

1.49  ±  0.0045 

1.38  ±  0.0064 

1.02  ±  5% 

5 

1.49  ±  0.0018 

1.38  ±  0.0034 

1.02  ±  5% 

10 

2.96  ±  0.010 

2.66  ±  0.036 

2.08  ±  4% 

10 

2.97  ±  0.016 

2.64  ±  0.024 

2.07  ±  3% 

2.97  ±  0.012 

2.72  ±  0.039 

2.09  ±  4% 

2.97  ±  0.0045 

2.67  ±  0.042 

2.08  ±  4% 

3.02  ±  0.044 

2.72  ±  0.0042 

2.06  ±  3% 

10 

2.97  ±  0.0055 

2.73  ±  0.0035 

2.10  ±  4% 

10 

2.96  ±  0.0042 

2.73  ±  0.0045 

2.04  ±  4% 

10 

2.96  ±  0.0026 

2.72  ±  0.0040 

2.03  ±  6% 

15 

4.31  ±  0.024 

3.90  ±  0.057 

3.07  ±  5% 

15 

4.30  ±  0.014 

3.91  ±  0.0046 

3.12  ±  5% 

15 

4.32  ±  0.030 

3.91  ±  0.0040 

3.10  ±  6% 

15 

4.31  ±  0.028 

3.83  ±  0.083 

3.06  ±  5% 

15 

4.37  ±  0.032 

3.92  ±  0.069 

3.00  ±  8% 

15 

4.35  ±  0.016 

3.98  ±  0.0065 

3.11  ±  6% 

4.35  ±  0.017 

3.98  ±  0.0047 

3.12  ±  4% 

4.35  ±  0.012 

3.98  ±  0.0085 

3.06  ±  5% 

■RiH 

5.75  ±  0.076 

5.18  ±  0.11 

4.05  ±  6% 

20 

5.70  ±  0.0051 

5.22  ±  0.14 

4.05  ±  7% 

20 

5.75  ±  0.093 

5.10  ±  0.0032 

4.08  ±  7% 

20 

5.76  ±  0.10 

5.16  ±  0.11 

4.15  ±  7% 

20 

5.74  ±  0.0038 

5.22  ±  0.0059 

4.16  ±  7% 

20 

5.73  ±  0.0046 

5.22  ±  0.0040 

4.16  ±  7% 

20 

5.73  ±  0.0057 

5.22  ±  0.0038 

4.10  ±  6% 

20 

5.73  ±  0.0049 

5.22  ±  0.0056 

4.18  ±  8% 
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Table  11:  Laser  calibration  data,  continued. 


LCS  (%) 


7.08  ±  0.12 


7.06  ±  0.0063 


7.07  ±  0.0058 


7.10  ±  0.039 


7.10  ±  0.040 


7.09  ±  0.020 


7.11  ±  0.039 


7.15  ±  0.061 


8.65  ±  0.028 


8.67  ±  0.029 


8.68  ±  0.038 


8.69  ±  0.035 


8.62  ±  0.052 


8.63  ±  0.036 


8.61  ±  0.0050 


8.64  ±  0.034 


0.306  ±  0.0015 


0.312  ±  0.0011 


0.406  ±  0.0057 


0.429  ±  0.0011 


0.549  ±  0.0012 


0.546  ±  0.0012 


0.903  ±  0.0030 


0.909  ±  0.0071 


1.17  ±  0.0032 


1.17  ±  0.0030 


1.52  ±  0.0021 


1.50  ±  0.0024 


3.00  ±  0.0074 


3.02  ±  0.0082 


4.35  ±  0.0048 


4.34  ±  0.0050 


5.81  ±  0.0066 


5.77  ±  0.010 


7.14  ±  0.044 


7.14  ±  0.039 


6.36  ±  0.14 


6.26  ±  0.039 


6.29  ±  0.014 


6.28  ±  0.0068 


6.32  ±  0.052 


6.30  ±  0.035 


6.40  ±  0.0089 


6.41  ±  0.0052 


7.68  ±  0.0091 


7.78  ±  0.099 


7.88  ±  0.025 


7.86  ±  0.0078 


7.67  ±  0.023 


7.67  ±  0.028 


7.69  ±  0.048 


7.72  ±  0.060 


4.88  ±  6% 


4.81  ±  7% 


4.88  ±  8% 


4.91  ±  5% 


4.99  ±  6% 


4.85  ±  7% 


4.92  ±  6% 


5.02  ±  6% 


6.13  ±  7% 


6.23  ±  6% 


6.30  ±  8% 


6.16  ±  8% 


6.30  ±  5% 


6.21  ±  6% 


6.27  ±  6% 


6.14  ±  6% 
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Table  12:  Thermocouple  calibration  equations:  Tkctuai  =  niTreading  +  b  (Valid  ranges:  TCs 
201-210,  10°C  <  Treading  <  95°C;  TCs  211  and  212,  10°C  <  Treading  <  40°C). 


TC  Number 

m 

6(°C) 

Max  Error  (°C) 

201 

0.9952 

1.2794 

0.30 

202 

0.9962 

1.0023 

■BiHi 

0.9956 

0.9031 

0.9952 

0.6704 

0.30 

205 

0.9958 

0.4097 

0.30 

206 

0.9965 

0.1150 

0.25 

207 

0.9935 

0.1410 

0.40 

208 

0.9968 

0.0079 

IIIIIIIIIIIIIIQQIII^^ 

209 

0.9965 

-0.0507 

210 

0.9964 

0.0122 

0.15 

211 

0.9936 

0.1318 

-0.093 

212 

0.9947 

-0.0919 

-0.054 

4.4  Nomenclature 


Cp 

LCS 

m 

m 

Q 

Qin 

Qout 

Q-windaw 

i'^th.out 

t 

Tin 

Tout 

T 

v,e 

T 

v  .r 


specific  heat,  J/(kg-K) 
laser  controller  setting,  % 
mass,  kg 

mass  flow  rate,  kg/s 

rate  of  heat  transfer,  W 

laser  power  incident  on  the  blank  wafer,  W 

power  removed  by  the  calorimeter  mounted  to  the  blank  wafer,  W 
laser  power  incident  on  the  window  of  the  vacuum  chamber,  W 
thermal  resistance  based  on  incident  power,  K/W 
thermal  resistance  based  on  power  removed  by  the  calorimeter,  K/W 
time,  s 

calorimeter  inlet  temperature,  K 
calorimeter  outlet  temperature,  K 
vapor  fine  temperature  at  the  evaporator,  K 
vapor  fine  temperature  at  the  condenser,  K 


no 


Serpentine 

Condenser 


Grooved 

Condenser 


Figure  39:  Variations  of  the  basic  /.iCPL  desigTi. 
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Figure  40:  Schematic  of  ^-CPL  test  setup. 


Figure  41:  Moimting  fixture  for  the  //CPL. 
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Figiire  42:  Photograph  of  the  copper  calorimeter. 


BEST  AVAILABLE  COPY 
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Figure  43:  Photograph  showing  thermocouples  adhered  to  the  back  of  the  /zCPL  wafer  with 
vacuum  epoxy.  The  stainless  steel  fill  tube  is  on  the  left-hand  side. 


BEST  AVAILABLE  COPY 


Time  (s) 


Laser  Controller  Setting  (%) 


Figure  44;  CO2  laser  calibration;  (o)  Typical  calorimeter  temperature  increase  trace  (LCS 
=  10%);  (6)  Typical  laser  power  measurement  traces  (LCS  =  10%);  (c)  Heat  rates  versus 
LCS  (Data  points  in  square  brackets  were  taken  when  the  borosilicate  glass  melted  on  the 
wafer). 
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Time  (s) 


Figure  45:  Temperature  responses  of  a  typical  thermocouple,  the  resistance  temperature 
detector,  and  the  recirculating  chiller  bath. 
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Figure  46:  Temperature  reading  obtained  during  calibration  (Thermocouple  201):  (a)  Tem¬ 
perature  obtained  by  RTD  vs.  temperature  reading  from  TC201;  (6)  Error  obtained  when 
substituting  raw  data  into  the  linear  best-fit  equation  for  TC201. 
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5  Design  and  Testing  of  a  Thermodynamic  Filling  Station  for  Miniature  Heat 

Pipes 

5.1  Abstract 

An  experimental  investigation  was  conducted  to  determine  the  feasibility  of  a  thermody¬ 
namic  process  to  fill  miniature  heat  pipes.  A  heat  pipe  filling  station  has  been  designed  and 
constructed,  which  consists  of  a  vacuum  pump,  pressure  gauges,  thermocouples,  a  manifold 
of  valves,  a  large  evacuated  volume,  and  an  environmental  test  chamber.  A  general  filling 
procedure  has  been  formulated.  The  filling  station  has  been  tested  and  preliminary  results 
have  been  obtained  using  a  heat  pipe.  The  target  range  for  the  working  fiuid  fill  amount  for 
the  present  project  is  0.01  to  0.5  g  ±  10%.  Prefiminary  results  presented  herein  show  that 
the  current  process  is  capable  of  filfing  miniature  heat  pipes  with  a  mass  of  approximately 
0.7  g  with  a  repeatability  of  ±  12%.  An  analytical  model  that  includes  the  effect  of  air  within 
the  fill  station  is  presented.  The  analytical  model  overpredicts  the  experimental  results  by 
approximately  26%. 

5.2  Introduction 

Suitably  filling  a  heat  pipe  with  a  working  fluid  is  a  critical  process  in  ensuring  proper 
performance.  The  steps  involved  in  any  filfing  process  include  thoroughly  cleaning  the  heat 
pipe,  evacuating  the  heat  pipe  to  remove  all  non-condensable  gases,  dispensing  the  proper 
amount  of  high  purity  working  fluid,  and  finally,  sealing  the  heat  pipe.  Variations  firom  the 
optimal  charge  can  adversely  affect  the  performance  of  the  heat  pipe  (Castle  et  al.,  2000). 

The  performance  of  the  heat  pipe  is  very  sensitive  to  the  quantity  of  working  fluid  charged. 
A  small  variation  in  the  fill  amount  can  degrade  the  heat  pipe  performance,  possibly  pre¬ 
venting  it  from  functioning.  Undercharging  can  cause  several  effects,  one  of  the  more  serious 
being  premature  dry-out  at  the  evaporator  section  due  to  an  unsaturated  wick  structure. 
This  dry-out  can  cause  an  increase  in  local  temperature,  resulting  in  failure  of  the  com¬ 
ponent  being  cooled  due  to  inadequate  heat  removal  (Gao  et  al.,  2000).  It  is  common  to 
spend  a  considerable  amount  of  time  calculating  the  amount  of  charge  required,  then  delib¬ 
erately  overcharging  by  10-20%  to  prevent  the  possibility  of  undercharging  (Peterson,  1994). 
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Overcharging  the  heat  pipe  will  also  reduce  performance  by  partially  blocking  the  condenser 
section  due  to  the  presence  of  a  puddle  of  liquid. 

There  are  several  methods  for  filling  heat  pipes.  However,  few  of  these  common  practices 
are  feasible  for  the  0.01  to  0.5  g  fill  quantities  required  for  miniature  heat  pipes.  One 
traditional  method  requires  initially  overfilfing  the  heat  pipe,  heating  it  up  to  a  desired 
operating  temperature,  venting  off  non-condensable  gases,  and  then  sealing  the  heat  pipe 
once  isothermal  operation  is  attained.  With  miniature  heat  pipes,  the  fill  quantities  are  too 
small  and  all  the  working  fluid  would  probably  be  lost  in  this  process. 

Another  method  involves  an  evacuation  and  back-filMng  procedure  (Castle  et  al.,  2000). 
The  heat  pipe  is  attached  to  the  heat  pipe  fill  station,  and  both  are  then  evacuated.  The 
working  fluid  is  then  degassed  using  a  freeze-thaw  cycle  or  by  boiling.  Finally,  the  desired 
amount  of  working  fluid  is  drawn  into  the  heat  pipe  via  a  dispensing  burette,  and  the  heat 
pipe  is  sealed.  The  estimated  uncertainties  are  well  over  50%  for  this  process  appfied  to 
the  small  fill  quantities  of  miniature  heat  pipes.  The  fill  range  for  this  process  as  tested  by 
Castle  et  al.  was  1.47  to  4.38  g  ±  5%. 

Gao  et  al.  (2000)  developed  three  different  methods  for  filling  miniature  heat  pipes. 
The  first  process  proposed  was  a  micro  syringe  method  similar  to  the  evacuation/back-filling 
procedure  mentioned  above.  This  method  requires  a  micro  syringe,  a  T-junction,  and  a 
vacuum  pump.  The  syringe  is  filled  with  working  fluid.  The  filhng  needle  is  then  replaced 
with  a  new  dry  needle.  The  syringe  and  heat  pipe  are  connected  to  the  T-junction  and  the 
vacuum  system.  The  heat  pipe  is  then  evacuated.  Lastly,  the  syringe  plunger  is  depressed, 
displacing  all  working  fluid  into  the  heat  pipe.  The  main  drawback  of  this  process  is  that  a 
certain  amount  of  air  inevitably  enters  the  heat  pipe. 

The  second  and  third  methods  are  quite  similar  in  theory  and  process.  The  thermody¬ 
namic  equilibrium  and  capillary  tubing  method  both  include  a  vacuum  system,  working  fluid 
source,  and  either  a  vacuum  chamber  (thermodynamic)  or  capillary  tubing  (capillary).  The 
process  starts  by  evacuating  the  system  and  all  components.  Next,  the  heat  pipe  is  com¬ 
pletely  filled  with  working  fluid  and  then  attached  to  the  filhng  station.  A  valve  is  opened 
and  the  heat  pipe  is  exposed  to  the  chamber/tubing,  which  results  in  the  evaporation  of 
some  of  the  working  fluid.  A  sudden  drop  in  temperature  occurs  due  to  the  evaporation 
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process.  The  pressure  then  rises  and  an  equilibrium  point  is  achieved  with  the  ambient.  The 
amount  of  liquid  evaporated  is  determined  by  the  chamber /tubing  volume.  In  the  capillary 
tubing  method,  the  liquid  removal  is  aided  by  pushing  the  extra  liquid  into  the  capillary 
tubing,  but  not  by  flashing.  Fill  ranges  for  the  chamber  process  were  given  for  two  chambers 
sizes:  Chamber  1  (1850  ml)  0.080  to  0.117  g  (-16.8  to  21.0%);  Chamber  2  (500  ml)  0.132 
to  0.165  g  (-11.1  to  11.2%). 

The  thermodynamic  method  mentioned  by  Gao  et  al.  (2000)  is  very  similar  to  the  process 
chosen  and  implemented  in  the  present  research  and  experimentation.  The  objectives  of  the 
present  experiment  are  to  reduce  the  minimum  fill  amount  beyond  that  given  by  Gao  et  al. 
by  an  order  of  magnitude  (while  maintaining  an  uncertainty  of  approximately  10%),  and  to 
compare  the  results  of  the  present  experiment  with  an  analytical  model  of  the  thermodynamic 
process. 


5.3  Analytical  Model 


The  current  heat  pipe  fill  station  is  shown  in  Fig.  47(o),  which  can  be  divided  into 
three  volumes  for  analysis.  The  vacuum  chamber  volume  extends  to  valves  VI  and  V2, 
and  includes  the  tubing  for  the  dial  pressure  gage,  the  thermocouple  pressure  gage,  and 
the  thermocouple  probe.  The  manifold  volume  includes  tubing  from  valves  V2  and  V3,  the 
critical  orifice,  and  the  tubing  and  fittings  up  to  the  point  at  which  the  heat  pipe  is  attached. 
The  heat  pipe  volume  includes  the  heat  pipe  envelope  and  the  fill  tube.  During  the  filling 
process,  the  heat  pipe  is  charged  with  liquid  and  attached  to  the  fill  station.  The  vacuum 
chamber  is  evacuated,  but  the  manifold  and  heat  pipe  are  at  atmospheric  pressure.  Valve 
V2  is  then  opened,  exposing  the  liquid  in  the  heat  pipe  to  the  vacuum  chamber.  After  a 
portion  of  the  liquid  evaporates,  the  system  reaches  a  steady  state  condition.  At  this  point, 
the  heat  pipe  fiU  tube  is  sealed.  The  objective  of  this  analysis  is  to  determine  the  amount 
of  liquid  in  the  heat  pipe  at  the  final  state  by  accounting  for  the  presence  of  air  within  the 
system. 

The  mass  of  air  initially  within  the  vacuum  chamber  is  given  by  the  perfect  gas  relation. 


P  ,1/ 

_  VCfl  ^  VC 

TTT'a.vc  —  rp^ 


(121) 
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The  initial  pressure  within  the  manifold  is  atmospheric  since  it  is  open  to  ambient  conditions 
when  the  heat  pipe  is  attached  to  the  fill  station.  The  mass  of  air  within  the  manifold  is 


m„ 


PatmVm 

RaTl 


(122) 


In  general,  the  heat  pipe  volume  at  the  initial  state  is  taken  up  by  the  working  fluid  (liquid 
and  vapor)  at  saturated  conditions,  which  is  in  equilibrium  with  any  air  that  may  be  present 
within  the  heat  pipe.  The  mass  of  air  within  the  heat  pipe  is  dependent  upon  the  method 
used  for  introducing  the  initial  liquid  charge.  If  the  Hquid  is  installed  using  a  syringe,  for 
instance,  the  amount  of  air  is  unknown.  Therefore,  the  mass  of  air  originally  in  the  heat 
pipe  will  be  varied  from  zero  to  the  maximum  value  possible  to  determine  the  effect  on  the 
final  conditions  within  the  fill  station.  The  maximum  mass  of  air  possible  within  the  heat 
pipe  can  be  determined  by  assuming  that  no  vapor  is  present 

Patm  {Vhp  —Vl^\) 


fT^a,hp,max  — 


RaTl 


(123) 


where  the  volume  of  the  liquid  is  based  on  the  initial  amount  metered  into  the  heat  pipe. 


Vi,i  =  mi^iVf 


(124) 


The  total  mass  of  air  within  the  system  is 


T  '^a,7n  T  '^a,hp  (1^^) 

After  the  heat  pipe  is  connected  to  the  fill  station,  valve  V2  is  opened  and  the  evaporation 
process  begins.  It  is  assumed  that  the  evaporation  of  the  working  fluid  sweeps  all  of  the  air 
out  of  the  heat  pipe  and  into  the  manifold  and  vacuum  chamber  volumes.  It  is  also  assumed 
that  the  wick  structure  and  the  vertical  orientation  of  the  heat  pipe  prevent  Hquid  droplets 
from  being  ejected  from  the  heat  pipe  when  it  is  mounted  to  the  All  station.  Finally,  it  is 
assumed  that  temperature  for  State  2  is  equal  to  the  initial  temperature  (prior  to  opening 
the  valve  V2).  The  final  pressure  in  the  system  is  the  partial  pressure  of  the  working  fluid 
at  the  ambient  temperature.  The  final  volume  of  air  is 


The  volume  occupied  by  the  working  fluid  is 


K,/.2  =  K  -  14.2 

Using  the  definition  of  quality,  the  mass  of  working  fluid  in  the  vapor  state  is 

V4/,2  -  mi,iVf 


(127) 


mv,2  = 


(128) 


The  mass  of  working  fluid  in  the  Uquid  state  is  determined  using  the  continuity  equation 

'Vwf,2-mi,iVf\ 


mi, 2  =  mil  -  =  mu 

5.4  Experimental  Setup 


(129) 


The  current  filling  station  is  shown  schematically  in  Fig.  47(a)  and  Fig.  48.  The 
station  consists  of  an  aluminum  vacuum  chamber,  a  stainless  steel  manifold  of  valves  and 
interconnecting  tubing,  a  critical  orifice  plate  and  a  mechanical  vacuum  pump.  Initially, 
the  open  heat  pipe  container  is  charged  with  liquid  working  fluid,  which  is  measured  using 
a  precision  scale  (Ohaus  Explorer  Model  E01140).  The  heat  pipe  is  then  attached  to  the 
filling  station  as  shown  in  Fig.  47(a).  The  vacuum  chamber  is  evacuated  and  subsequently 
isolated  from  the  vacuum  pump.  Valve  V2  is  then  opened  to  expose  the  working  fluid  in  the 
heat  pipe  to  the  reduced  pressure  within  the  vacuum  chamber.  This  initiates  evaporation  of 
the  liquid  working  fluid,  which  halts  when  the  system  reaches  a  steady  state  condition  with 
the  ambient.  The  ambient  temperature  is  controlled  by  placing  the  entire  filling  station  into 
an  environmental  chamber  (Thermotron,  AToo  1.0  °C).  Temperatures  within  the  vacuum 
chamber  and  on  the  surface  of  the  heat  pipe  were  measured  using  Type  T  thermocouples 
and  a  data  acquisition  system  (Fluke  2286A  Data  Logging  System).  The  pressure  within 
the  vacuum  chamber  was  tracked  using  a  thermocouple  vacuum  gauge  and  reader  (Varian 
Model  0531  and  Model  804-A)  and  a  dial  absolute  pressure  gauge  (Wallace  and  Tiernan).  A 
mechanical  vacuum  pump  (Duo  Seal  Model  1397)  was  employed  to  evacuate  the  system. 

As  the  preliminary  testing  of  the  system  progressed,  it  was  found  that  a  device  was  needed 
to  control  the  flow  of  working  fluid  from  the  heat  pipe  to  the  vacuum  chamber.  When  the 
valve  V2  is  initially  opened,  a  portion  of  the  liquid  could  flash  off  into  the  vacuum  chamber, 
which  could  result  in  a  slug  of  liquid  being  ejected  from  the  heat  pipe  container.  A  critical 
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orifice  design  was  implemented  to  slow  the  evaporation  process  by  restricting  the  maximum 
mass  flow  rate  of  the  vapor  from  the  heat  pipe.  The  cross-sectional  area  of  the  hole  in  the 
orifice  plate  was  calculated  using  the  following  expression,  which  is  based  on  the  assumption 
that  the  flow  is  choked  at  the  throat. 


Once  the  area  was  known,  the  orifice  plate  was  constructed  using  copper  plate  and  placed  in 
between  two  flanges  (Fig.  49(a)).  Three  orifice  plates  were  made  with  three  different  orifice 
diameters  for  experimentation;  Do  =  0.203  mm  (0.008  in.),  0.381  mm  (0.015  in.)  and  0.762 
mm  (0.0300  in.).  A  photograph  of  the  orifice  plate  installed  in  the  filfing  station  is  shown  in 
Fig.  49(6). 

It  was  necessary  to  accurately  determine  aU  of  the  volumes  of  the  fill  station,  which 
included  all  of  the  components,  the  chamber,  and  the  heat  pipe,  for  the  calculations  made  in 
the  mathematical  model.  AU  volumes  were  found  either  by  calculations  or  by  given  values 
in  the  manufacturer  manuals.  Since  aU  of  the  volumes  were  cylindrical,  the  uncertainty  of 
the  volume  AV  was  found  by  using  the  root-sum-square  method 


AV  = 


(131) 


=  y/ (27rr/iAr)^  -I-  (nr^Ah)'^  (132) 

where  the  uncertainties  of  the  measured  values,  Ar  and  Ah,  were  comprised  of  the  readability 
or  the  uncertainty  of  the  measurement  tool.  The  dial  caliper  used  had  a  readability  of  0.0254 
mm  (0.001  in.),  and  a  readability  per  measurement  of  0.0254  mm  (0.001  in.).  The  scale  used 
had  a  readability  of  0.396  mm  (1/64  in.)  and  a  readability  of  0.793  mm  (2/64  in)  per 
measurement.  The  measured  values  for  aU  of  the  components  and  the  total  volumes  are 
given  in  Table  13. 

5.5  Results  and  Discussion 

Results  of  preliminary  tests  are  shown  in  Fig.  50,  where  the  experimental  data  are  com¬ 
pared  to  the  analytical  model.  While  the  model  significantly  over-predicts  the  experimental 
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data  (approximately  26%),  the  model  matches  the  slope  of  the  linear  best-fit  curve  for  the 
experimental  data  fairly  well.  The  curve  labelled  “Analytical,  No  Air”  shows  the  results 
when  no  air  is  present  in  the  system;  i.e.,  the  vacuum  chamber,  manifold  and  heat  pipe  are 
completely  evacuated.  Under  this  condition,  the  amount  of  liquid  left  in  the  heat  pipe  is  less 
than  with  air  present,  as  expected.  It  was  determined  that  assuming  the  heat  pipe  did  not 
contain  any  air  at  the  start  of  the  process  had  a  neghgible  effect  on  the  results.  The  data 
points  presented  fit  the  finear  curve  to  within  ±12%,  but  more  data  points  are  needed  to 
fully  evaluate  the  repeatability  of  the  system. 

5.6  Conclusions  and  Future  Directions 

Initial  testing  of  the  heat  pipe  filling  station  has  verified  the  viability  of  the  thermody¬ 
namic  process.  Inconsistencies  between  the  analytical  model  and  the  experimental  data  will 
require  further  investigations  into  methods  of  improving  the  agreement.  In  order  to  reduce 
the  amount  of  air  within  the  system,  a  combined  filling  station  is  proposed  as  shown  in  Fig. 
47(6),  where  a  conventional  back-fiUing  station  is  attached  to  the  current  thermodynamic 
filling  station.  The  merits  and  drawbacks  of  this  combined  station  still  need  to  be  evaluated. 
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5.8 


Nomenclature 


A* 

h 

k 


'^a,hp 

f^a,hp,max 

'^a,m 

■ma,t 

'^a,vc 

mi, I 

mi, 2 


m^,2 

m 


P itm 

Pv,2 

Pvc,l 

Po 


To 

Ti 

Vf 

Vg 

y 

Va,2 

Vhp 

^.1 

14 


cross-sectional  area  of  the  critical  orifice  throat, 

cylinder  height,  m 

ratio  of  specific  heats 

mass  of  air  in  the  heat  pipe,  kg 

maximum  possible  mass  of  air  in  the  heat  pipe,  kg 

mass  of  air  in  the  manifold,  kg 

total  mass  of  air  in  the  filling  station,  kg 

mass  of  air  in  the  vacuum  chamber,  kg 

initial  mass  of  liquid  working  fluid  in  the  heat  pipe,  kg 

mass  of  saturated  working  fluid  liquid  in  the  heat  pipe  at  State  2,  kg 

mass  of  saturated  working  fluid  vapor  in  the  filling  station  at  State  2,  kg 

mass  flow  rate,  kg/s 

atmospheric  pressure.  Pa 

partial  pressure  of  the  saturated  working  fluid  vapor  at  State  2,  Pa 
air  pressure  in  the  vacuum  chamber  at  State  1,  Pa 
stagnation  pressure.  Pa 
cylinder  radius,  m 

particular  gas  constant  of  air,  J/(kg-K) 
stagnation  temperature,  K 
initial  system  temperature,  K 
specific  volume  of  saturated  liquid,  m^/kg 
Vg-Vf,  mVkg 

specific  volume  of  saturated  vapor,  m^/kg 
volume,  m^ 

volume  of  air  at  State  2,  m^ 
heat  pipe  volume,  m^ 

initial  volume  of  liquid  working  fluid  in  the  heat  pipe,  m^ 
manifold  volume,  m^ 
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Vt  total  volume  of  the  filling  station  and  heat  pipe, 

Kc  vacuum  chamber  volume,  m^ 

Ku/,2  volume  occupied  by  the  working  fluid  at  State  2,  m^ 
A/i  uncertainty  in  measured  height,  m 

Ar  uncertainty  in  measured  radius,  m 

Al^  uncertainty  in  measured  volume,  m 
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Figure  48:  Thermodynamic  heat,  pipe  filling  station. 
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Wz,l  (g) 


Figure  50:  Mass  of  liquid  in  the  heat  pipe  at  State  2  versus  initial  liquid  charge  {ma,hp  —  0, 
P,,,a  =  30  mTorr,  Ti  =  25°C,  P^tm  =  10"  Pa). 
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6  Acceleration  Testing  of  Three  Raytheon  Heat  Pipes  using  the  AFRL/PRPS 
Centrifuge  Table 

6.1  Abstract 

The  purpose  of  this  experiment  was  to  determine  the  effect  of  adverse  longitudinal  accel¬ 
eration  on  the  operation  of  three  heat  pipes  provided  by  Raytheon.  The  heat  pipe  assembly 
(as  dehvered)  was  mounted  to  a  centrifuge  table  such  that  the  evaporator  sections  were  closer 
to  the  center  of  the  table  than  the  condenser  sections,  which  resulted  in  an  adverse  accelera¬ 
tion  condition.  The  start-up  test  is  described  as  follows:  Rotate  the  centrifuge  table  for  five 
min.  at  the  prescribed  acceleration  level  (5-g  or  10-g),  stop  the  rotation,  wait  for  5  s,  then 
apply  a  fixed  power  input  to  the  evaporator.  The  burst  test  consists  of  the  following:  Apply 
a  fixed  power  input  to  the  evaporator,  allow  the  heat  pipe  to  reach  steady-state  operation, 
rotate  the  centrifuge  table  at  the  prescribed  acceleration  level  (5-g  or  10-g)  for  5  s,  then 
stop  the  rotation.  During  the  start-up  tests,  it  was  found  that  the  3-inch  and  6-inch  heat 
pipes  (3IHP  and  6IHP)  were  not  affected,  but  the  9-inch  heat  pipe  (9IHP)  showed  signs  of 
a  deprime  condition.  The  burst  tests  resulted  in  the  3IHP  and  9IHP  becoming  deprimed  to 
a  significant  degree.  The  response  time  of  the  9IHP  to  the  initial  start-up  was  significantly 
longer  than  the  other  heat  pipes. 

6.2  Experimental  Setup 

The  heat  pipe  experiments  were  carried  out  on  an  eight-foot  diameter  horizontal  cen¬ 
trifuge  table  at  Wright-Patterson  Air  Force  Base  (AFRL/PRPS).  This  table  was  driven  by 
a  20-hp  DC  motor  (General  Electric)  and  was  controlled  remotely  for  safety.  A  data  acqui¬ 
sition  unit  (Keithley  Model  500A)  with  a  personal  computer  running  data  logging  software 
(Viewdac)  was  used  to  collect  temperature  and  acceleration  profiles.  The  heat  pipe  assem¬ 
bly  was  attached  to  the  centrifuge  table  as  shown  in  Pigs.  51  and  52.  An  accelerometer 
(Columbia  Research  Laboratories  Model  SA-307IIPTX)  was  used  to  measure  the  accelera¬ 
tion  at  the  lengthwise  center  of  the  heat  pipes.  Temperature  signals  were  conditioned  and 
amplified  on  the  centrifuge  table.  These  signals  were  transferred  off  the  table  through  the 
instrumentation  shp  ring  assembly,  which  was  completely  separate  from  the  power  slip  ring 
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Name 

r  (in) 

h  (in) 

V  (in^) 

Large  Cylinder 

4.516 

8.031 

514.6 

Cylinder  Weld 

4.516 

0.312 

-1.385 

Hole 

0.233 

0.300 

0.051 

Hole 

0.213 

0.150 

0.021 

Hole 

0.200 

0.285 

0.036 

Hole 

0.218 

0.160 

0.024 

Hole 

0.225 

0.145 

0.023 

Hole 

0.183 

0.170 

0.018 

Hole 

0.168 

0.340 

0.030 

Stainless  Tubing 

0.125 

1.375 

0.068 

Stainless  Tubing 

0.125 

3.25 

0.160 

Dial  Gauge  Tubing 

0.063 

21.0 

0.258 

Male  Extension 

0.093 

1.325 

0.036 

TC  Vacuum  Gauge 

0.132 

2.00 

0.110 

TC  Probe  Fitting 

0.125 

0.250 

0.012 

Reducing  Bushing 

0.140 

0.188 

0.012 

Swagelok 

0.125 

0.50 

0.025 

Swagelok  to  Pipe 

0.125 

1.00 

0.049 

Bellows  Valve 

N/A 

N/A 

0.153 

Stainless  Tubing 

0.125 

1.50 

0.074 

0-ring  Fitting 

0.093 

0.375 

0.010 

Mini  T  Jen 

0.313 

3.00 

0.920 

Mini  T  Jen 

0.313 

1.125 

0.345 

Varian  Valve 

0.188 

1.50 

0.166 

Bellows  Valve 

N/A 

N/A 

0.076 

Heat  Pipe 

N/A 

N/A 

0.047 

AV  (in^)  Percent  Error 


0.013 


0.0005 


0.0002 


0.0004 


0.0003 


0.0003 


0.0002 


0.0004 


0.0019 


0.003 


0.0083 


0.0011 


0.0024 


0.0015 


0.0019 


0.0016 


0.0017 


0.0015 


0.0019 


0.0009 


0.0113 


0.0098 


0.0039 


0.0015 


0.001 


Table  13;  Heat  pipe  filling  station  volume  data  (t4c  =  514.3  in^  =  8.428  x  10  ^  m^,  K 
1.591  in^  =  2.607  x  10“®  =  0.047  in^  =  7.702  x  10"^  m^  K  =  8.455  x  10-^  m^) 


assembly  to  reduce  electronic  noise.  Conditioning  the  temperature  signals  on  the  centrifuge 
table  eliminated  difficulties  associated  with  creating  additional  junctions  within  the  sHp  ring 
assembly. 

Power  was  supplied  to  the  heat  pipe  evaporator  section  by  a  precision  power  supply 
(Kepco  Model  ATE  150-7M)  through  power  sfip  rings  to  the  table.  The  input  power  was 
calculated  using  the  current  and  voltage  readings.  While  the  current  reading  could  be  made 
directly  using  a  precision  ammeter,  the  voltage  across  the  electric  heater  had  to  be  mea¬ 
sured  on  the  rotating  table  because  voltage  drops  were  noted  between  the  control  room  and 
the  table,  regardless  of  the  size  of  wire  used.  Therefore,  the  voltage  at  the  heater  was  ob¬ 
tained  through  the  instrumentation  slip  ring  assembly  and  measured  by  a  precision  voltmeter 
(Hewlett  Packard  Model  3478A). 

6.3  Experimental  Procedure 

Two  separate  procedures  were  developed  to  examine  the  performance  of  the  heat  pipes 
under  adverse  acceleration  conditions.  The  start-up  test  and  the  burst  test  are  discussed 
in  this  section.  In  all  of  the  experiments,  the  table  was  allowed  to  rotate  slowly  to  prevent 
damage  to  the  instrumentation  and  power  slip  rings.  This  baseline  rotational  speed  was  cu  = 
0.0329  rad/s,  which  imposed  a  radial  acceleration  of  Or  =  1.24  x  10~^.-g.  While  a  baseplate 
heater  was  present  on  the  heat  pipe  assembly,  it  was  not  used  during  these  tests.  In  addition 
to  recording  the  temperatures  and  accelerations  to  computer  data  files  continuously  during 
the  test,  the  temperature,  acceleration,  and  heater  power  information  (voltage  and  current) 
were  recorded  into  a  logbook  at  intervals  during  the  testing.  For  all  tests,  the  electrical 
power  input  to  the  evaporator  heater  was  approximately  17  W.  To  prevent  damage  to  the 
heater  mounted  to  the  heat  pipe  evaporator  section,  a  maximum  evaporator  temperature 
of  115°C  was  imposed.  Each  test  was  run  three  times  to  evaluate  the  repeatability  of  the 
experimental  procedure. 

Burst  test  procedure 

1.  Turn  on  the  data  logging  software  to  record  basehne  readings. 

2.  Turn  on  the  power  to  the  heater. 
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3.  Allow  the  heat  pipe  to  reach  steady-state  operation. 

4.  After  about  25  min.  of  steady-state  operation,  bring  the  table  rapidly  to  the  desired 
radial  acceleration  by  increasing  the  rotational  speed. 

5.  Hold  this  acceleration  for  5  s. 

6.  Rapidly  decelerate  the  table  back  to  the  baseline  rotational  speed. 

7.  Monitor  the  temperatures  in  case  of  evaporator  dry-out. 

8.  If  the  pipe  returns  to  steady-state  operation,  record  data  for  about  20  min. 

9.  Repeat  steps  4  through  7,  using  the  same  acceleration. 

10.  If  the  heat  pipe  returns  to  steady-state  operation  after  the  second  burst,  record  data 
for  about  20  min. 

11.  Shut  off  the  heater  power. 

12.  Record  the  cool  down  process  for  several  minutes. 

13.  Turn  off  the  data  acquisition  software  and  download  the  data  files. 

Start-up  test  procedure 

1.  Start  the  data  logging  software  to  record  basehne  data. 

2.  Increase  the  rotational  speed  of  the  table  to  reach  the  desired  radial  acceleration.  Allow 
the  table  to  rotate  at  this  speed  for  5  min. 

3.  Bring  the  table  rapidly  back  down  to  its  basehne  speed. 

4.  Wait  5  s  and  then  turn  on  the  evaporator  heater. 

5.  Allow  the  heat  pipe  to  reach  steady-state  operation. 

6.  Run  the  heat  pipe  at  steady-state  condition  for  about  30  min. 

7.  Turn  off  the  power  to  the  heater. 
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8.  Record  the  cool  down  process  for  several  minutes. 


9.  Stop  the  data  logging  software  and  download  the  data  files. 

For  each  test,  plots  were  generated  of  temperature  vs.  time  and  radial  acceleration  vs. 
time. 

6.4  Results  and  Discussion 

Three  heat  pipes  have  been  tested  to  determine  their  thermal  performance  characteristics 
when  subjected  to  adverse  acceleration  conditions.  The  heat  pipe  assembly  was  mounted  to 
a  centrifuge  table  at  AFRL/PRPS  to  impose  a  radial  acceleration  along  the  axial  direction 
of  the  heat  pipes.  Included  in  this  section  are  typical  temperature-time  traces,  steady-state 
measurements  of  the  temperature  difference  between  the  evaporator  and  middle  (adiabatic) 
thermocouples,  and  initial  start-up  times,  both  for  stationary  start-ups  and  start-ups  follow¬ 
ing  long  periods  of  acceleration. 

Figure  53  presents  typical  temperature-time  traces  during  a  5-g  start-up  test.  All  three 
heat  pipes  reached  the  steady  state  in  a  relatively  short  time.  This  steady-state  condition 
was  maintained  for  approximately  30  min.,  and  then  the  heater  power  was  turned  off.  For 
the  3IHP  and  6IHP,  the  temperature  difference  between  the  evaporator  and  adiabatic  sec¬ 
tions  (ATea)  was  lower  than  that  for  the  9IHP.  This  temperature  difference  was  chosen  as 
an  indicator  of  wick  priming  based  on  previous  experience  with  heat  pipes  ([l]-[8]).  The 
evaporator-adiabatic  temperature  differences  for  all  of  the  start-up  tests  are  shown  in  Fig. 
54,  where  ATea  was  significantly  higher  and  more  scattered  for  the  9IHP.  Both  of  these  points 
indicate  that  the  wick  was  poorly  primed  prior  to  start-up  for  the  9IHP. 

Figure  55  shows  the  results  for  four  burst  tests.  In  Fig.  55  (o),  the  3IHP  was  allowed  to 
reach  a  steady-state  condition  while  the  centrifuge  table  was  rotating  at  its  baseline  speed. 
At  approximately  t  =  2500  s,  the  rotational  speed  was  increased  such  that  =  5.0-g  for 
5  s,  and  then  the  speed  was  reduced  to  the  baseline.  While  the  adiabatic  and  condenser 
temperatures  remained  fairly  constant,  the  evaporator  temperature  increased  significantly, 
indicating  a  partial  dry-out  event.  The  heat  pipe  was  again  allowed  to  reach  another  steady- 
state  condition,  and  another  burst  test  was  carried  out  at  approximately  t  =  4500  s.  The 
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evaporator  temperature  increased  further,  and  the  power  to  the  heater  was  shut  off  due  to 
a  run-away  dry-out  situation. 

Figures  55(6)  and  55(c)  compare  two  burst  tests  for  the  6IHP.  Of  particular  interest  is 
the  initial  start-up  period  shown  in  Fig.  55(6),  where  a  wick  priming  event  occurred  at 
approximately  t  —  600  s.  In  this  event,  the  temperature  difference  between  the  evaporator 
and  adiabatic  sections  was  initially  high  {AT^a  —  23°C).  After  the  wick  primed,  the  adiabatic 
temperature  increased  dramatically,  which  decreased  the  temperature  difference  to  ATea  — 
11°C.  Two  burst  tests  were  carried  out  in  Fig.  55(6)  at  t  =  2500  and  4000  s,  and  the  6IHP 
did  not  deprime  in  either  case.  In  Fig.  55(c),  the  wick  did  not  prime  during  start-up,  and  the 
temperature  difference  remained  at  AT^a  —  22°C  prior  to  the  acceleration  burst.  Again,  two 
burst  tests  were  performed  (t  =  2500  and  4500  s).  In  each  of  these  tests,  the  temperature 
difference  actually  decreased  sUghtly,  showing  that  the  heat  pipe  wick  structure  was  partially 
primed  by  the  acceleration  burst. 

Figure  55(d)  shows  the  performance  of  the  9IHP  during  a  burst  test  with  bursts  at 
t  =  2500  and  4000  s.  Initially,  the  heat  pipe  was  primed  with  a  very  low  temperature 
difference  {ATea  —  7°C).  After  the  first  burst  test,  the  adiabatic  temperature  decreased 
significantly  which  resulted  in  ATea  —  25°  C.  The  second  burst  test  resulted  in  a  run-away 
dry-out  situation,  so  the  test  was  halted. 

Figure  56  shows  the  steady-state  data  for  all  of  the  burst  tests,  except  for  the  run-away 
dry-out  situations,  since  steady-state  data  could  not  be  collected.  This  test  had  a  detrimental 
impact  on  the  wick  priming  for  the  3IHP  and  the  9IHP,  where  ATea  increased  after  the  burst 
test.  The  6IHP  was  less  sensitive  to  the  burst-test  conditions  in  comparison  to  the  other 
heat  pipes. 

The  response  times  of  the  heat  pipes  to  the  initial  application  of  the  heat  load  are 
displayed  in  Fig.  57.  These  times  were  calculated  using  the  63%  time-constant  method, 
where  the  time  constant  is  defined  as  the  time  required  for  the  evaporator  temperature  to 
reach  63%  of  its  steady-state  value.  These  data  include  the  initial  start-up  responses  for  both 
the  start-up  tests  and  the  burst  tests,  where  the  power  to  the  evaporator  is  first  energized. 
The  responses  due  to  the  actual  burst  tests  (immediately  after  the  radial  acceleration  was 
increased  and  then  decreased)  were  not  included  due  to  difficulties  in  defining  an  appropriate 
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time  constant.  The  time  constants  for  all  three  heat  pipes  were  very  repeatable,  with  the 
time  constant  for  the  9IHP  being  more  than  twice  as  long  as  those  of  the  other  heat  pipes. 
It  is  believed  that  the  3IHP  was  not  properly  charged  based  on  the  following  observations. 

1.  During  the  start-up  tests,  ATga  was  higher  for  the  3IHP  than  the  6IHP  (Fig.  54). 

2.  For  the  burst  tests,  AT^a  for  the  3IHP  was  greater  than  both  the  6IHP  and  the  9IHP 
(Fig.  56). 

3.  The  response  time  for  the  3IHP  (52  s)  was  higher  than  the  6IHP  (44  s)  (Fig.  57). 

It  was  expected  that  the  3IHP  would  perform  with  a  lower  AT^a  and  with  a  shorter  response 
time  than  the  6IHP  due  to  the  shorter  effective  heat  pipe  length  of  the  3IHP. 

6.5  Conclusions 

Based  on  the  experiments  described  above,  the  following  conclusions  have  been  made 
concerning  the  performance  of  the  three  heat  pipes  under  adverse  longitudinal  acceleration 
conditions: 

1.  During  the  start-up  tests,  the  9IHP  showed  signs  of  a  deprimed  wick  structure. 

2.  The  burst  tests  caused  the  3IHP  and  the  9IHP  to  become  deprimed. 

3.  The  start-up  response  time  for  the  9IHP  was  significantly  higher  than  the  other  heat 
pipes  (more  than  double). 

4.  The  3IHP  may  not  have  been  properly  charged  with  working  fluid. 
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Figiire  51:  Close-up  view  of  the  heat  pipe  assembly  mounted  to  the  centrifuge  t,able. 
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Q3  13 


Centrifuge  Table 
(Seen  from  above) 


Heat  Pipe 
Assembly 


Accelerometer 


[•'igui’c  52:  Placement  of  the  heat  pipe  assembly  on  the  centrifuge  table  (7?  —  44.5  inches). 
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Figure  53:  Typical  start-up  test  results  for  ttr  =  5.0-g:  (a)  3IHP  (RAY33);  (6)  6IHP  (RAY22); 
(c)  9IHP  (RAY9). 


Figure  54:  Temperature  differences  between  the  evaporator  and  adiabatic  sections  for  the 
start-up  tests. 
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Figure  55:  Typical  burst  test  results  for  Ur  =  5.0-g:  (a)  3IHP  (RAY28);  (6)  6IHP  (RAY16); 
(c)  6fflP  (RAY17);  (d)  9IHP  (RAY2). 
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Figure  56:  Temperature  differences  between  the  evaporator  and  adiabatic  sections  for  the 
burst  tests. 
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Heat  Pipe  Length  (inches) 


Figure  57:  Start-up  response  characteristics. 
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Table  14:  Temperature  differences  between  the  evaporator  and  adiabatic  sections  for  the 
start-up  tests  (Fig.  54). 


Heat  Pipe  Length 
(inches) 

(°c) 

ttr  =  5.0-g 

Or  =  10.0-g 

3 

14.4 

13.7 

3 

15.2 

13.7 

3 

13.9 

14.8 

6 

9.6 

9.8 

6 

10.0 

9.2 

6 

9.6 

10.2 

9 

22.7 

20.8 

9 

23.1 

24.7 

9 

24.5 

27.0 

Table  16:  Start-up  response  characteristics  (Fig.  57). 


Heat  Pipe  Length 
(inches) 

Start-Up  Time  Constant  (s) 

Stationary 

ttr  =  5.0-g 

3 

X 

50 

X 

3 

52 

52 

X 

3 

50 

52 

X 

3 

56 

X 

54 

3 

48 

X 

54 

3 

46 

X 

56 

6 

X 

44 

X 

6 

44 

44 

X 

6 

40 

44 

X 

6 

42 

X 

42 

6 

X 

44 

6 

42 

X 

44 

9 

X 

122 

X 

9 

138 

120 

X 

9 

124 

128 

X 

9 

128 

X 

134 

9 

118 

X 

118 

9 

118 

X 

128 
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